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Preface 


Beginning with Volume XX, the Deep Space Network Progress Report changed from 
the Technical Report 32- series to the Progress Report 42- series. The volume number 
continues the sequence of the preceding issues. Thus, Progress Report 42-20 is the 
twentieth volume of the Deep Space Network series, and is an uninterrupted follow-on to 
Technical Report 32-1526, Volume XIX. 

This report presents DSN progress in flight project support, tracking and data 
acquisition (TDA) research and technology, network engineering, hardware and software 
implementation, and operations. Each issue presents material in some, but not all, of the 
following categories in the order indicated. 

Description of the DSN 

Mission Support 

Ongoing Planetary /Interplanetary Right Projects 
Advanced Right Projects 

Radio Astronomy 

Special Projects 

Supporting Research and Technology 
Tracking and Ground-Based Navigation 
Communications-Spacecraft /Ground 
Station Control and Operations Technology 
Network Control and Data Processing 

Network and Facility Engineering and Implementation 
Network 

Network Operations Control Center 
Ground Communications 
Deep Space Stations 
Quality Assurance 

Operations 

Network Operations 
Network Operations Control Center 
Ground Communications 
Deep Space Stations 

Program Planning 
TDA Planning 

In each issue, the part entitled '^Description of the DSN'* describes the functions and 
facilities of the DSN and may report the current configuration of one of the seven DSN 
systems (Tracking, Telemetry, Command, Monitor & Control, Test & Training, Radio 
Science, and Very Long Baseline Interferometry). 

The work described in this report series is either performed or managed by the 
Tracking and Data Acquisition organization of JPL for NASA. 
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N79«16001 


Network Functions and Facilities 

N. A. Renzetti 

Office of T racking and Data Acquisition 


The obje> Hves, functions, and organization of the Deep Space Network are 
summarized; deep space station, ground communication, and network operations control 
capabilities are described. 


The Deep Space Network was (established by the National 
Aeronautics and Space Administration (NASA) Office of 
Space Tracking and Data Systems and is under tl!e system 
management and teclinical direction of the Jet Propulsion 
Laboratory (JPL). The network is designed for two-way 
communications with unmanned spacecraft traveling approxi- 
mately 16,000 km (10,000 miles) from Earth to the farthest 
planets and to the edge of our solar system. It has provided 
tracking and data acquisition support for the following NASA 
deep space exploration projects: Ranger, Surveyor, Mariner 
Venus 1962, Mariner Mars 1964, Mariner Venus 1967. Mariner 
Mars 1969, Mariner Mars 1971, and Mariner Venus-Mercury 
1973, for which JPL has been responsible for the project 
management, the development of the spacecraft, and the 
conduct of mission operations; Lunar Orbiter, for which the 
Langley Research Center carried out the project management, 
spacecraft development, and conduct of mission operations; 
Pioneer, for which Ames Research Center carried out the 
project management, spacecraft development, and conduct of 
mission operations; and Apollo, for which the Lyndon B. 
Johnson Space Center was the project center and the Deep 
Space Network supplemented the Manned Space Flight Net- 
work, which was managed by the Goddard Space Flight 
Center. The network is currently providing tracking and data 
acquisition support for Helios, a joint U.S./West German 
project; Viking, for which Langley Research Center provided 
the Lander spacecraft and project management until May, 


1978, at which time project management and mission opera- 
tions were transferred to JPL, and for which JPL provided the 
Orbiter spacecraft; Voyager, for which JPL provides project 
management, spacecraft development, and is conducting 
mission operations; and Pioneers, for which the Ames 
Research Center provides project management, spacecraft 
development, and conduct of mission operations. The network 
is adding new capability to meet the requirements of the 
Galileo mission to Jupiter, for which JPL is providing the 
Orbiter spacecraft, and the Ames Research Center the 
probe. In addition, JPL will carry out the project management 
and the conduct of mission operations. 

The Deep Space Network (DSN) is one of two NASA 
networks. The other, the Spaceflight Tracking and Data 
Network (STDN), is under the system management and 
technical direction of the Goddard Space Flight Center 
(GSFC). Its function is to support manned and unmanned 
Earth-orbiting satellites. The Deep Space Network supports 
lunar, planetary, and interplanetary flight projects. 

From its inception, NASA has had the objective of 
conducting scientific investigations throughout the solar sys- 
tem. It was recognized that in order to meet this objective, 
significant supporting research and advanced technology devel- 
opment must be conducted in order to provide deep space 
telecommunications for science data return in a cost effective 



manner. Therefore, the Network is continually evolved lo keep 
pace with the state of the art of telecommunications and data 
handling. It was also recognized early that close coordination 
would be needed between the requirements of the flight 
projects for data return and the capabilities needed in the 
Network. This close collaboration was effected by the appoint- 
ment of a 1 racking and Data Systems Manager as part of the 
flight project team from the initiation of the project to the 
end of the mission. By this process, requirements were 
identified ea’-ly enough to provide funding and implementa- 
tion in time for use by the flight project in its flight phase. 

As of July 1972, NASA undertook a change in the interface 
between the Network and the flight projects. Prior to that 
time, since 1 January 1964, in addition to consisting of the 
Deep Space Stations and the Ground Communications 
Facility, the Network had also included the mission control 
and computing facilities and provided the equipment in the 
mission support areas for the conduct of mission operations. 
The latter facilities were housed in a building at JPL known as 
the Space Flight Operations Facility (SFOF). The interface 
change was to accommodate a hardware interface between the 
support o( the network operations control functions and those 
of the mission control and computing functions. This resulted 
in the flight projects assuming the cognizance of the large 
general-purpose digital computers which were used for both 
network processing and mission data processing. They also 
assumed cognizance of all of the equipment in the flight 
operations facility for display and communications necessary 
for the conduct of mission operations. The Network then 
undertook the development of hardware and computer soft- 
ware necessary to do its network operations control and 
monitor functions in separate computers. A characteristic of 
the new interface is that the Network provides direct data flow 
to and from the static rs; namely, metric data, science and 
engineering telemetry, atid such network monitor data as are 
useful to the flight project. This is done via appropriate ground 
communication equipment to mission operations centers, 
wherever they may be. 

The principal deliverables to the users of the Network are 
carried out by data system configurations as follows; 

(1) The DSN Tracking System generates radio metric data, 
i.e., angles, one- and two-way doppler and range, and 
transmits raw data to Mission Control. 

(2) The DSN Telemetry System receives, decodes, records, 
and retransmits engineering and scientific data gener- 
ated in the spacecraft to Mission Control. 

(3) The DSN Command System accepts spacecraft com- 
mands from Mission Control and transmits the com- 
mands via the Ground Communication Facility to a 


Deep Space Station. The commands are then radiated 
to the spacecraft in order to initiate spacecraft func- 
tions in flight. 

(4) The DSN Radio Science System generates radio science 
data, i.e., the frequency and amplitude of spacecraft 
transmitted signals affected by passage through media 
such as the solar corona, planetary atmospheres, and 
plantetary rings, and transmits this data to Mission 
Control. 

(5) The DSN Very Long Baseline Interferometry System 
generates time and frequency data to synchronize the 
clocks among the three Deep Space Communications 
complexes. It will generate universal time and polar 
motion and relative Deep Space Station locations as 
by-products of the primary data delivery function. 

The data system configurations supporting testing, training, 
and network operations control functions are as follows: 

(1) The DSN Monitor and Control System instruments, 
transmits, records, and displays those parameters of the 
DSN necessary to verify configuration and validate the 
Network. It provides the tools necessary for Network 
Operations personnel to control and monitor the Net- 
work and interface wiih flight project mission control 
personnel. 

(2) The DSN Test and Training System generates and 
controls simulated data to support development, test, 
training and fault isolation within the DSN. It partici- 
pates in mission simulation with flight projects. 

The capabilities needed to carry out the above functions 
have evolved in the following technical areas: 

(1) The Deep Space Stations, which are distributed around 
Earth and v'hich, prior to 1964, formed part of the 
Deep Space Instrumentation Facility. The technology 
involved in equipping these stations is strongly related 
to the state of the art of telecommunications and 
flight-ground design considerations, and is almost com- 
pletely multimission in character. 

(2) The Ground Communications Facility provides the 
capability required ir., the transmission, reception, and 
monitoring of Earth-based, point-to-point communica- 
tions between the stations and the Network Operations 
Control Center at JPL., Pasadena, and to the JPL Mis- 
sion Operations Centers. Four communications dis- 
ciplines are provided: teletype, voice, high-speed, and 
wideband. The Ground Communications Facility uses 
the capabilities provided by common carriers through- 
out the world, engineered into an integrated system by 
Goddard Space Flight Center, and controlled from the 
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communications Center located in the Space Flight 
Operations Facility (Building 230) at JPL. 

The Network Operations Control Center is the functional 
entity for centralized operational control of the Network and 
interfaces with the users. It has two separable functional 
elements; namely. Network Operations Control and Network 
Data Processing. The functions of the Network Operations 
Control are: 


> (1) Control and coordination of Network support to meet 

i commitments to Network users. 


i, (2) Utilization of the Network data processing computing 

I capability to generate all standards and limits required 

fe for Network operations. 

^ (3) Utilization of Network data processing computing 

I: capability to analyze and validate the performance of 

i all Network systems. 


The personnel who carry out the above functions are located 
in the Space Flight Operations Facility, where mission opera- 


tions functions arc carried out by certain flight projects. Net- 
work persormel are directed by an Operations Control Chief. 
The functions of the Network Data Processing arc: 

(1) Processing of data used by Network Operations Control 
for control and analysis of the Network. 

(2) Display in the Network Operations Control Area of 
data processed in the Network Data Processing Area. 

(3) Interface with conununications circuits for input to 
and output from the Network Data Processing Area. 

(4) Data logging and production of the intermediate data 
records. 

The personnel who cany out these functions are located 
approximately 200 meters from the Space Flight Operations 
Facility. The equipment consists of minicomputers for real- 
time data system monitoring, two XDS Sigma 5s, display, 
magnetic tape recorders, and appropriate interface equipment 
with the ground data communications. 
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DSN Monitor and Control System, Mark 111-78 

W. G. Stinnett 

TDA Engineering 


A description of the DSN Monitor and Control System, Mark III- 78, is discussed. The 
major implementation required to evolve from the Mark IIT75 to the Mark IIT78 
configuration is also discussed. 


I. Introduction 

The DSN Monitor and Control System, as well as the other 
DSN systems, has undeigone a significant evolution from late 
1976 through the major portion of 1978. This evolution has 
resulted in the Mark 111-78 configuration. Major implementa- 
tion occurred in each of the three DSN facilities: Deep Sp :;- 1 
Stations (DSS), Ground Communications Facility (GCF), a:ic’ 
the Network Ciperations Control Center (NOCC). At the DSS, 
implementation of new computers for the Telemetry, Com- 
mand, Tracking, and Radio Science Systems significantly 
changed the DSN Monitor and Control System. In the GCF, 
implementation has occurred at tlie stations ajiu at JPL in tlie 
central communications terminal. A major software upgrade 
occurred within the NOCC due to the changes at the stations 
and the Ground Communications Facility. 

The DSN Monitor and Control System, Mark III-78 
configuration is multimission and is utilized by the DSN to 
support the Pioneer 6-11, Helios 1 and 2, Viking, Voyage , and 
Pioneer Venus-1978 missions. Figure 1 provides a block 
diagram of the DSN Moidtor and Control System, Mark 111-78 
configuration. 

II. Dmp space station 

As part of tlie DSN Monitor and Control System, the DSS 
Monitor and Control Subsystem (DMC) provides two primary 


functions: (l)the DMC provides the central control and 
morutor capability for the stations, and (2) the DMC receives 
control data from and furnishes monitor data to the N(XC for 
Network Central Control and Monitor purposes. 

There are three major assemblies of the DMC that accom- 
plish the above two functions. The three major assemblies are: 

(1) An XDS 920 computer. 

(2) A station monitor and control console. 

(3) A special input-output terminal. 


A. XDS 920 Computer 

The XDS 920 computer has been a part of the DMe for 
many years, hs designation is the Digital Instrumentation 
Subsystem (DIS). It has been utilized as a central collection 
point for station monitor data and interfacing with the NOCC. 
For collecting station monitor data, the DIS has interfaces 
ivith all of the station subsystems. The station monitor data 
are accumulated by the DIS and displayed to the station 
operator located at the station monitor a^id control console. 
Additionally, a subset of this station monitor data is sent to 
the NOCC for Network central monitoring purposes. A 
significant implementation effort was undertaken for the DIS 
to convert from the Mark 1II-7S to the Mark 111-78 system. 


i 
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Major interface changes occurred with the implementation 
of new computers in the DSS Command, Telemetry, Tracking, 
and Radio Science Subsystems. T!ie interface with these 
subtystem computers is via the DSN standard computcr-to- 
computer intercommunication device, a star switch controller. 
Significant implementation was required for the DIS to 
operate in this computer network. The new generation 
computers are approximately 10 times as fast for memory 
cycle time and are 16-bit word machines versus 24-bit word 
for the DiS (XDS 920). A special hardware interface device 
was installed between the DIS and the star switch controller to 
convert from 24-bit to i 6-bit words, and nee versa, and to 
speed up access time for interfacing ;vith the star switch 
controller. To accommodate interfaces with the new com- 
puters and new piotocol and timing requirements with the star 
switch controllers, significant software changes were required 
in the DIS. Additionally, the stmctaie of the monitor data 
parameters was either new or significantly modified from that 
v^ich previously existed. 

Another interface change, from tlK Mark HI-75 to Mark 
iII-78 configuraticn, was implemented for the DIS. The 
interface with the GCF high-speed daia circuit, for receiving 
and sending data from and to the NOCC, changed. The DIS to 
GCF interface for the Mark 111-78 sysrem is vLi the star switch 
controUer. The GCF has implemented a computer, the 
Communication Monitor Formatter (CMF), at each DSS. Tlie 
DIS generates high-speed data blocks, consisting of station 
monitor data, and proviaes rhe data to the CMF via the sta*- 
switch controUer. The data are i^.ultiplexed with other systems 
data (e.g., telemetry) onto vhe high-speed data line for 
transmission to f 3CC. TTie DIS receives control data from the 
NOCC via this same interface. 

B. Station Monitor and Control Consok 

The station monitor and control (SMC) console has been 
provided as part of the Mark 111-78 system configuration. The 
SMC console provides a central coi.trol and monitor operator's 
position for some of the subsystems at the DSS (other 
subsystems are controlled via another input-output device - 
see Subsection C, below) The SMC consoles are currently in 
limited use in the Network; however, work is in process to 
make them a more useful operational tool. Additional station 
subsystems are being modified to permit central station 
control from the SMC consoles. 

C. Special Input-Output Tetminal 

With the implementation of new computers a the Com- 
mand, Telemetry, Tracking, and Radio Science Subsystems 
(and the CMF at the DSS), a special terminal was installed as a 
central DSS input-oatput device. It consists of four 'elements. 


two kcyboard-CRTs, a hardcopy printer, and “host” software 
capability. The keyboard-CRTs and printers interface with the 
host software (\^ch runs in any of the new compiitere), and 
the host software then interfaces with all of the nev/ 
computers via tlie star switch controUer. In this manner, the 
new computers can be centraUy controlied and monitored. 


III. Ground Communlcaitions Facility 

Significant implementation in the GCF has occurred over 
the last two years. Concurrent with these changes, the GCF 
Monitor and Control Subsystem (GMC) has undergone major 
changes. At J PL in the GCF central communications terminal, 
the Central Communications Monitor (CCM) computer has 
keen implemented to provide central facility monitor capabU- 
ity. The CCM receives GCF monitor data from two sources: 
(1) from each Deep Space Station CMF, and (2) from equip- 
ment located at the central communication terminal at JPL. 

At the DSS, each CMF monitors the status of die GCF 
equipment and provides periodic status messages, via high- 
speed data circuit, to the central communications tenninal 
CCM. It should be noted that the control of the CMF is from 
the central input-output termin?! at the DSS (see Subscc 
tion II-C, above). 

Within the central communications terminal, the CCM 
monitors the status of the GCF equipment. Major implementa- 
tion has occurred within the central communications terminal 
over the past two years. Computers, as part of the GCF 
liigh-speed subsystem, have been implemented :o allow mes- 
sage switching and error detection-correction. These functions 
are being accomplished in the high-speed switch (HSW) and 
the error detection-correction (EDC) assemblies. Tlie new 
computers interface to the CCM via the st:jidard DSN 
computei interface — the star switch controller. 


IV. Network Operations Control Center 

The NOCC Monitor and Control Subsystem (NMC) pro- 
vides the processing capability for central network monitor 
and control. Although no major hardware changes have 
occurred in the NMC, significant software changes were 
required due to major modifications in the DSS and GCF. The 
high-speed data messages received from the DSS underwent 
significant modification due to the implementation that 
occurred there. Still to be implemented in the NMC is tlie 
ability to process GCF monitor data from the CCM. When this 
has been accomplished, the DSN Monitor and Control System, 
Mark III-78 configuration will be complete. 
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Viking Continuation Mission Support 

R. L. Gillette 

Deep Space Network Operations Section 


This article covers the period from I July through 31 August 1978, It reports on DSN 
support of Viking spacecraft activities during the period and continues reporting on the 
DSN Viking Command and Tracking support 


I. Viking Operations 

A. Status 

The Viking Orbiter 2 (VO-2) mission ended during this 
reporting period after a severe gas leak depleted its remaining 
supply of attitude stabilizing gas. On 20 July, the VO-2 
spacecraft was brought out of its roU-drift mode of operation 
to collect its last sequence of scientific observations of Mars. 
On 24 July the VO-2 spacecraft developed a severe gas leak in 
the positive roll valve, and on 25 July drifted off sun 
acquisition. 

The Mars mission officially ended for VO-2 at 11:02 p.m. 
on 25 July when radio signals from the spacecraft cea. L 
Launched from Earth on 20 August 1975, VO-2 will now 
slowly tumble about Mars for at least 50 years before crashing 
to the surface. 

The Viking Orbiter 1 (VO-1) spacecraft continued to 
operate normally during this reporting period as it collected 
and returned to Earth weather data and Mars photos as well as 
relaying to Earth data from the two Viking Landers. 

The Viking Landers continue to operate in an automatic 
mode of operation as they collect and relay to Earth, via 
VO-1, Mars weather information and photos. 


B. Radio Science 

Radio Science activities and experimentation continued 
during July and August. These activities include rear- 
simultaneous Lander /Orbit ranging, VO-1 Earth occuKation 
coverage and the Gravity Field Experiment. 

II. Network Support 

Table 1 shows the DSN tracking support for January 
through August 1978. A significant decrease in tracking 
support is notable in August. This is primarily due to the 
increase in tracking support required from higher prionty 
protects such as the Voyager and Pioneer Venus Missions. The 
Viking Project will continue to have major conflicts with other 
projects for tracking support during the Continuation Mission. 

Table 2 gives the total number of commands transmitted by 
the DSN for Viking Project during January through August 
1978. The reduction in the number of commands beginning in 
June reflects the reduction in Viking operation activity 
anticipated during the Continuation Mission. 

Figure 1 shows the performance of the DSN delivery of the 
telemetry Intermediate Data Record (IDR) during January 
1977 through August 1978. The IDR is a magnetic tape, 
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produced by the DSN, of digital telemetry data. As a goal, 
within existing time constraints, the DSN attempts to provide 
100 percent of the data recorded at the Deep Space Station on 
the Digital Original Data Record (DODR) magnetic tape. 
IDR’s arc required to be delivered within 24 ho following 
the end of a station pass. 

During July and August 1978, the IDR delivery time 
significantly improved from the previous months. The delivery 
time dropped from 23 hours to under 7 hours. 


In May of 1978 a new position was created titled DSN Data 
Records Engineer. The primary function of the Data Records 
Engineer is to coordinate the data record (IDR) production 
conflicts among the many DSN projects, including Viking. The 
improvement in IDR delivery times can be attiibuted to the 
coordination efforts of the DSN Data Records Engineer. 

The percent of telemetry data delivered to the project also 
increased during July and August with an average of 99.1 
percent of telemetry data delivered- 
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Tabtol. DSN tracking wipport 197$ 


DSS 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

1 

Aug 




19 


18 

1 

2 


11 


- 


_ 




_ 




167 


164 

1 

13 




4 

5 

1 





12 

- 




_ 

_ 

_ 

_ 



13 

22 

1 

8 






40 

25 

37 

47 

25 ! 


14 

13 

14 



1 







319 

218 

272 

365 

231 


89 

85 


22 

23 



2 


2 

1 

42 






_ 




166 

76 

201 

192 

4 


1 9 

3 


57 

14 

15 

27 

40 

46 

51 

23 

43 










294 

207 

68 

163 

281 

346 

449 

135 



11 

6 


11 

3 


3 

44 

- 



- 



_ 



1 

51 

43 


77 

17 


11 



38 

25 

15 i 

1 2 



, 

61 






_ 

_ 



88 

243 

281 

159 

26 



5 


61 

2 

3 

3 

8 

3 

1 

1 

62 










27 

25 

22 

29 

75 

33 

7 

9 


13 

21 

30 

42 

26 

41 

40 

16 

63 










101 

181 

293 

440 

206 

343 

291 

130 


156 

138 



132 

128 

no 

58 

Total 










1067 

1014 



1064 

982 

858 

378 


^Number of tracks the summation of all Viking spacecraft tracked, 
^rack time: scheduled station support in hours. 
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IflDIv Z« liUlllDQr Of COnMIWnOO VWMIlllllM In VW VIKinQ mlMIOn 

during January through Auguat ig78 


DSS 

JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

11 

1947 



119 



1628 

6 

0 


12 

- 

1 

1 

374 

- 

_ 


— 

14 

4565 

1079 

1326 

3032 

2845 

1214 

870 

562 

42 

1447 

1305 

261 

1079 

0 

0 

0 

0 

43 

1593 

1732 

124 

275 

1326 

1567 

791 

7 

44 


255 

3 

- 

373 

0 

_ 

30 

61 

992 

3548 

1073 

1503 

0 

0 

— 

0 

62 

1 

1006 

461 

73 

0 

960 

0 

683 

63 

895 

128 

2597 

4005 

669 

2451 

2984 

2073 

TOTAL 

11440 

9054 

5965 

10341 

6841 

6198 

4645 

3355 
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Pioneor Venus Mission Support 

A. C. Bouck 
0$N Operations Section 


This article reports on activities within the Deep Space Network (DSN) to prepare 
for the Pioneer Venus Multiprobe Venusian Encounter. 


I. Mission Description 

The Pioneer Venus Multiprobe Mission is designed to inves- 
tigate the Venusian atmosphere down to the lowest scale 
height above the surl'jce (probe survival after impact is not a 
mission requirement). The mission consists of a Large Probe, 
thiee Small Probes and the Bus. The Bus is the mother space- 
craft that carried and supported the probes from launch until 
probe separation. 

The four probes are mdividual spacecraft and transmit 
their data directly to Earth. The Large Probe carries seven 
scientific instruments, a radio receiver for two-way doppler 
tracking only, and transmits data at 256 or 128 bits per 
second. The Small Probes are identical. Each cairies three 
scientific instruments and transmits data at 64 or 16 biis per 
second. Their reference frequency for one-way doppler track- 
ing is provided by stable oscillators. Each probe will begin 
transmission about 22 minutes prior to entering the atmo- 
sphere. 

II. DSN Entry Support 

The DSN will meet one of its greatest challenges on Decem- 
ber 9, 1978. In a brief 2 hours, the 64-meter Deep Space 


Stations (DSS) at Goldstone, California (DSS 14) and Can- 
beira, Australia (DSS 43) must acquire and collect data in real 
time from four probes that have never been acquit jd in flight. 
After the probes reach the surface, DSSs 14 and 43 must sup- 
port the Bus entry along with any surviving probes. The Bus 
will be supported by two 26-meter DSSs during the probe 
descent phase. 

III. DSN Implementation at DSSs 14 and 43 
for Multip.wbe Entry 

The following special implementation was required for the 
Pioneer Venus Multiprobe entry: 

(1) Spectral Signd Indicators (SSI): 

Three special spectrum indicators (SSI) with micro- 
processor controls were designed and installed to 
rapidly identify the probe signals and provide accurate 
frequencies to acquire lock on the DSN closed loop 
receivers. The SSI is the key to performing the mission. 

(2) Multiprobe Simulator: 

The simulator was required for mission test and train- 
ing. It generates the four Probers S-band signals, and 
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simulates the entry doppler profile and telemetry bit 
rate changes in accordance with the mission timeline. 

(3) Open Loop Receivers (OLR) 

Two CLRs were added making a total of four 300- 
kHz bandwidth OLRs. They are the signal input to the 
SSI and provide predetected signals to the analog 
recorders for telemetry data recovery. 

(4) Venus Wind Experiment Receiver/Recorder: 

A special 2-mHz open loop receiver and a 12-megabit- 
per-second digital recording assembly that simultan- 
eously records the signals from the four Probes and Bus 
plus two calibration tones were installed to support 
measurement of the Venus winds using long baseline 
interferometry techniques. 

(5) A fifth closed loop receiver (RCVR 7) and subcarrier 
demodulator assembly (SDA 7), and three Program- 
mable Oscillator Control Assemblies (POCA) were also 
installed. 

(6) The telemetry processor software was modified to 
accommodate dual-channel convolutional decoding in a 
single computer. This added the capability of process- 
ing the four Probe data streams simultaneously. 


IV. Entry Procedure Development 

Entry procedures were developed by the Pioneer Network 
Operations Project Engineer (NOPE) and three DSS shift 
supervisors who will support the actual entry event. With a 


preliminai>' version of the procedure complete, the authors 
went to DSS 14 and performed the sequence using the proce- 
dure with regular station personnel. Mtcr three days of test- 
ing and rewriting, the procedures were ready for use in train- 
ing tests. One major discovery was that station operations 
during the entry event required a crew complement of 16 to 
17 operators as compared to a normal crew complement of 
5 to 7 operators It was decided that entry would be sup- 
ported by combining two crews plus additional personnel from 
other crews. 

V. Unique Entry Support Features 

A. Mission Control 

The probe entry phase, by its very nature, precludes any 
project direction. The DSS must therefore support this phase 
0 .heir own. 

B. Dual Uplinks to the Large Probe 

To maximize the probability of acquiring portentry two- 
way lock with the Large Probe in a minimum amount of time, 
both DSS 14 and 43 will simultaneously transmit to the 
probe. The real-time doppler data will be utilized to determine 
which DSS has captured the probe receiver. 

VI. Entry Sequence of Events 

In Table 1, the uncertainties in probe event times are 
(1) ±22 seconds for preentry events, and (2) ±60 seconds for 
entry 'and postentry events. 
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r ^ 


Time 

(GMD DSS Event 


1 . 1530 Large Probe receiver on. 


Time 

(GMT) DSS Event 


2. 

1600 

14 

Antenna on point-Bus in view. 

3. 

1618 

14 

Perform Large Probe uplink acquisition. 

4. 

1730 

43 

Venus rise. 

5. 

1748 

14/43 

Perform minus 60-min checklist. 

6. 

1825 

14/43 

Start prifiie analog (PPR) and digital 
(DRA) recorders Put doppler data 
to line. 

7. 

182613 


Large Probe carrier on. 



14 

Acquire on SSI, lock receiver 3 
(2-way doppler). 



43 

Acquire on SSI, lock receiver 3 
(3-way doppler). 

8. 

1830 

14/43 

Start backup PPR and DRA recorders. 

9. 

183113 


Small Probe 2 carrier on. 



14/43 

Acquire on SSI, lock receiver 2. 

10. 

183113 

14/43 

Large Probe subcarrier on; lock telemetry. 

11. 

183401 


Small Prube 3 canier on. 



14/43 

Acquire on SSI, lock receher 7 

12. 

183613 

14/43 

Small Probe 2 subcarriei on; lock 
telemetry. 

13. 

183708 


Small Probe 1 carrier on. 



14/43 

Acquire on SSI, lock receiver 1 

14. 

183901 

14/43 

Small Probe 3 subcarrier on; lock 
telemetry 

15. 

184013 

43 

Reconfigure for Large Probe bit rate 
change. 

16. 

184313 

43 

Large Probe bit rate change; lock 
telemetry. 

17. 

1845 

43 

Transmitter on-start Large Probe 
postentry uplink sweep. 

18. 

1846 

14 

Start Large Probe postentry uplink 
sweep. 

19. 

184813 

14/43 

Large Probe entry; receiver 3 out of lock. 

20. 

184825 

14/43 

Acquire on SSI, lock receiver 3 
(1-way doppler). 

21, 

184830 

14/43 

Large Probe bit rate change; lock 
telemetry. 

22. 

1851 

14/43 

Receiver 3 out of lock; acquire on SSL 
lock receiver 3 (2-way doppler). 

23. 

185130 

TRK 

Inform DSS of correct doppler node; 
3-way DSS turns rff transmitter. 


24. 

185313 

14/43 

Small Probe 2 entry; receiver 2 out of 
lock. 

25. 

185326 

14/43 

Acquire on SSI, lock receiver 2. 

26. 

185601 

14/43 

Small Probe 3 entry; receiver 7 out of 
lock. 

27. 

185614 

14/43 

Acquire on SSI, lock receiver 7. 

28. 

185918 

14/43 

Small Probe 1 entry; receiver 1 out of 
lock. 

29. 

185931 

14/43 

Acquire on SSI, lock receiver 1. 

30. 

190413 

14/43 

Large Probe parachute jettison. Possible 
receiver 3 out of lock due to doppler 
shift. 

31. 

1906 

43 

Reconrigure for Small Probe 2 bit rate 
change. 

32. 

1909 

43 

Reconfigure for Small Probe 3 bit rate 
change. 

33. 

190913 

14/43 

Small Probe 2 bit rate change; lock 
telemetry. 

34. 

1912 

43 

Reconfigure for Small Probe 1 bit rate 
change. 

35. 

191201 

14/43 

Small Probe 3 bit rate change; lock 
telemetry. 

36. 

191518 

14/43 

Small Probe 1 bit rate change; lock 
telemetry. 

37. 

1925 

14/43 

Prime PPR tape change. 

38. 

1930 

14/43 

Backup PPR tape change. 

39. 

1943 

14/43 

Large Probe impact. 

40. 

1945 

14/43 

Prime DRA tape change. 

41. 

1948 

14 

Reconfigure receiver 3 for Bus; lock 
telemetry. 

42. 

1949 

14/43 

Small Probe 2 impact. 

43. 

1950 

14/43 

Backup DRA tape change. 

44. 

1953 

14 

Bus bit rate change; lock telemetry. 

45. 

1953 

14/43 

Small Probe 3 impact. 

46. 

1954 

43 

Reconfigure receiver 2 for Bus - lock 
telemetry. 

47. 

1955 

14/43 

Small Probe 1 impcct. 

48. 

2024 

14/43 

Bus entry. 

49. 

2025 

14/43 

Prime PPR tape change. 

50. 

2028 

14/43 

Bus burns up; loss of signal. 
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This article reports on activities of the DSN Network Operations Organization in 
support of the Helios Project from 15 Augi s: 1978 through 15 October 1978. 


I. Introduction 

This article is the twenty-fourth in a continuing series of 
reports that discuss Deep Space Network suppoit of Helios 
Mission Operations. Included in this article is information on 
the upcoming eighth perihelion of Helios- 1 and the sixth 
perihelion of Helios-2, science experiments, 22-bii error 
polynomial code (EPC) testing and other mission-related 
activities. 


II. MIsskKi Operations and Status 

On 7 September 1978, commands were sent to Helios- 1 to 
turn on Experiment 9. Experiment 9 had been off since 2 July 
1978 due to a limitation in available regulator output power. 
The expenment was fully configured and science data had 
been received for more than one hour when DSS 44 (Honey- 
suckle Creek, Australia), which was tracking the spacecraft, 
reported a loss of downlink at 22:50:25 Universal Time 
Coordinated (UTC) on 7 September. As a result, a spacecraft 
emergency was declared at 23:10 UTC and 64-meter coverage 
requested. On 8 September at 00:39:37 UTC, OSS 43 
(Tidbinbilla, Australia) acquired the downlink. After evaluat- 


ing the first telemetry data, it was discovered the spacecraft’s 
tenth regulator switch had occurred, which caused the turn-off 
of vhe non-essential bus, changed numerous status bits and off 
pointed the High-Gain .Antenna. After sending several hundred 
commands, the spacecraft and its experiments, except for 
Experiment 9, was reconfigured with the High-Gain Antenna 
back on Earth point. On 8 September 1978 at 09:00 UTC, the 
emergency was lifted. Two possible caases foi this regulator 
switch have been identified: 

(1) Th^' switch might have been triggered by a solar event 
that was reported at the Space Environments Services 
Center, Boulder, Colorado, at approximately the same 
time p^rioQ. 

(2) Experiment 9 does not have a current limiter for 
protection and a stuck stepper motor might have 
cause j an overload condition. Experiment 9 will be left 
off until the circumstances can be clarified with the 
experimenter. 

On 26 October 1978, the Helios-1 spacecraft will enter its 
eighth perihelion period, which will last until 3 November. 



Mere information concerning this event will be recorded in 
future articles. 

The Hclios-2 spacecraft has performed nominally for this 
peifod with no signiflean: occurrences or problems. The 
spacecraft entered its sixth perihelion period on 8 October 
1978. Helios>2 entered blackout on 9 October while being 
tracked by DSS 67/68 (Weilheim, Germany). This earlier han 
expected blackout v'as probably due to high solar activity. The 
spacecraft was at a Sun-Earth*Probe (SEP) angle of 4.1 degrees 
at blackout. Assuming a symmetrical exit angle, the spacecraft 
was configured for a 13-da> unattended operation period with 
memory storage covering. Helios-2 was configured in data 
mode 4, format 3, bit rale of 8 bits per second (bps) and se- 
quence 17. Blackout exit is scheduled to occur on 22 October 
1978. As with Helios-1, more information will be available for 
a future aiticle. 

Helios-2 also passed its 1001st day in orbit on 12 October 
1978, 


Overall coverage for Helios- 1 and Helios-2 for this period is 
listed in Table 1. 

IN. Special Activities 

A. Support of On-Board and Ground Expet^ments 

In conjunction with the perihelion and superior conjunc- 
tion periods, science data will be collected during times 
discussed earlier. The prime data types will be faraday rotation 
data and ellipticity data for Experiment 12. Also solar wind 
data will be collected during the same period. 


B. German Space Operations Center (G&'OC) 
Converaion to the 22-Bft Error Polyiiomial Code 
(EPC) and the Mark III Command System 

As reported in the last article (Ref. 1), testing has con- 
tinued on the 22-bit EPC interfaces and Mark 111 Command 
System throughout this period. Overall, the testing has gone 
well. On 3 October, the first demonstration track was per- 
formed by DSS 44 while tracking Helios- 1 to verify that GSOC 
was capable of processing telemetry, monitor, arid command 
data while in the 22-bit EPC mode. The monitor software at 
the Network Oper^^tions Control Center (NOCC) was, unfor- 
tunately, not available in time for this test, and the Mark III 
Command System performance was hampered by certain 
anomalies discovered in the station’s Command Processor 
Assembly’s test software. These anomalies affected command 
transmission to the spacecraft and are being corrected. The 
interfaces between GSOC, JPL, and DSS 44 worked well in the 
22-bit EPC mode, and telemetry and command data flowed 
without difficulty. Another demonstration was conducted on 
5 October with basically the same results. On 8 October, 
demonstration tracks with DSS 62 (Cabreros, Spain) tracking 
Helios-1 and DSS 63 (Madrid, Spain) tracking Helios-2 were 
conducted. The purpose here was to verify that GSOC could 
process both Helios- 1 and HePos-2 telemetry data, while in the 
22-bit EPC configuration. This test was successful with no 
problems encountered. The most recent test was a telemetry 
and monitor test with DSS 44 tracking Helios- 1 and DSS 14 
(Goldstone, California) tracking Helios-2. This test followed 
the same foimat as the telemetry test above, but with the 
addition of one monitor data stream from each station. This 
test proved successful and GSOC processed both telemetry and 
monitor data from each station. Future tests (demonstration 
tracks) are scheduled to check out the Mark III Command 
System folloving the correction of DSN software. These tests 
Vriil be reported a mture r-port. 


References 


Goodwin, P. S., Jensen, W. N., and Rockwell G. M., “Helios Mission Support,” in The 
Deep Space Network Progress Report 42-47 y pp. 26-28, Jet Propulsion Laboratory, 
Pasadena, CA, 15 October 1978. 



Tflbtel. HsHm tracking 


Number 


Month 

Spacecraft 

Station 

type 

of 

tracks 

Tracking time 
(h:min) 

August 

Helios-1 

26 meter 

31 

147:13 



64 meter 

0 

00:00 


Helios-2 

26 meter 

4 

20:47 



64 meter 

25 

84:41 

September 

Heiiosrl 

26 meter 

27 

139:40 



64 meter 

2 

8:25 


Helios-2 

26 meter 

0 

00:00 



64 meter 

21 

91:02 
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‘ 77iS arti^e presents results of simulation of several shbrt-cdnstrdinMength convolu- 
tional codes using a noisy symbol stream obtained via the turnaround ranging channels of 
the MVM*73 spacecraft First operational use of this coding technique is on the Voyager 
mission. The relative performance of these codes in this environment is as previously 
predicted from computer-based simulations. 


This short article presents the results of simulations of 
several short-constraint-length convolutional codes using 70 
min of noisy 5.8-kbps symbol stream recorded by the 
930 computer at DSS 14 on Jan 15, 1974. This data was 
obtained as part of an experiment intended to define the 
performance to be expected from X-band telemetry on the 
Voyager and subsequent missions. The simulated telemetry 
stream Was originated on the ground and returned vii the 
turnaround ranging channels of the MVM’73 spacecraft at 
both X- and S-band (Ref. 1). The same experiment was later 
performed witli Viking (Ref, 2). 

The data stream consisted of the repeated sequence 
111010, modulated on a subcarrier. This sequence is a 
codeword of the K=7, rate 1/2 convolutional code now flying 
on Voyager, and also of two rate 1/3 codes at K=4, and K=5. 
The signal was processed at the time of reception as an 
uncoded bit stream and also via a hardware Viterbi decoder for 
the K=7, rate 1/2 code. It had been initially intended to also 
operate the JPL-designed K=5, rate 1/3 decoder (Ref. 3) in 
real-time, but this was not in fact pursued. The data recorded 
by the on-site SDS-930 computer as 6-bit quantized symbols 
was subsequently converted to a stream of 4-bit quantized 


channel-noise samples by extraction of the known data stream. 
This conversion was done to simplify processing in software 
through a variety of short-constraint-length convolutional 
codes. The use of 4-bit quantization compressed this data to 
two virtually full magnetic tapes, one each for S-band and 
X-band data. 

Decoding from these tapes was performed using a software 
decoder in the Sigma 5 computer, at constraint lengths K-3 to 
K=9, rates 1/2 and 1/3 (Ref. 4 describes a 930 version of this 
decoder). The data was blocked into 5-min segments of 
approximately 2X10^ symbols to obtain statistical averages 
for the coded without masking variations that had been 
observed in the signal-to-noise ratio during the real-time 
experiment. Figure 1 shows a time-trace of the decoding 
results for the S-band data at rates 1/2 and 1/3. Figure 2 shows 
a similar time trace for the X-band data, decoded at rate 1/2 
only. The SNR at S-band was too high tc permit statistically 
significant results for the rate 1/3 codes. 

The traces for several of the codes are not complete, work 
having been stopped pending resolution of the difference 
between hardware and software decoder performance. This 
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problem has since been resolved (Ref. 5) and the software 
results shown to be c )rfC4.t, at least at K=7. rate 1/2, X-band. 

Figure 3 summarizes the coded bit error probability seen 
for these codes rs a “function” of the signal-to-noise ratio 


estimated by way of the (uncoded) symbol error rate. There 
are no real surprises in this figure. The performance of all 
codes tested, either rate 1/2 or rate 1/3, is comparable to their 
performance as seen in earlier simulation with computer- 
generated noise (Ref. 6). 
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Microwave Radiometer Measurement of Water Vapor 
Path Delay: Data Reduction Techniques 

E. S. Claflin, S. C. Wu, and G. M. Resch 
Tracking Systems and Applications Section 


Unmodeled electrical path delay from atmospheric water vapor is a limiting error 
source in geodetic measurements made with very long baseline interferometry and in 
radio ranging to spacecraft A dual channel water vapor radiometer, operating near the 
22,235-GHz water vapor line, is capable of measuring water vapor-induced delay with 
good accuracy under most weather conditions. Theory shows that water vapor path delay 
AL^ is proportional to a linear combination of saturation-corrected sky brightness 
temperatures, measured on and off the water vapor line. The second, off-line, channel 
removes the effects of emission from liquid water droplets in clouds as well as most of the 
oxygen emission. Tipping curves remove instrumental error. Sky brightness temperatures 
are saturation-corrected or **linearized ** using estimates of effective sky temperature made 
from surface temperature. Coefficients in the expression for path delay AL^ are functions 
of surface temperature, pressure, and water vapor density, allowing use of our data 
reduction algorithm at any altitude and in any climate. Coefficients are found by two 
methods: (l)from a regression analysis of measured brightness temperatures versus 
radiosonde measured delay, and (2) from a regression of theoretical brightness tempera- 
tures versus radiosonde measured delay. Regression solutions are constrained to remove 
liquid water contributions and to ^ve the correct slope (i.e., one) for radiometer versus 
radiosonde path delay. 


I. Introduction 

Very long baseline interferometry (VLBl) and other extra- 
terrestrial, microwave, geodetic techniques have potential 
accuracies at the few centimeter level (Ref. 1). Propagation 
media effects constitute a limiting error source, however, and 
must be compensated for to achieve the highest accuracy 
otherwise inherent in present and platmed systems. In particu- 
lar, one must devise a way to compensate for tropospheric 


water vapor. Water vapor induced path delay in tiie zenith 
direction can vary from less than 1 to more than 30 centimeters 
depending on climate and weatlier. 

Under favorable conditions, water vapor delay can be 
solved for as part of tlie final fitting of residuals in the VLBI 
data reduction process. If baseline geometry is not favorable, 
the solve-for procedure will degrade accuracy appreciably. 
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Moreover, the solve-for solu ion breaks down if changes in 
water vapor delay occur rapidly or if the atmosphere has large 
horizontal inhomogeneities. An independent way to measure 
delay from tropospheric water vapor is needed. 

From among a number of candidate measuring techniques, 
we have found tliat the most satisfactory is passive microwave 
radiometry. A microwave radiometer will operate day or night 
and under all weather conditions except rain or heavy snow. 
Spatially, measurements can be made in all directions that are 
more than a few beamwidths from the horizon and from the 
sun. Microwave radiometry will, in the future, also be 
applicable to VLBI geodesy with satellite radio sources and to 
radio ranging for navigation of interplanetary spacecraft. Hie 
purpose of this paper is to present techniques for conversion 
of radiometer data to water vapor path delay. A discussion of 
errors will be presented in a later paper. 


II. Background 

A solution of the radiative transfer equation governing 
propagation of microwave radiation downward througli a 
nonscattering atmosphere in thermodynamic equilibrium is 
(Ref. 2) 



is the optical depth, cc is the absorption coefficient 
(distance'^), s is distance upward along the signal path, T is 
physical temperature of the air, Tq is the ^ky brightness 
temperature, and is the cosmic biackbody temperature of 
about 2.8 K. Typical half-power antenna beamwidths of 5 to 
10 degrees are large enough that all celestial sources of 
radiation are negligible except for the sun. 

Important contributions to a are made by atmospheric 
water vapor, oxygen, and the liquid water droplets in clouds or 
fog. Water vapor has a spectral line at 22.235 Ghz. Sensitivity 
to water vapor is achieved by choosing a primary frequency on 
or near this line. The liquid water contribution from clouds is 
removed by using a second frequency channel and then 
linearly combining sky temperatures from both channels. We 
find that it is sufficiently accurate to assume an P frequency 


dependence for the liquid \^*ater absorption coefficient,^ The 
oxygen contribution is from a complex of spectral lines near 
60 GHz and has a frequency dependence which is also 
approximately p. About 82 percent of tlie oxygen contribu- 
tion is removed by the p dependence assumed for liquid 
emission. This is fortuitous because oxygen emission also 
depends on atmospheric tempeiature, a quantity subject to 
estimation error. The remaining 18 percent of the oxygen 
contribution is removed by a small term proportional to air 
mass in the expression for water vapor path delay. At one air 
mass this term is approximately -0.4 cm path delay. 

A related matter that has already been investigated in some 
detail is optimal choice of frequency pair for the two 
radiometer channels. To retrieve water vapor path delay with 
highest accuracy, it is necessary to use a frequency pair which 
minimizes error from all sources. No single frequency pair is 
ideal under all conditions; however, we recognize that a very 
important consideration is that a frequency pair be chosen 
which minimizes variation with altitude of the weighting 
function l/(s) (Ref. 3): 


W\s) 


p [/? 


( 2 ) 


where is the absorption coefficient for water vapor and p is 
water vapor density. The penalty for a nonconstant Eq. (2) 
along the observing path is appreciable error whenever the 
water vapor altitude distribution varies markedly from some 
calibration distribution. Such variations appear to be most 
common in coastal regions, where layering of marine and 
continental air is a frequent occurrence. Optimal choice of 
frequency pairs has been discussed by Wu (Ref. 3). Frequency 
pairs of 20.3/31.4, 20.0/26.5, and 24.5/31.4 GHz are all near 
optimal for minimizing variations in weighting function 
(£q. 2) with altitude. 


^ Liquid water absorption coefficient frequency dependence was calcu- 
lated for the case of drcplets small compared with wavelength 
(Rayleigh scattering). The dependence is very nearly f\ but is slightly 
affected by temperature (Refs. 4, 5). For a 22.235/31.4-GHz fre- 
quency pair, deviations from p behavior will cause an error in path 
delay of less than 0.7 cm even under the extreme conditions of 
1000 urn liquid water, 20-degree elevation angle, and temperature of 
(fC. Note also that path delay from liquid water is negligible; 
1000 ^m of vertical columnar liquid water at 20-degree elevation angle 
causes path delay of less than 0.5 cm (Ref. 6). For a next generation 
system, these effects could be estimated to refine the total delay 
correction. 
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III. Outline of Technique 

The data reduction system developed here can be separated 
into two parts or activities. The first part comprises techniques 
for converting raw radiometer data and surface temperature 
into linearized sky brightness temperatures. An important 
element of this procedure is the use of tipping curves to 
provide an instrumental calibration. 

The second part is calculation of suitable coefficients for 
use in finding path delay. Two methods will be discussed for 
calculating coefficients^ both using a regression analysis. In 
each case the coefficients are adjusted for differing climatic 
conditions by measurements of surface temperature, pressure, 
and relative humidity. 


IV. Linearized Brightness Temperature and 
Water Vapor Path Deiay 

Presence of the exponential teims in Eq. (1) prevents sky 
brightness temperature from increasing linearly with increasing 
water vapor path delay. In the limit of small opacities, tlie 
relation is linear. In the limit of extremely high opacities, sky 
brightness temperature becomes equal to surface atmospheric 
temperature (saturates) and is completely insensitive to 
changes in water vapor path delay. Between these limits it is 
possible to "‘linearize” (saturation-correct) brightness tempera- 
ture. The purpose of linearization is to transform the 
calculation of path delay into the familiar problem of solving 
two linear equations in two unknowns. (The two unknowns 
are water vapor path delay and integrated liquid water.) 

If the linearized brightness temperature 7^ is defmed by 


Tads 


(T-T^)ads (3) 



n = r ri-r(oo] 


-)]+ r 

*'0 


with 



( 5 ) 


( 6 ) 


For small opacities, * T»fr For large zenith opacities 
or small elevation angles, Since atmospheric 

temperature is readily accessible only at the ground, surface 
values must be used to estimate 2 A simple expression 




T* 


kj 


(7) 


(where T is absolute surface temperature) is sufficient for dry 
climates, is a function of frequency and can be found by 
integrating radiosonde profiles with Eq. (6). Typically 
0.92 < kg <0.95. More accurate, but more complex, schemes 
for estimating and 7gyy can be devised. 

7^, the linearized brightness temperature, is not very 
sensitive to Jeff and T'eff: Figure 1 shows the linearization 
correction 7^ - 7^ as a function of Tq for tliree different 
estimates of Jff~ Jfr The correction becomes a significant 
error source only at very high sky temperatures. Observed 
zenith sky brightness temperatures at 22.235 GHz have been 
in the range 13 to 100 K. At 31 .4 GHz the range has been 1 1 
to 45 K. The corresponding range of zenith path delays 
covered by these observations was from 2 to 30 cm. 

Wu (Ref. 3) has shown that for a dual channel radiometer, 
path delay is given by 


then linearized briglitness temperature may be calculated from 
measured brightness temperature 7^ using the following 
relation derived from Eqs. (1) and (3) (Ref. 3): 




where 


( 8 ) 







T,)ln 
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(9) 

m(PjPf (f/r)‘ ®* 

(10) 





“ constants 

P = surface pressure: force/area 
T = surface temperature 

j and Tq 2 are derived from measured brightness 
temperatures at the two frequencies using Eq. (4). M/'(0) is the 
weighting function of Eq. (2) evaluated at the surface to adjust 
path delay for varying values of surface temperature, pressure, 
and water vapor density. The variable G adjusts for changes in 
air mass m and for changes in the oxygen absorption 
efficient with changes in surface temperature and pressure. 
T and P are nominal surface values, but are best chosen to be 
representative of the calibrating or operating conditions. 
Expressions for the water vapor absorption coefficient can 
be found in Waters (Ref. 2) and Westwater (Ref. 7). Tlie 
coe'^'i'-ients by, b^ can be determined either from theory 
and a suitable ensemble of radiosonde profiles or from 
regression of measured sky temperatures against path delays 
from simultaneous radiosonde data. These methods will be 
discussed in a later section. 

For convenient routine application, Eqs. (8), (9), and (10) 
can be written as 


^0'^ ~Tc^ ‘‘l Tb,\ '*'‘*2 Tg2 

(11) 



” h^'(O) 

(12) 

y. ^ C 1 

(13) 



= bJW'iO) 

(14) 

flj = bJW\Q) 

(15) 


and a table of coefficients ^2 prepared for suitable 

choices of P, T, and p values. We have found that sufficient 
resolution is provided by increments of 20 mb (2000 N/m^) in 
P, 10®C in r, and relative humidities of 20, 50 and 80 percent. 

V. Measurement of Sky Brightness 
Temperature 

Radiometer output is an analog signal which varies linearly 
with input signal power. Since the constant of proportionality 


is gain dependent, every radiomc;ter contains two reference 
loads at carefully measured physical temperatures to establish 
a brightness temperature calibration. If the sky brightness 
temperature is Tq, the ambient or base load temperature 
and the hot load temperature then the corresponding 
number of counts (from an analog to digital 

converter) are related to bri^tness temperature by 

T, = T^HT„-T^yr (16) 

y-iNg-N^)HN„-N^) (17) 


In practice Eq. (16) is imperfect because of losses and 
reflections in tlie horn, waveguide, and swi*.jhing system. 
Figure 2 shows typical switches and signal paths for a 
radiometer. Tb’ are apparent brightness temperatures 

at the common input. We assume a priori that switch isolation 
is very good (typically the isolation is ~30 dB). Tlien 


Tb ~ Tb^b^ Tgi 

(18) 

Th ~ 

(19) 

II 

(20) 


are expressions for the input temperatures for a simple linear 
model in which the input signal suffers attenuation and 
receives augmentation Til- The and Til^ are temperature 
dependent. Temperature stabilization of the instrument should 
remove most of the variation in these parameters. Equa- 
tion (20) assumes that file switch and waveguide temperatures 
are nearly the same as tlie ambient load temperature . 

Solving Eqs. (16), (18), (19), and (20) for Tg, 

h = fs [<^ - ■ r<>7] (21) 

Determination of the four instrumental constants 

and Thl is not straightforward either theoretically or 
experimentally. A simplified version of Eq. (21) can be written 
by assuming ^ 1 and Tbl=^^> which implies negligible 
losses and reflections between antenna and receiver input. The 
instrumental constants are then replaced by the “hot load 
correction,” AT// = ’ 1 ) + so that 

+ ( 22 ) 
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The single constant ATf^ can be found by performing 
“tipping curves'* as outlined in the following section. Some 
variation in AT^r will occur with different brightness and 
physical temperatures because £q. (22) is only an approxima- 
tion to Eq. (21). 


VI. Tipping Curves and the Hot Load 
Correction 

Tipping curves are performed by taking measurements of 
sky brightness temperature at several elevation angles. When a 
St aight line is fitted to linearized briglitness temperatures 
from Eqs. (22) and (4) as a function of air mass, the expected 
intercept at zero air mass is Suppose the intercept has a 
different value r^(0). AT^ must be adjusted to achieve the 
desired intercept. Figure i sliows a straight line fit before 
(dashed) and after (solid) adjustment of AT//. (The two lines 
intersect off the graph at Tq = .) 


If the hot load correction is initially taken to be zero, then 
one can show that 




(r;(o)- 


(23) 


Experience with an instrument will show to what degree 
ATff is constant for each channel One can then estimate error 
incurred by using nominal values of to replace tipping 
curves. 

In principle, the most meaningful delay is found by 
measurements along the line-ofsight to the radio source being 
observed. For an automated system, a horizon mask for f <ch 
site and a solar ephemeris are needed to avoid pointing the 
instrument toward contaminating radiation. Under dry condi- 
tions, pointing along the line-of-siglit is not likely to give 
measurable improvement over extrapolated zenith values. In 
fact, error may be introduced by sidelobe pickup and poorer 
performance of the data reduction algorithm at low elevation 
angles. We expect pointing to improve VLBI results under 
very humid, cloudy conditions. 


VII. Coeffsient Determination: Regression 
of Sky Brightness Temperatures Versus 
Radiosonde Path Delay 

Water vapor delay is computed from a radiosonde profile 
with the expression (Refs. 8 and 9). 


ALj = (1.723 X 10*^ K-m^/g) I p/Tds (25) 

•'o 


is a good estimate of the required hot load correction. 
Corresponding sky brightness temperature correct ^''n is found 
from Eq. (22) to be 



New sky brightness temperatures + ATq are linearized 
with Eq. (4) and fitted again to a straight line. The process is 
repeated^ as many times as necessary to adjust Tj^(O), the 
intercept, to within about 0.1 K of Normally two or three 
iterations are sufficient. We generally use three points at 1, 
1.5, and 2 air masses (90-, 42-, and 30-degree elevation angles) 
for tipping curves. Repeating the measurement at four points 
of the compass (different azimuths) provides redundancy in 
the data and averages over spatial gradients. 

Various problems witlt the measurement system can cause a 
poor tipping curve. Quality criteria for a tipping cur*e are the 
value of the linear correlation coefficient, the size of the hoi 
load correction and its agreement with a “historical** estimate, 
and the agreement of the estimated zenith paili delay taith 
those derived from tipping curves at other azimuths. 


where the prime indicates a radiosonde measurement. 


^Equations (23) and (24) must be generi/lized slightly to take into 
account each new adjustment of the hot load temperature t ' h - 
If 


^ Hi ' ^ Hi 

i 




pine 

*Hi ’ ^ “^/// 
/=0 


and 


then 








Vh-'^aI 


(23.) 


(24a> 


where /■ 0, 1, 2, ... is the iteration number and **inc** is incremental 
change. 
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SuppoM a suitaHy representative group of radiosonde 
launches has betn made, and sky brightness temperatures 
measured with a radiometer at the time of each launch. The 
coefficients b^, b^, b^ of Eq. (8) and Eq. (9) may be found by 
minimizing: 

( 26 ) 


with respect to b^, and b^y but subject to several 
constraints. 

The fust constraint ensures that background contributions 
are eliminated: 

6A“-(V-Ci)* (22) 

Two more constraints enforce unity slope and zero intercept 
on the best-fit straight line of radiometer- measured to radio- 
sonde-measured path delays:’ 

• Z“’., < 2 «) 


= E(^w)* (29) 

/ / 

If these constraints are not used, slope and intercept may 
vary appreciably from 1 andO. This unphysical result is a 
necessary consequence of random measurement error in sky 
brightness temperatures and in radiosonde-measured padi 
delays: Scatter in the residuals of Eq. (26) will be minimized 
by values of b^ and which are too small, offset by a value of 
b^ wiiich is too large. Figures 4A and 4B illustrate fits obtained 
with and without constraints on slope and intercept. The 
coefficients of Fig. 4B will predict too much water vapor path 
delay under dry conditions and too little under wet conditions. 


^One may show that only the constraint on slope is actually required. 
In practice it is easier to constrain both slope and intercept. 


Lagrange multipliers may be used to minimize Eq. (26 j with 
respect to, b^, bj, bj, but subject to the constraints of 
Eqs. (27), (28), (29). The problem is thereby simplified to a 
set of six linear equations in six unknowns. The regression may 
ako be carried out without actual measurement of sky 
brightness temperatures. Linearizec brightness temperatures 
can be calculated from Eq. (3) usin)' the water vapor absorp- 
tion coefficient of Waters (Ref. 2) or Westwater (Ref. 7) and 
the oxygen absorption coefficient of Westwaier (Ref. 7).^ As 
before, one uses a representative sample of radiosonde profiles. 

There is sufficient uncertainty in the water vapor absorp- 
tion coefficient to cast some doubt on the accuracy of results 
obtained by the theoretical method. Nevertheless, we have 
found that the theoretical method gives coefficients in fairly 
good agreement with those obtained by regression against 
measured brightness temperatures. 


VIII. Summary 

A systematic method for reducing dual-frequency water 
vapor radiometer data to water vapor path delay has been 
developed. Linearized sky brightness tenip^ratures and path 
delay coefficients are computed in separate procedures and 
then combined to calculate path delay. Sky brightness 
temperatures may be measured by using either tipping curves 
directly or by using nominal hot load corrections based on 
tipping curves. Coefficients may be derived from brightness 
temperature data using a constrained regression analysis 
against radiosonde-measured path delays, or from a con- 
strained regression analysis of tlieoretically calculated sky 
brightness temperatures against radiosonde measure ' pjth 
delays. Hie coefficients are functions of surface temperature, 
pressure, and water vapor density, to allow use of the data 
reduction algorithm at any altitude and in any climate. 


^ Three corrections should be made to Wcstw?f#«r*s oxygen absorption 
coefficient. The pressure should be partial pressure of dry air only, 
except in the line width where it is total pressure. The coefficient of 
the line width should be changed to 0.025 from 0.049. The powei law 
temperature dependence of the line width should be changed to 0.85 
from 0. 75 (Westwater, E. R., private communication). 
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Rg. 3. Tipping ctirvM with no hot food correction (dashed line) 
and with hot ioaa correction required to give desired zero air mass 
intercept (solid line); hypothetical data points for each case are 
also shown 
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Fig. 4. Path delay ALr from radiometer measurements vs path delay m; from radiosondes 
using regression: (A) with constraints on slope and intercept, (B) without constraints on slope 
and intercept A circle Is used to Indicate two overlapping points; the straight lines are least 
squares fits to the poinU generated by each regression 
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Low-Noise Receivers: S-Band Parametric 
Upconverter Development 

S. Petty, D. Neff, and D. Norris 
Radio Frequency and Microwave Subsystems Section 


The combination of a cryogenicaUy<ooled parametric upconverter and a higher 
frequency maser post amplifier has been proposed as a method of achieving maserdike 
receiver noise temperatures over much larger instantaneous bandwidths and tuning ranges 
than are presently obtainable with masers in the range of I to 18 GHz. An experimental 
2.0- to 2.5-GHz parametric upconverter Imaser system has been developed to explore 
these possibilities. Initial tests of this system have resulted in an effective input noise 
temperature of 3.1 K at 2295 MHz and 3.2 K at 2388MHz. The parametric upconverter 
has logged over 1500 hours at 4.5 K and has undergone 5 thermal cycles (300 K to 


4.5 K to 300 K) without degradation. 

I. Introduction 

Recent advances in many areas of technology have made 
possible the development of an ultralow-noise receiving system 
that can exhibit roise temperatures similar to those of micro- 
wave maser amplifiers, while providing instantaneous band- 
widths and tuning ranges many times greater than those of 
present maser amplifiers in the 1- to l8-GHz range. This 
system is comprised of a cryogenically-cooled, upper-sideband 
parametric upconverter (Ref. 1) followed b> a maser amplifier 
operating at a much higher frequency where wide bandwidth 
and large tuning ranges have beeu achieved. Two recent devel- 
opments in particular have made these performance levels 
possible: the first is the availability of very high-quality 
gallium arsenide varactor diodes with low package parasitic 
reactances, and the second is the recent development of a 


wideband microwave maser amplifier (Ref. 2), which tunes 
over the range of 19 to 26 GHz. 

To evaluate the feasibility of this type of system for future 
use in the DSN, an S-band-to-K-band cryogenically cooled 
parametric upconverter^ and a K-band maser have been 
developed, assembled, and tested at JPL. A block diagram of 
the experimental system is shown in Fig. 1 , and a photograph 
of the completed package is shown in Fig. 2. Following 
sections of this report will describe each of the system compo- 


^The parametric upconverter described in this report was developed by 
AIL, a division of Cutler-Hammer, Melville, L.I., N.Y„ under cont jct 
to JPL. 


nents, provide a theoretical noise temperature analysis, and 
present performance measurements obtained to date. 

II. S-Band Parametric Upconverter 

The ^)-band-to K-band parametric upconverter (Ref. 3), 
shown in Fig. 3, has a coaxial signal input circuit and a varac- 
tor mount, pump circuit, and (sum) output circuit in K-band 
waveguide. The input circuit consists of a Type N input con- 
nector, a bias tee, which provides diode bias voltage via the 
mput coaxial center conductor, and a matching and tuning 
section. The varactor diode mount contains a pair of matched 
diodes mounted across a K-band v/aveguide in a balanced 
configuration. This configuration isolates the input circuit 
from pump, output, and image frequency components, there- 
fore eliminating the need for input circuit filtering with the 
attendant insertion loss. Separation of pump and output 
frequency components is accomplished by two bandpass 
filters. Waveguide spacers and reduced-height sections (not 
labeled in Fig. 3) provide impedance transformation and 
tuning at the output, pump, and image frequencies. 

The upconverter measures 8 dB of net gain with a 3-dB 
instantaneous bandwidth of approximately 500 MHz when 
operating at a physical temperature of 4.5 K. Input return loss 
measures -7 dB or better from 1950 to 2450 MHz, with the 
region of best input match (return loss >-10 dB) between 
2250 and 2420 MHz. Note that these parameters, which vary 
considerably with bias current level, were measured at a 
varactor bias current of 1 to 2 juA. Subsequent noise tem- 
perature measurements (described in Section VI of this article) 
have demonstrated that the best midband noise performance 
occurs with zero bias current. Therefore, gain, bandwidth, 
and other RF parameters must be examined after optimum 
bias settings (for best noise performance) have been 
determined. 

The maser pump rejection filter (Fig. 1) located at the out- 
put port of the upconverter is necessary to prevent maser 
pump power leakage from combining with upconverter pump 
power in the varactor diodes to produce spurious in-band 
interfering signals. The present filter is of the “waffle-iron” 
type and exhibits sufficient attenuation at the maser pump 
frequency (and its harmonics) to permit noise temperature 
measurements. An operational upconverter/maser system for 
the DSN (where received signal levels can be as low as -180 
dBm) would require an interstage filter with much higher 
attenuation at these frequencies. 

III. K-Band Reflected Wave Maser (Post 
Amplifier) 

The K-band reflected wave maser is a four-stage, ruby- 


fiUed waveguide traveling wave structure operating at 4.5 K 
(Figs. 4 and 6). A circulator is used to direct input signals 
into each ruby-filled waveguide section where maser amplifi- 
cation occurs; on reaching the end of the filled waveguide 
section, the signal is reflected back to the input and is ampli- 
fied a second time before passing through the input circulator 
and onto the next stage. Additional circulators (Fig. 5) are 
used to increase interstage isolation. The three sections shown 
in Fig. 4 are combined into the single assembly shown in 
Fig. 6 with the superconducting magnet required for maser 
operation at the left. 

When the maser is tuned to the upconverter output fre- 
quency, the maser pump frequency is modulated to cover a 
1-GHz range centered near 51 GHz. Pump power (75 to 
100 mW) is supplied witli a Siemens Corp. model number 
RW060 backward wave oscillator. Overall maser gain is 30 to 
40 dB with 60 to 250 MHz of instantaneous bandwidth 
depending on the frequency of operation. The limited instan- 
taneous bandwidth available at present is due to magnetic 
field nonlinearities in the superconducting magnet. This 
condition does not affect the validity of upconverter/maser 
system noise temperature measurements, but it does require 
that the maser be tuned to the specific frequency of interest 
within the upconverter bandwidth so that maser gain and 
effective mput noise temperature (4.7 ± 0.5 K) are maintained. 

This particular maser is a copy of the K-band reflected-wave 
maser developed by JPL and NRAO^; a complete technical 
description of which is soon to be published (Ref. 2). Con- 
tinuing K-band maser development work by NRAO has 
produced a 500-MHz bandwidth. This improvement was 
accomplished by using a larger superconducting maguc^ with 
the same maser structure as described here (Ref. 4). initial 
development work on the reflected-wave maser concept began 
in 1972 under a California Institute of Technology President’s 
grant to the University of California at San Diego and the Jet 
Propulsion Laboratory (Ref. 5). 

IV. Cicsed-Cycle Refrigerator 

A standard 1-W JPL4.5-K closed-cycle refrigerator (CCR) 
(Refs. 6, 7) is used to house the upconverter/maser assembly. 
Figure 2 shows the packaged system with attached pump 
sources. Figure 7 shows the CCR with vacuum housing, input 
transmission line, and radiation shields removed. 


^The development of this maser design was the result of a joint research 
effort between the National Radio Astronomy Observatory, operated 
by Associated Universities, Inc., under contract with the National 
Science Foundation, and the JPL/CIT under contract NAS7-100 
sponsored by NASA. 
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The S-band input transmission tine (Fig. 8) utilizes an 
existing waveguide-to-coax design (Ref. 8) that cools the 
entire length of coaxial center conductor to 4.5 K. The noise 
temperature contribution from this transmission line is calcu- 
lated to be 0.5 K. A length of 3.5-nun (0.141-inch) O.D. 
semirigid cable, ail at a 4.5-K physical temperature, couples 
the input signal from the CCK input transmission line to the 
upconverter input connector. Two liquid-helium heat 
exchangers are required to properly maintain the upconverter 
diode me mt, superconducting magnet, maser circulator 
assembly, maser structure, and related interconnecting wave- 
guide components at 4.5 K. 

Pump inputs in WR 42 and WR 19 waveguide (Figs. 2, 7) 
i are used for exitation of the parametric upconverter (21.76 

GHz) and the maser ruby structure (51 GHz), respectively. A 
WR 42 waveguide is used for the maser output at 24 GHz. 
Other inputs to the 4.5 K station include the magnet charging 
system and normal instrumentation wires. 

New, larger heat shields were required >>ecause of the 
increased size of the 4.5 K assembly. The increased thermal 
radiation load from this source, combined with the total input 
pump power (120 to 150 mW) and the heat conduction 
associated with the many input connections to the 4.5 K 
station, have severely loaded the 1-W capacity of the CCR. 
Operation in a 35®C ambient environment and rapidly 
repeated retuning of the superconducting magnet have resulted 
in several CCR warmups. The system, while functional for 
testing, is very marginal. The development of CCRs with 
greater 4.5 K cooling capacity is planned. 

Tlie large mass kg) attached to the 4.5-K station 
required changes in the support structure and precool system. 
Mechanically, the upconverter/maser assembly and 4.5-K 
station are supported by the input waveguides and a solid 
copper block connected to the 15-K station via the hydrogen 
thermal switch (Fig. 7). A larger nitrogen precool heat 
exchanger and copper sSieath were attached directly to the 
superconducting magnet assembly since heat transfer through 
the Hyperco 27 magnet material is low. Additional copper 
straps are used to aid heat transfer. These changes have 
reduced the precool time (ambient to 80 K) from three hours 
to one hour; a total of eight hours is required to cool down to 
4.5 K. The 15-K station appears to have additional cooling 
capacity to decrease the cooldown time; a larger capacity 
hydrogen thermal switch should be developed to improve 
^ heat transfer from the 4.5-K station to the 15-K station during 

j the cooldown process. A new digital thermometer^ is being 


? ^ Model DRC-70C, Digital Therm ometer/Controller manufactured by 

Lakeshore Cryotronics, Inc.. Columbus, Ohio. 



used to monitor CCR temperatures replacing the traditional 
thermocouple sensors. The direct readout feature from 2 to 
300 K enables more convenient and more accurate measure- 
ments than has previously been achieved with thermocouple 
instrumentation. 


V. Theoretical Noise Temperature Estimate 

Table 1 itemizes the calculated contribution to the effec- 
tive input noise temperature (at the room temperature wave- 
guide input flange) of each component in the upconverter/ 
maser system as diagrammed on the left side of the table. 

These calculations do not include the effect of input or output 
circuit mismatch. Therefore, the results are applicable only to 
midband frequencies where upconverter gain is maximum and 
input and output reflection losses are small. At 2295 MHz, 
upconverter gain is within 1 dB of maximum and input return 
loss between -10 and -13 dB (depending on bias setting), 
and it is valid to compare noise temperature measurements 
obtained at this frequency with the theoretical values in 
Table 1. 

As indicated in the table notes, the parameter values used 
in the calculation are best estimates based on handbook values, 
previously measured data on similar parts, or data from the 
various references. Upconverter input and output losses 
and are estimated from (1) measurements reported in 
Ref. 3, and (2) the difference between theoretical gain (assum- 
ing no circuit losses but including calculated diode losses) and 
actual measured gain. 

The physical temperature of the diode junctions has a sub- 
stantial effect on the overall system effective noise tempera- 
ture according to Table 1 : 2.5 K if the diode junctions experi- 
ence no pump heating, and 3.5 K if the diode junctions are 
heated to 20 K. The uncertainties in the estimates of the 
various contributions listed in Table 1 have not been rigor- 
ously analyzed. However, it is believed that the sum of the 
worst case errors for all other contributions does not exceed 
the ±0.5 K resulting from the two different diode junction 
temperature assumptions discussed above. Workers in the ^ 

field have estimated this pump heating effect to be anywhere 
between 0.5 K and 30 K above the ambient temperature. 

Careful system noise temperature measurements as a function ^ 

of bias level and upconverter pump power level may shed 
some light on this important question of pump diode heating. : 

Upconverter diode gain (see Table 1, note f)is approxi- 
mately proportional to the ratio of output frequency to input 
frequency. For the S-band to K-band upconverter, this ratio is 
high, and therefore noise contributions from upconverter 
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output circuit losses, interstage losses, and maser input noise 
temperature are small. For this reason, the system input noise 
temperature is, in fact, lower than the 4.7-K maser input noise 
temperature. Consider similar varactor diodes used in a 
S.S^Hz to 24-GHz upconverter. Gq for this case would be 
approximately 3 dB and system input noise contributions 
from £ 5 , and Tm would be substantially larger. If maxi* 
mum performance is to be obtained from upconverter/maser 
systems at X-band or Ky*band, careful effort will be needed to 
mininuze the loss of each component in the system, and the 
benefits of using a higher frequency maser postamplifier 
should be explored. 


VI. Noise Temperature Measurement s 

Noise temperature measurements of the upconverter- 
maser system were perfonned at JPL. A feed horn was 
attached to the input waveguide flange and a large piece of 
microwave absorber was used as an ambient temperature 
termination for the horn (see Fig. 9). System noise tempera- 
ture data are obtained by covering the feed horn with the 
absorber and then removing the absorber, allowing the horn to 
view the cold sky. Total operating system noise temperature is 
obtained by making precision power-level measurements with 
and without the absorber in place. A detailed description and 
analysis of noise temperature calibrations using ambient 
terminations has been published by C. T. Stelzried. (Ref. 10) 

Initial measurements to date have resulted in a total system 
operating noise temperature of 9.4 K at 2295 MHz and 9.5 K 
at 2388 MHz. Best estimates of noise contributions for the 
parts of the system are given in Table 2. The noise contribu- 
tion estimates for the sky and horn are those used in previous 
S-band maser noise calibrations (Ref. 8 ). The follow-up 
receiver contribution is measured by making a system power 
measurement with the maser pump source on and off. When 
these three contributions are subtracted from the total sy tem 
noise temperature, the effective input noise temperature 
for the upconverter/maser system is found to be 3.1 K at 


2295 MHz and 3.2 K at 2388 MHz. These values faB within 
the theoretical estimate of 2.5 to 3.5 K obtained in Table 1 . 

The upconverter noise temperature can be estimated by 
subtracting component noi«e contribution estimates (Table 1 ) 
from the system effective input noise temperatures reported 
above. This calculation results in an input noise temperature 
for the upconverter of less than 2 K. 

VII. Reliability 

The ability of the parametric upconverter to withstand 
repeated temperature cycling from room temperature to 4.5 K 
and to operate reliably at this temperature for extended 
periods of time is of importance. If a cryogenic component is 
to become an operational reality in the DSN, it must demon- 
strate high performance in the above areas as well as in areas of 
RF performarice. To date, the upconverter itself has logged 
over 1500 hours at 4.5 K and has undergone 5 full tempera- 
ture cycles of 300 K - 4.5 K - 300 K without degradation. 

VIII. Conclusions 

The combination of an experimental cryogenically-cooled 
S*band upconverter and a K-band maser postamplifier has 
demonstrated an effective input noise temperature of 3.1 K 
at 2295 MHz and 3.2 K at 2388 MHz. These initial results 
compare favorably with input noise temperatures of 2.0 to 
2.1 K achieved with the best S-band maser amplifiers (Ref, 8 ). 
These figures represent a remarkable achievement in noise 
performance for a system which has a potential instantaneous 
bandwidth of 400 to 500 MHz at S-band frequencies. 

Testing of this system will be continued. Feed horn con- 
figurations that are suitable for other frequencies within the 
2.0 to 2.5-GHz range will be obtained so that system noise 
temperature can be >.v.-iuated across the entire upconverter 
bandwidth. Other R1 parameters will be investigated and 
reliability in the cryogenic environment will be monitored. 
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Tabto 1. TtiMTOtical noiat liin5ianrture aatiiiMte. 2.3 GHz 


System schematic diagram 


Component 

description 


Component 
loss I, or 
gai T, or 
input noise 
temperature T 


Component contribution 
to effective input noise 
Value of temperature at S-band 

L^G^oiT waveguide input port 

Contnbution®**’ Va 


INPUT PORT*- 


S-band input 0.05 dB^ 

transmission line 


s r 


1 


0.141 semirigid Z .2 

cable 


Upconverter 
input circuit loss 


Varactor diodes 


Upconverter 
output circuit loss 


Interstage loss, 
including pump 
rejection filter 


0.03 dB** ^'1 


0.4 dB* Ijlj 0.44 K 


Gjj = 9 dB^ 
rD = 0.3K8-h 
Tq = 1.2 K"'* 


0.6 dB« 


0.9 dB** 


L,L.L.(L.-\) 

C Tccr 0.09 K 


V — ■Tccr 01’ k 


^1^2^3^4^5 


Total system effective input noise temperature 


2.5 K* 

3.5 K‘ 




Tabto 1. Theoretical noise temperature estimate, 3 = 2.3 GHz (contd) 


*CCR physical temperature * 4.5 K. 

^In the equations, losses L and gains G are dimensionless ratios. 


^Temperature along this line varies from 300 K to 4.5 K. Q measurements, insertion loss, a.id temperature gradient calculations were used to 
determine 0.5 K contribution. 


^Estimated from published values, previous measurements of similar components, and measured variations of insertion loss vs temperature. 

^See Section V. 

^For minimum noise tuning: 

Go (Ref. 10) 

^ M 

where 

Gjy - conversion gain of varactor diodes 

M ~ varactor diode figure of merit m|f^ 

= 80 GHz (Ref. 3) 

= output frequency 
= 24 GHz 
= input frequency 
= 2.3 GHz 

Substituting: 

Gj^ = 8.0 = 9 dB 

*If diode junction temperature is 4.5 K, i.e., no pump heating. 

^For mifiimum noise tuning: 

24 

T’d ^ (Ref. 10) 

O.aKforTy- ,.5K 
« UKforr^=20K 

where 


Tq = effective input noise temperature of /aractor diodes due to series resistance. 
T'j ' physical temperature of diode junction. 

‘If diode junction temperature is 20 K. 

^Calculation of maser effective input noise temperature gives * 4.7 K ± 0.5 K (Ref. 2) 
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imOm dm n o twm conuwvuotm fOr In© upcoffvvfiwr / nfWMr sytvvfn 


Part of «y$tcm 

Noise contribution, K 

2295 MHz 

2388 MHz 

Sky (includes atmosphere 
and cosmic bickground) 

4.9 

4.9 

Horn (includes mode 
generator and transition) 

1.2 

1.2 

Upconverter/Maser System 

3.1 

3.2 

Follow-up receiver 

0.2 

0.2 

Total operating system 
noise temperature T^p 

9.4 

9.5 



r 


CLOSEO-OrCLE HELIUM KEFR1GERA70R 4.5 K 


n 



K-BAND 
OUTPUT 
23.76 TO 
24.26 GHt 


Fig. 1. Block diagram of upconverter/maaar/CCR system 
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Fig. 2. Packaged upconvartar/matar/CCR system 
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Fig. 3. S-band parametric upconverter 
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Fig. 4. K-band maser structure 
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Fig. 5. K-band maser circulator and ruby waveguide detail 
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Fig. 6. K-band ma«er with superconducting magn«t 
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Simultaneous Dual-Frequency, Round-Trip Calibration of 
Doppler Data With Application to 
Radio Science Experiments 

A. L, Berman 
TDA Engineering Office 


Simultaneous dual- frequency, round-trip (uplink and downlink) calibration of Doppler 
data is expected to be a requirement of several radio science experiments being planned 
for the next decade, and such (calibration) capability is expected to be achie^^ed by the 
mid-1980s. Simultaneous dual- frequency , round-trip calibration would be straightforward 
except for the condition of unequal spacecraft turnaround ratios at S- and X-band. This 
article discusses the impact of unequal turnaround ratios on calibration accuracies in the 
specific cases of the Gravitational Wave Detection Experiment and the Solar Gravitational 


Quadrupolc Moment Experiment. 

I. introduction 

B> the early to inid-1980^, the Deep Space Network (DSN) 
IS expected to ha\e augmented the current S-band uplink 
capabihtv with X-band uplink capability. With additional 
niodificaiion to the spacecraft and ground tracking system, 
one can readily foresee the provision of the capability for 
simultaneous round-trip (uplink and downlink) transmission of 
S- and X-band frequencies. Simultaneous (and independent) 
transmission of S- and X-band frequency provides the impor- 
tant capability of high-precision calibration of charged particle 
effects on ladio metric (i.e.. Doppler and range) data. This new 
capability will be crucial to certain radio science experiments, 
and may ,id utionally prove to be beneficial in regard to navi- 
gational usage of radio metric data. Such round-trip, dual- 
frequency calibrations would be (ch c to) exact and trivial to 
apply, except for the fact that the spacecraft turnaround 
ratios* at S- and X-band are, according to the current design, 
expected to be (slightly) unequal. 


Simultaneous dual-frequency^ calibration of radio metric 


data is expected to be irr*’^^ 
experiments: 

"tant to the following radio science 


Sun-Earth- 

Radio metric data 

Radio science experiment 

probe angle, deg 

type of interest 

Gravitational Wave Detection 

- 180 

Diippler 

Relativity Tests 

- 0 

Range 

Solai Quadrupolc Measurement ~ 0 

Doppler 


Callalian (Ref. 1 ) has treated simultaneous dual-frequency 
calibration of range data. Although tlie mechanics of the 
calibration are similar for both the Doppler and range applica- 


* Ratio of spacecraft transmitted frequency to spacecraft received 
frequency for a single-frequency band. 

“Simultaneous dual-frequency" will specifically imply “ iniulianeous. 
independent, round-trip, dual-frequency." t 
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tions, the final efficacy of the calibration in any given applica- 
tion is highly dependent upon the particular circumstances of 
the application. 

This article will derive and discuss simultaneous dual- 
frequency calibration of Doppler data. Specifically empha- 
sized will be the application of simultaneous dual-frequency 
calibration of Doppler data in regard to the Gravitational Wave 
Detection Experiment (Ref. 2) and the Solar Quadrupole 
Moment Experiment (Ref. 3). 

II. Simultaneous Dual-Frequency, 

Round-Trip Calibration of Doppler Data 

It was noted in Section I that the n;.'iin difficulty with 
simultaneous dual-frequency calibration capability is that the 
spacecraft turnaround ratios at S- and X-band are expected to 
be (slightly) unequal. Figure 1 presents a simplified block dia- 
gram of dual-frequency transmission; the relevant parameters 
are defined as follows: 

= S-band uplink frequency 

^ 2,1 15 GHz 

= ratio of uplink X-band to uplink S-band 
^ 3.404 

= S-band spacecraft turnaround ratio 
1.086 

Cj = X-band spacecraft turnaround ratio 
^ 1.169 

One now makes the following additional definitions (where A0 
relates to columnar electron (number) density fluctuation): 

Uplink charged particle fluctuation effect (in cycles) 



Downlink charged particie fluctuation effect (in cycles) 



/ = frequency 

To isolate the frequency-dependent effects (i.e., charged parti- 
cle fluctuation effects), one scales the downlink X-band by the 
downlink S-X ratio {Cq/C^Kq) and differences this with the 


downlink S-band. The quantity being operated on is integrated 
(counted) Doppler frequency during a given (averaging time) 
interval: 

A0^ = measured charged particle fluctuation 
effect, in cycles (at S-band downlink 
frequency) 



where: 

= (integrated) S-band downlink frequency 
F^ = (integrated) X-band downlink frequency 

In tenns of A0^p and one expects for the measured 

charged particle fluctuation effect: 



By defining the difference between uplink and downlink 
charged particle fluctuation effect: 

e = A0 - A0 . 

^up ^dn 

one obtains the downlink charged particle fluctuation effect in 
terms of the measured charged particle fluctuation effect and 
the uplink-downlink difference: 



49 


One now writes the expression for the measured charged particle fluctuation effect in terms of the desired (scaled) calibration 
(first bracketed term) and a term that is dependent on the (unknown) uplink-downlink differences (second bracketed term): 



The appropriate (S-band) calibration is now identified as follows; 



Christensen (Ref. 4) first noted that dual-frequency calibration of either frequency is (essentially) equivalent; the X-band calibra- 
tion corresponding to the above (S-band) calibration is: 



These expressions reduce to: 

cl A<t>^^ a + 6.40 X 10’^ e (S-band) 

^1 ^ ^ e (X-band) 

The S-band result is similar to that obtained by Callahan 
(Ref. 1). To utilize these relationships, one needs to be able to 
relate e to this can only be done in the context of 
specific experimental conditions. Applications of these results 
to specific radio science experiments will be discussed in 
Sections III and IV to foUow. 


III. Gravitational Wave Detection 
Experiment Application 

The gravitational wave detection experiment utilizing ultra- 
precise Doppler data has been described in detail in Ref. 2. 
The experiment is expected to be performed under the 
following conditions: 

(1) SEP 180 deg (specifically to reduce solar wind 
fluctuations). 

(2) Long round-trip light times (RTLTs). 

Under such conditions, most of the electron density (and 
density fluctuation) is near-Earth, so that the correlation time 
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of interest between uplink and downlink solar wind fluctua- 
tion is essentially the RTLT. The Doppler sample interval of 
interest will be significantly less than the RTLT (Ref. 2). 
Under these circumstances, one can expect little correlation 
between and A reasonable and quite conservative 
assumption to make under such conditions is that the 
ifference between A0^p and is of the same order as the 
quantities themselves, or, very approximately: 

€ - Ad - A<t>. ^ A<t> 

^up ^dn 

For the X-band calibration, one thus expects an error of 
approximately 0.1 (^.17 X 10“^ simul- 

taneous dual-frequency calibration with the currently 
expected unequal turnaround ratios. Since the quantity of 
interest to the gravitational wave detection experiment is 
fractional frequency fluctuation: 

g{AF/F) 

one expects a sample order-of-magnitude improvement over 
the performance for uncalibrated two-way X-band (uplink and 
downlink). Figures 2 and 3 update previous graphs (Ref. 5) of 
(with = averaging time): 


a(AF/F) versus 

under the assumption that simultaneous dual-frequency cali- 
brations result in an order-of-magnitude improvement for the 
gravitational wave detection experiment conditions. 


IV. Solar Quadrupole Moment Experiment 
Application 

The experiment to measure the solar gravitational quad- 
rapoie moment via a spacecraft passing very close to the Sun 
(for instance, at four solar radii as in the proposed Solar Probe 
Mission) ?s described by Anderson, et al., in Ref. 3. In this 
experiment, Doppler measurements at an approximate 
10-minute sample interval are the quantity of interest. Most 
pertinent to the consideration of simultaneous dual-frequency 
calibration is that the spacecraft is (essentially) at the closest 
approach distance during the time of minimum SEP (^1 deg 
for the Solar Probe Mission), hence the spacecraft is located at 
the point of maximum signal path electron density (and 
density fluctuation). The importance of this is that for any 
Doppler averaging times greater than \ few seconds the 


and A^^j, will be extremely well correlated. Allowing the very 
conservative approximation (for a 10-minute averaging time): 

results in an error of approximately 0.01 (~9.17X 10"^ 
e/A0^), hence one can expect a two orders-of-magnitude 
improvement via simultaneous dual-frequency calibration of 
X-band as compared to uncalibrated two-way X-band. From 
Ref. 3, the uncalibrated two-way X-band error at four solar 
radii is (solar wind fluctuation error only): 

la - 2.5mm/s 

and therefore, with simultaneous dual-frequency calibration of 
X-hand, one can expect: 

la = 2.5 X 10“^ mm/s 

At the current time, the ground tracking system (rubidium 
frequency standard) has a frequency fluctuation of approxi- 
mately (with = 10 minutes): 

la ^ 1.0 X 10“^ mm/s 

so that even with unequal spacecraft turnaround ratios, 
simultaneous dual-frequency calibrations appear to he able to 
reduce solar wind fluctuation in this experiment to below the 
current (system) inherent limitations. 

V. Discussion and Summary 

Simultaneous dual-frequency, round-trip (uplink and down- 
link) capability provides the crucial ability to be able to 
calibrate solar wind fluctuation in radio science experiments 
that require highly precise Doppler measurements. The calibra- 
tion process is complicated by the usage (.as per current design) 
of unequal spacecraft turnaround ratios at S- and X-band. In 
the case of the Gravitational Wave Detection Experiment, 
the unequal turnaround ratios limit the calibration to an order- 
of-magnitude improvement over the uncalibrated two-way 
(uplink and downlink) X-band; eventually this may result in 
solar wind fluctuation being the limiting factor in achieving 
increased system sensitivity. For the Solar Cr'^vitational 
Quadrupole Moment Experiment, on the other i.^uid, the 
unequal turnaround ratios will basically not compromise the 
calibration accuracy of simultaneous dual-frequency, round- 
trip capability. 
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Fig. 1. Simplified functional block diagram of aimuttaneous dual* 
frequency upiink and downiink tranamission (Fq 2.1 IS GHz, 
Ko ^ 3.404, Co 1.086, C^ ^ 1.196) 
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A Solar Wind Turbulence Event During the Voyager 1978 
Solar Conjunction Profiled via a New DSN Radio Science 

Data Capability 

A. L. Berman and A. D. Conteas 
TOA Engineering Office 


The Deep Space Network (DSN) has implemented a new radio science data capability 
within the DSN Tracking System - routine provision of phase fluctuation data concur- 
rently computea over several different time scales. This new capability has been used to 
observe phase fluctuation spectral characteristics during a rapid increase in solar wind 
iurbulence that occurred during a July 23, 1978, track of ihe Voyager 1 spacecraft by 
Deep Space Station (DSS^ 1 1, This article suggests that the new capability will prove 
quite useful in future studies of variations in Solar Wind phase fluctuation spectral 
characteristics with, for instance, parameters such as the solar (sunspot) cycle and radial 
distance. 


I. Introduction 

During 1977, the Deep Space Network (DSN) iinplemeiited 
a new radio science data capability within the DSN Tracking 
System. This new capability was described in detail in Berman, 
Ref. 1 . Esseniially, the Tracking System routinely (i.e„ auto- 
matically whenever a spacecraft is being tracked) computes 
rms phase fluctuation (0) over the fluctuation frequency (v) 
range (where v is related to the Doppler averaging time (r^) by 
the approximate relationship u ^ (36Tj)“ * ): 

2.8 X 10’^ Hz<j;<2.8X 10“' Hz 

The method of computation is briefly summarized as follows. 
Doppler samples (accumulated phase in cycles) are provided 


concurrently at four sample intervals: 0.1 s, 1,0 s, 10.0 s, and 
100.0 s. In each case, 18 Doppler frequency samples (differ- 
enced Doppler accumulated phase divided by averaging time) 
are differenced with the equivalent predicted Doppler fre- 
quency. These 18 differenced (actual minus predicted) fre- 
quency samples are then fit with a least squares linear curve fit 
to remove low-frequency trajectory errors (and in the process, 
low-frequency solar wind induced phase fluctuations). The 
(scaled) rms phase fluctuation where rms 

frequency fluctuation, is then computed relative to the least 
squares linear curve fit. Subsequent to the fit, 100 contiguous 
rms phase fluctuation computations for the 0.1 -s Doppler 
sample interval data, and 10 contiguous rms phase fluctuation 
computations for the 1.0-s Doppler sample interval data are 
averaged The final provision of data is summarized in Table 1 . 
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This article presents and briefly describes phase fluctuation !i!. ColUfTindr Soldf Wind PhaSG Fluctuation 
data generated during a solar wind turbulence “event” so as to Spectral Index 

illustrate the potentiality of this new DSN capability. 

The relationship between the columnar (two-dimensional) 
spectral index is simply (from Berman, Ref. 1, or Douglass, 

II. The Data 

On July 23, 1978, at about 21:00 Greenwich Mean Time 
(GMT), the Voyager 1 (ground) observed two-way Doppler 
data generated by Deep Space Station (DSS) 1 1 began to 
indicate that a solar wind flow of (relatively ) high turbulence 
was interdicting the Voyager 1 signal path. At the time. 

Voyager 1 was at a Sun-Earth-probe angle (SEP) of approxi- 
mately 6.4 deg. Figure 1 presents the phase fluctuation data 
for this period. The data presented in Fig. 1 were computed 
from two-way Doppler data generated by DSS 1 1 . Tlie vertical 
scale in Fig. 1 is Doppler ims phase fluctuation (0, deg), wliich 
is related to rms frequency fluctuation (o^) by the approxi- 
mate relationship (Ref. 1 ): 

0 ( 0 I 

For the 0.1-, 1.0-, and 10.0-s Doppler sample interval calcula- 
tions in Fig. 1 , ten 3-min phase fluctuation computations (as 
described in Section I) have been averaged to produce each 
30-min (interval) data point. This explains the considerable 
“smoothness” of the 0.1-, 1.0, and 10.0-s Doppler sample 
interval data as compared to the 100.0-s Doppler sample 
interval data. As can readily be seen, the behaviour of the 
phase fluctuation data in all Doppler sample interval regimes is 
quite similar. The increase in the 0.1 -s Doppler sample interval 
phase fluctuation data is relatively less than the longer Doppler 
sample interval data because the 0.1 -s Doppler sample interval 
data contain phase fluctuation components of both the solar 
wind and the ground tracking system. For the 1.0-, 10.0-, and 
100.0-s Doppler sample intervals, the phase fluctuation data 
are essentially all a result of Solar Wind turbulence. 


References 

1. Berman, A. L„ “Phase Fluctuation Spectra: New Radio Science Information To 
Become Available in the DSN Tracking System Mark III-77”, in 77/e Deep Space 
Network Progress Report ^2-40, Jet Propulsion Laboratory, Pasadena. California, 
August 15, 1977. 

2. Douglass, D. H., “Detection of Gravitational Radiation and Oscillations of the Sun via 
Doppler Tracking of Spacecraft”, in A Close-up of the Sun, edited by Neugebauer, M., 
and Davies, R. W., JPL Publication 78-70, Jet Propulsion Laboratory Pasadena, 
California, September 1, 1978. 

56 



2.8 X 10' - Mz<t/<2.8X lO"^ Hz 
2.8 X lO'"* Hz < 1 / < 2.8 X 10’^ Hz 

IV. Discussion and Summary 

The DSN has implemented a new radio science data capabil- 
ity within the DSN Tracking System - routine provision of 
phase fluctuation data concurrently over several different time 
scales. The columnar solar wind phase fluctuatioii spectral 
index can be easily extracted from such data, and hence tliis 
new capability should prove quite valuable in studies such as 
the radial and solar cycle variations of the spectral index. 


where 

0 = RMS phase fluctuation 
V = tluctuaiion frequency 
Kq = spectral index 

Figure 1 additionally provides spectral index data computed 
from the phase fluctuation data. Inherent to these data is the 
assumption that the fluctuation spectrum is power-law with 
frequency. Spectral index data at 30-min intervals were com- 
puted for the following frequency ranges: 



TaM 1. Proviston of phase fluctuation data 


Doppler 
sample 
interval, s 

Number of 
Doppler 
samples 

Number of 
estimates 
averaged 

Print-out 

•nterval, 

min 

0.1 

18 

100 

3 

1.0 

18 

10 

3 

10.0 

18 

1 

3 

100.0 

18 

1 

30 
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Radial and Solar Cycle Variations in the Solar Wind Phase 
Fluctuation Spectral Index as Determined From Voyager 

1978 Solar Conjunction Data 

A. L. Berman and A. D. Conteas 
TDA Engineering Office 


Of current interest is the value of and possible variations in the solar wind phase flue- 
tuation spectral index. This article presents columnar spectral index information that has 
been extracted from a sizable volume of Voyager 1978 solar conjunction doppler phase 
fluctuation data. The Voyager 1978 results, when compared to similar information 
derived from the 1976 Helios and Viking Solar Conjunctions, lead to the following 
inferences: (1) there has been a significant change in the spectral index from 1976 to 
1978; 12) there is continuing evidence that favors a slight (positive) correlation between 
the spectral index and the solar (sunspot) cycle; (3) there is little or no evidence in sup- 
port of a radial variation of the spectral index. 


I. Introduction 

During the last decade, much work has been performed to 
determine the form of the solar wind fluctuation spectrum. 
During this period, and through analysis of both spacecraft 
in*situ and radio scattering measurements, it has become 
evident that tJ'e form of the fluctuation spectrum is best 
described as power-law with fluctuation frequency (»/); see, 
for example, Intriligator, et al., Ref. 1; Goldstein, et al., 
Ref. 2; Unti, et al., Ref. 3; Woo, et al.. Ref. 4; and Chang, 
Ref. 5. It is currently of considerable interest to determine the 
numerical value of the (power-law) spectral index*, and. 


^The tenn “spectral index'’ wiU specifically imply the columnar (two- 
dimensional) spectral index. 


additionally, any possible functional dependencies of the 
spectral index on: 

(1) Fluctuation frequency regime. 

(2) Radial distance. 

(3) Solar cycle variations. 

Recently, Coles, et al., Ref. 6, have determined spectral 
index variations with both radial distance and fluctuation 
frequency regime through analysis of intensity scintillation of 
incoheicnt radio sources. In contrast to these results, Berman, 
(Ref. 7). analyzing Viking Doppler phase fluctuation data 
generated b\ ihe Deep Space Network (DSN) during the 1976 
solar conjuncuon, found no indication of a radial dependence 
of the spectral index. It should be noted, however, that these 
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(differing) results were obtained in different regimes of the 
fluctuation spectrum; whereas the scintillation results from 
Ref. 6 are for fluctuation frequencies (v) greater than 0.1 Hz, 
the Doppler results are for fluctuation frequencies less than 
0.03 Hz. 

To date, little work has been performed on determining 
possible variations of the spectral index with solar cycle. 
Berman recently noted (Ref. 8) that, based on the totality of 
spectral index measurements to date (1977), the (available) 
evidence pointed to, at most, only a slight variation of spec- 
tral index with solar cycle. Further, Ref. 8 indicated that if 
indeed correlation did exist between the spectral index and 
solar cycle, the correlation appeared to be positive, i.e., the 
fluctuation spectrum appeared to be steeper during solar 
cycle (sunspot) maximum. 

This article presents the results of extracting spectral index 
information from a large volume of Doppler phase fluctuation 
data generated by the DSN during the Voyager 1978 solar 
conjunction. Analysis of these data indicates: 

(1) There vvas a significant increase in the spectral index 
between solar cycle minimum (1976) and (near) solar 
cycle mjiximum (1978). 

(2) There is .ontinued evidence for slight (positive) corre- 
lation between the spectral index and the Solar Cycle. 

(3) There is little evidence of a significant variation of the 
spectral index with radial distance (at least in the low 
frequency portion of the fluctuation spectrum). 

II. Helios and Viking Spectral Index Data 
During 1976 

Berman (Refs, 7 and 9) has analyzed a large volume of 
Doppler phase fluctuation data generated by the DSN during 
the 1976 solar conjunctions of Helios and Viking. Analysis 
of the spectral index information aerived from the Helios and 
Viking data yielded the following conclusions (Ref. 7): 

(1) The mean value of the spectral index during 1976 
was 2.42. 

(2) There was little or no indication of a spectral index 
radial dependence. 

Figure 1 (from Ref. 7) presents the spectral index, as com- 
puted from Viking Doppler phase fluctuation data, versus 
Sun-Earth-probe (SEP) angle. The absence of any significant 
spectral index radial dependence is obvious from inspection 
of Fig. 1. 


III. Voyager Spectral Index Data During 1978 

In 1977, the DSN implemented a new capability to allow 
the convenient extraction of spectral index information from 
Doppler phase fluctuation data; this capability is described in 
detail in Refs. 7 and 10- Using this new capability, spectral 
index information has been computed from Voyager 1 Dop- 
pler phase fluctuation data. Very briefly, the spectral index 
extraction process is predicated on the relationship (Ref. 10): 


where: 

<t> = rms Ooppler phase fluctuation 
V = fluctuation frequency 
= spectral index 


The Voyager 1 two-way Doppler phase fluctuation data 
were analyzed during the periods June 17, 1978, to July 7, 
1978, and July 20, 1978, to Augusts, 1978. During these 
periods, the SEP varied between 3.6 deg and 17.9 deg. The 
fluctuation frequency (v) range over which the spectral index 
was computed for the Voyager 1 data was: 


2.8X 10'^ Hz<j^<2.8X 10“^ Hz 


The spectral index information determined from the Voy- 
ager 1 Doppler phase fluctuation data is presented in Figs. 2 
and 3. Figure 2 presents the spectral index vs day of year 
(DOY), while Fig. 3 presents the spectral index vs SEP. 

The mean value of the spectral index as computed from 
the Voyager 1 data is = 2.67, as compared to a spectral 
index mean value of Kq = 2.42 as determined from the Viking 
and Helios data. This difference is considered statistically 
significant in light of the large volume of data analyzed in 
both cases. For the 1976 mean spectral index, 42 pass average 
(*^3 hours of Doppler data) spectral index measurements were 
available, while in the case of the 1978 mean spectral index, 
71 pass average measurements were utilized. Since 1976 
constituted solar cycle 20 minimum, and mid-1978 is near 
(according to Sunspot number) solar cycle 21 maximum, 
these results can be considered perhaps the strongest evidence 
to date of (positive) correlation between spectral index and 
solar (sunspot) cycle. 
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^ Figures presents the spectral index vs SEP. Although there (1) There has been a significant change in the spectral 

appears to be a slight negative slope to the data, it is con- index from solar cycle minimum (1976; Aq = 2.42) to 

sidered that this is not significant, particularly when con- (near) solar cycle maximum (1978; Aq = 2.67). 

sidered in context of the (smaller) SEP data of Fig. 1 , which 

in no way shows any continuation of such a trend (i.e., a (2) There continues to be evidence for a slight (positive) 

spectral index (negative) correlation with SEP). When con- correlation between spectral index and solar (sunspot) 

sidered jointly, it is suggested that Figs. 1 and 3 provide little cycle, 

indication of any significant spectral index radial variation. 

(3) There continues to be little or no evidence for a signi- 
ficant variation of spectral index with radial distance. 

IV. Summary and Discussion 

A new DSN capability has been utilized to obtain a sizeable This new DSN capability is expected to be exercised in 
volume of spectral index information during the Voyager 1978 future years to derive spectral index information from Voy- 

solar conjunction. Major conclusions derived from a compara- ager solar conjunction Doppler phase fluct-.?.tion da*.a, and 

I tive study of these data with similar data obtained during 1976 hence to allow the continued monitoring of spectral index 

solar conjunctions of Helios and Viking are: variations during the remainder of solar cycle 21 . 
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Fig. 2. Columnar solar wind phase fluctuation spectral index vs day 









N79 


16012 




DSN Progress Report 42-48 


September arxJ October 1 978 


A Simplified Algorithm for Correcting Both Errors and 

Erasures of R-S Codes 

I. S. Reed 

University of Southern California 

T. K. Truong 
TDA Engineering Office 


Using the finite field transform and continued fractions, a simplified algorithm for 
decoding Reed-Solomon codes is developed to correct erasi^res caused by other codes as 
well as errors over the finite field GF(q^ }, where q is a prime and m is an integer. Such 
an R-S decoder can be faster and simpler than a decoder dutt uses more conventional 
methods. 


I. Introduction 

Fast real-valued transforms over the group (ZjX" were used 
first by Green (Ref. 1 ) to decode the (32, 6) Reed-MuUer code 
(Ref. 2) used by the Jet Propulsion Laboratory (JPL) in 
Mariner and Viking space probes. Gore (Ref. 3) extended 
Mandelbautn’s methods (Ref. 4) for (Recoding R-S codes 
(Ref. 5). He proposed to decode R-S codes with a finite field 
transform over Gf(2'"), where m is an inte^r. Miche'-on 
(Ref. 6) implemented Mandelbaum’s algorithm and showed 
that the decoder, using the transform over GFil”*), requires 
substantially fewer multiplications than the standard decoder 
(Ref. 7). Recently, it was shown (Ref. 8) that R-S cod'.s can 
be decoded efficiently with a combination of the fi»',t trans- 
form method and continued fractions. 

For a space communication link, it was shown in (Ref. 9) 
that an R-S code that is conc-itenated with a Viterbi-decoded 
convolutional code can be us id to reduce the signal-to-noise 
ratio required to meet a :?'<'ified bit error rate. Such a 


concatenated R-S code is implemented presently on the Voy- 
ager spacecraft. 

In such a concatenated code, the inner convolutional de- 
coder is sometimes able to find only two or more equally 
probable R-S symbols. Then it is possible to declare an erasure 
of the R-S symbol, i.e., mark the location of a possible error. 
The outer R-S code can take advantage of the knowledge of 
these erasure locations if erasures are reported to the R-S 
decoder by the Viterbi decoder. Erasures are not reported by 
the MCD as presently implemented in the DSN. 

In this article, a simplified decoding algorithm is developed 
to correct erasures and errors of R-S cod.i using both a finite 
transform and continued fractions. This decoding algorithm is 
based on the algorithm originally invented by Forney (Refs. 10 
and 11). An important ajvantage of the present decoding 
techmque over previous metho 1» is that a Chien-type search 
(Refs. 7 and 1 11 for the roots of thr er’or locator polynomial 
is completely avoided. 
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In general, let 


II. On the Decoding of Erasures and Errors 

Let n be the block length of an R-S code in GF ((?'"). Also 
let d be the minimum distance of the code where d- P- 1. 

Then n = I, where P is the number of parity symbols and / Y.X^ + ^ for all /c (lb) 

is the number of information symbols. /= > /- * 


Define the following five different vectors: 

(Cfl, Cl , j ) = c, code vector 

(^o» ''i ’ “ *» '‘n- 1 ) ^ received vector 
(Mq, M 1 , 1 ) = u, erasure vector 

(Cq, Cl , • *, 1 ) = c. error vector 

(u^, u 1 , * • *, u^_ I ) = IT, new erasure vector 
These vectors are related by u = e + u and r = c + u + e. 

Suppose that t errors and s erasures occur in the received 
vector r and tliat s^ 2t<d - 1 . The syndromes are given by 

s,=K«*)= E 

i =0 

* E for 1 < A: < d - 1 

1=0 i =0 1=0 


>^cre Ty is the /-th error amplitude, is the /-th error 
location, Wj is the /-th erasure amplitude, is the y-th known 
erasure location, and = 5^ is known for . <k<d- 1 . 

Following Forney, let the erasure locator polynomial be 
defined by 


r(x)=n(*-V = E(-l)V' 

/=! /=0 

This implies 

s 

(-l)^r.Z*-r = 0 forJfc=l,2. -.J (3) 

;=o 

where Tq = 1 and t^’s are known functions of Z| , Zj, * * *, Z, 
for 1 </ < j. The Forney syndrome (Ref. 10) is defined by 


= E for 1 < I <d - 1 - I (4a) 

/=0 


but 

n- I 

2 fori <*<</-! 

MO 

Thus 


From Eqs. (la) and (3), we know that fo' / = 1 , 2, — *, s and 
5;^ for A: * 1, 2, * 1 is known. Thus the F/s for 1 <i < 

d- 1 ’ s are known. In general, let 


T,‘ £ 

/=0 


fot all i 


(4b) 


= E E 


where 1 i + s - / < d - 1. Then the 

Foiney s^-ndrome is known for I < i - 1 - s. 


■E'r; forl<^<£/ 1 

/^i /'I 

(la) 


Now if one substitutes Eq. (lb) into Fq. (4b), 

I 

r* ^ (•'D'r E 
;=0 


4 
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(7) 


/=o L*=i *=* J 


t •' 


■ E E •v')' 


lc=l /=! 


-E E <-»''/ ■*r' 


k=l ;=0 


A 


= V z) 


+ E ^'k^k E 

;=0 




From Eq. (3) one observes 


where 


<Kjc) = 11^^“ ^k) 


k=l 


E K^i E 

Jk=0 /=0 

Thus, 


is the error locator polynomial and degp(x) < deg a(x). Since 
is known for/= 1,2, • • • , J - 1 - 5 , where r <</ - 1, then, 
by Theorem 2 in Ref. 8, the error locator o(x) can be obtained 
by using continued fractions. The roots of a(x) arc the loca« 
tions of the t errors. 


7;"E r,ylZ<r‘^’i’yr' 

k^l /=0 

1 

*=1 

^ere the quantities 

^■^E 

/=o 

for « 1, 2, • • •, r ue noi known. 

Now let TQc) be a formal power aeries in the indeterminate 
X, defined by 

r(x) = r,x" * + 2'jX~* + TjX"’ + • • • 

Then substituting Eq. (S) into Eq. (6), gives 

r(x) = f; v' = 2 ( E 


Since the locations of ihe t errors and the s erasures are 
now known, we may assume that we have the problem of only 
s erasures where s = / + s. That is, the only unknowns are the 
“erasure” amplitudes •••, the amplitudes of 

both the error and erasure vectors. The corresponding known 
locators are 2^, • - *, Only the case of erasures need 
be considered. 

Suppose there are s' < J - 1 erasures. Then r = c + u, where 
u = u + e is a new erasure vector in which s' erasures occur 
^ere the received vector is r and c is the code vector. The 
syndromes are 


3^ = r(a*) = ^ c^a*' + “/*** 

/=o f=o /• 0 

forifc=l,2,.-,rf- 1 
but 

lf-1 

y ! = 0 for *= 1 , 2 , •••,</- 1 

Thus 

/»0 


I 



for*» 1. 2, I 




/=! 


(8a) 


where Wj are the /-th amplitudes of the erasure and error 
vectors and ^ are the /-th known locations of the original 
erasure and error vecton. Note that actually equals 5^ in 
Eq. (la)for k= 1, 2, • • 1. In genu j,let 

(SI’) 


where ^ ^ » 2, • • % - 1 . The erasure polynomial 

of the original erasure and enor vector is given by 


/=! itO 

^^diere 7^^ 1 , s' = r + s, and t(x) and o(x) arc given in Eqs. (2) 
and (7), respectively. Hence 


where f/j, are knovm. From Eq. (11) one 

obtains the rest of the transform of u, i.e,, £/g for 0 < £ < 
n- 1, llie vector of amplitudes u is found by taking the 
inverse transform over (7F(q"*) of i/g for 8 = 0, 1, 2, * * •, 
n - 1 . Finally, the original n-tuple code vector can be obtained 
by subtracting u from the received vector r. 

Let us recapitulate the decoding of R-S codes for both 
errors and erasures using the transform over GF(q^) and 
continued fractions. This algorithm is composed of the follow- 
ing five steps: 

(a) Compute the transform over GF (q”*) of the received 

vector «-tuple, t'n- 1 ). 

1=0 

where r^e GF(q^X a is an element of order n, and is 
the minimum distance of the R>S code. 

(b) Compute for / =» 0, 1, 2, • • s from the erasure 
locator polynomial 


r(Zf) = £ = 0 

*=o 


fori = 1 , 2, • • • , s' 


( 9 ) 


T(.x) = (X - -Z/) = Z (- 


/-I 


/=o 


Multiplying Eq. (9) by 

+ Z (10) 


where s is the number of erasures in the received 
vector. The Zj*s for 1 < ;<s are the known erasure 
locations. Then compute the Forney syndrome T„ for 
1 <n - 1 - s from the equation 


Summing Eq. (10)on/fori = 2,2, • * ss', yields 

/=! 1=0 *=I 

Hius 

a=i 


s 

y, = Z Ki<d-l-s 

/=0 

where for 1 < / <s and £) = for I </ - I arc 

known. 

(c) Use continued fractions to determine the error locator 
polynomial o(jf) from the knowrr T/s for 1 < / < d - 
1-5. 

(d) Compute the erasure and error locator polynomials 
from the equation 


In general, 


*•' Z (■!)*% for«>t/-l (11) 

*=l 


r(x) » a(x)r(x) « £ (-l)^?*Jf**'‘^ 
*=0 


where a(jr) and r(Jt) aiC the known polynomial. *^en 


6 ? 



compute the rest of the transform of the erasure and 
error vector from the C 4 uation 


r = y' (-iVr.S 


/-O 


2 




s*t 


u. 




*=1 


= 5 ^ - 2S +95 

w+2 n+1 n 


fori < n l-j = 2 


where £/g = 5g for 1 < 5 - 1 . ^2 ^ ^h® 

(e) Invert the transform to recover the error and erasure 

vector, then obtain the corrected code vector. T(x) = T x" * + ^ • — + . 

' 1 2 a— 1 —I 


To illustrate the decoding procedure for conecting errors 
with erasures, a simple example of an R-S code over GF (17) is 
now presented. 


= 2x“* +4of”^ + ?x'^ + 


D(x) 

o(x) 


Example: Let GF(17) be the field of integers, modulo the 
Fermat prime = 22^ + 1 = 17. We consider the correction 
of one error and two erasures in an eighttuple R-S code over 
GF(17). For this case n=8, ^/~1 = F=4, / = l,s=2, and 
d- l = 2r + 5. 

Let c = (5, 2, 12, 15, 2, 3, 2, 1) be transmitted. Assume the 
erasure vector is u = (0, 0, -3, 0, 0, -2, C, 0) *»nd the error 
vector is e = (0, 2, 0, 0, 0, 0, 0, 0). Then u = u + e = (0, 2, -3, 
0, 0, -2, 0, 0) is a new erasure vector. Hence the received 
vcctorisr®c + u + € = (5,4,9,15,2, 1,2, 1). 


Since ^ 1 - s = 2r = 2, by Theorem 2 in Ref. 8, o(x) can be 
determined by the use of continued fractions. Thus, 7T(x) = 
2/(x - 2). Hence o(jc) = x - 2. 

The erasure and error locator polynomial is 

r(x) = a(x)r(x) = (x - 2%x^ - 2x + 9) = x^ - 4x* + 13x- 1 

Thus, ?Q =!,?,= 4, ?2 = -4, and = 1. The rest of the 
transform for o is given by 


Now the syndromes for r are 




£'. = S,='<2*)=i;r,2*' 

#=0 




= 2(2*)' - 3(2*)^ - 2(2*)* 
for I < A: < </ - 1 = 4 

That is. (/, * S, a -4, £/, » 5 , = 3 , (/, * 5, * 10, and = 
* -3. The erasure locator polynomial is given by 


r(x) =f[ (x-2.) = (x-2*)(x-2*) = (x-4)(x-32) 

/=i 

= X* - 2x + 9 

Thus we obtain = 1, Tj =2, and - 9. Then the Forney 
syndromes are 


Thus, f/j = -3, = 3, {/^ = -3, and The 

inverse transform over GF(17) of the for 0< it <7 is 
given by 

E i/*2-'* = (-2) 2 t/,2-'* 

*=2 *-0 

for/ = 0. 1.2.--.7 


Hence, u * (0, 3, -3, 0. 0,-3, 0, 0). The corrected code vector 
is thus 

c = r-u * (5.4.9, l.S,3.1, 3,1)- (0,3, -3,0, 0,-3, 0,0) 
= <5.2. 13. 15.3.3,3 ') 
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This article describes the effects of the DSN software methodology, as implemented 
under the DSN Programming System, on the DSN Mark III Data Subsystems Implementa- 
tion Project (MDS), The software meinodology is found to provide a markedly increased 
visibility to management, and to produce software of greater reliability at a small decrease 
in implementation cost. It is also projected that additional savings will result during the 
mainterumce phase. Documentation support is identified as an area that is receiving fur- 
ther attention. 


CMF Communications Monitor 

and Formatter 

DST Host 

The subsystcios arc broken down by size in Table 1 . 


I. Introduction 

This anicle describes the quantitative and qualitative effects 
of the software methodology (Ref. 1) as implemented under 
the JSN Programming System, on the DSN Mark III Data 
Subsystems Implementation Project (MDS). The period in the 
life cycle of the MDS that was studied extended from Sep* 
tember 1974 through November 1976. During this period, 
approximately 100,000 lines of application code were 
delivered as well as diagnostic and system software. 

The scope of this article covers the implementation of the 
appioximately 100,000 lines of appheation code that can be 
directly allocated to the following subsystems: 


DCD 

Command 

DTM 

Telemetry 

DIK 

Tracking 

I»4C 

Mordtor 

DTT 

Test and Training 


The effects of the methodology on a single subsystem were 
studied in detail and the results, with corroboration from a 
second subsystem, were extrapolated for the entire effort. In 
addition to this quantitative evaluation of the MDS imple- 
mentation, a qualitative evaluation by cognizant software 
engineers is included at the end of the report. 


At the start of the MDS P oject most of the Standard Prac- 
tices governing software implementation had been formulate 
but weie as yet unpublished. However, the concepts goverr 
the methodology were well known and interim documt. 
were published to guide software implementation. The princi- 
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pal concepts used in the MDS software effort and implied by 
the term “methodology” in this report are: 

(1) Structured Programming. 

(2) Top down Devcic;, 

(3) Visibility Reporting. 

(4) Management Reviews. 

(5) Documentation. 

ill. Impact of Methodology on the 
Command Subsystem 

The Deep Space Station Command Subsystem (DCD) was 
cnosen as a representative system for the MDS for the follow- 
ing reasons: 

(1) The DCD software was of the “average” size of the 
MDS subsystem software. 

(2) The DCD was the best documented subsystem in terms 
of statistics gathered during the project lifetime under 
study. 

(3) The DCD development team most closely adhered to 
the methodology. 

Several important milestones should be noted that occurred 
during the Command subsystem implementation and repre- 
sented application of the methodology. These milestones 


were: 

The Software Requirements Document 

(SRD) published for DCD 3/3/75 

The Software Definition Document 

(SDD) published for DCD 9/19/78 

DCD High-Level Design Review 11/13/75 

MDS System Level Review 3/13/76 

Woik Breakdown Structure 

(WBS) Methodology applied 4/1/76 

Monthly Management Reporting started 5/1/76 

Weekly Reports 9/8/76 


The impact of these milestones will be discussed. 

A. Visibility 

Application of the methodology resulted in a series of mile- 
stones that were designed tr increase management visibiliiy. 
Visibility was indeed increased; an example is found in the 
management of the subsystem delivery date. The SRD 


reported an estimated delivery date of 1 August 1976. By the 
time the SDD was published, six months later, the estimated 
delivery date had slipped one month. Six months later, at the 
MDS System Level Review, a further slip of two months for 
the estimated delivery date had occurrec\ The successive slips 
of estimated delivery date can be projected (see Fig. 1). This 
projection is nonlinear and. ds uie desired delivery date is 
approached, the number of slips ten:*" to increase. This slip 
in estimated delivery date was typical for all subsystems, but 
is illustrated only for the DCD. 

It was possible bee use of the visibility to determine seven 
months in advance of the need date that the actual delivery 
date of the software was slipping. Increased visibility plus the 
application of the WBS methodology allowed the splitting of 
the effort into two phases. The software necessary to meet 
first priority requirements was identified and that portion of 
the software was delivered as a first phase, thereby pulling the 
delivery date back to the need date. Without corrective action, 
delivery would have taken place four to five months late. 

B. Anomalies 

Anomalies were not originally reported for all subsystems 
and original published anomaly repoits do not refiect the DTT 
or DST subsystems. For the MDS, approximately 200 anoma- 
lies were reported at verification time and during acceptance 
testing at CTA 21. After acceptance testing, approximately 
90 more were uncovered. 

For the DCD, 77 anomalies were reported during verifica- 
tion testing and acceptance testing, and 43 after acceptance 
testing. While the anomaly curve for the DCD presents the 
typical curve exhibited for the MDS (Figs. 2 ar.J ^), it can 
easily be determined that the rate is approximately twice that 
for the project as a whole. 

Possible causes for this phenomenon could be: 

(1) Tire software was not ready for delivery. 

(2) Reporting of preacceptance test anomalies was more 
comprehensive than for other subsystems. 

(3) The software was more difficult. 

(4) The software was considered more crucial to flight 
operations and therefore more ngorously tested. 

(5) The methodology docs not affect the number of 
anomalies. 

Close investigation shows that reasons (2) and (4) are 
probably the true reasons for this phenomenon. It has already 
been mentioned that the DCD was the best documented sub- 
system and team members adhered most closely to the meth- 



odology, and therefore were more meticulous about reporting 
errors during the verification phase. Exte.isive performance 
test programs were written that exercised the software 
extensively. 

Even at the rate of 4.8 errors per 1000 lines of code prior 
to delivery, the DCD software compares favorably to the 
industry average of 18 per 1000 hnes of code (see Refs. 2, 

4). 

C. Implementation Coat 

One measure of the effectiveness of a programming meth- 
odology is the productivity rate in delivered lines of source 
code per man-month (LOC/mm). An industrial avt;rage pro- 
ductivity rate is available for projects produced without the 
elements of the software methodology. For a comparison, it 
is desirable to determine the corresponding productivity rate 
for the Command Subsystem. 

It is felt that for a fair comparison, only equivalent activi- 
ties should be included in the productivity measure. Thus, 
included in the Command manpower measure is the direct 
effort associated with technical manpower, derived not from 
the System for Resource Management (SRM), but as a direct 
translation of time expended by personnel in the various tasks 
associated with software implementation such as designing, 
coding, testing, and documentation. The time spent by engi- 
neering personnel on documentation is included but not the 
costs of support people, reproduction, and other costs 
involved in the production of the documents themselves. 
Time expended was derived primarily from the DCD Opera- 
tional Software DeveLpment Schedule and Manpower Alloca- 
tion Chart. 

Figure 4 represents the manpower loading of tiie DCD 
during the time span under discussion. Studies of available 
data from industry (Refs. 5 and 6) indicate that an average of 
3.2 persons per month (200 LOC/mm) should have been 
utilized to produce the necessary code. Instead only an average 
of 2.7 (237 LOC/mm) persons were employed. Over the 
25-month period, a cost savings of 1 person was achieved, or 
approximately $50,000 for one subsystem. 

A second subsystem was studied to corroborate the pro- 
ductivity gains found in the Command Subsystem. The Com- 
munications Monitor and Formatter (CMF) was chosen and 
the CMF man-power loading is shown in Fig. 5. Again, based 
on an industry productivity average of 200 lines of code per 
month (Refs. 5 and 6), the savings were 17.3 man months. 

The increase in productivity can be considered significant 
for both subsystems because it exceeds one standard devia- 
tion from the r^an. The result reinforces the belief that 


small but signitlcant direct implementation cost savings were 
obtained with the use of the software methodology. 

IV. Software Development Period 

The development period of 25 months used lo develop 
cost estimates is an extremely conservative one consisting of 
starting where the first known manpower could be attributed 
to the project. There are many valid and cogent reasons why 
a later date may be taken. The implementation schedule 
(Fig. 6) published at the High Level Design Review shows a 
start date for software implementation at approximately 
January 1975. if this date is taken as the start of the develop- 
ment period then the following productivity figures can be 
derived: 

DCD 256LOC/month 

CMF 274 LOC/month 

The savings would then be: 

DCD 1 7.4 man/months 

CMF 23.3 man/months 

or approximately $167, OC ) for just the two subsystems, 

V. Documentation 

As stated previously, the productivity measure includes 
engineering time spent on design, code, testing, and docu- 
mentation, but does not include special documentation 
support required by the software methodology. For the MDS, 
documentation support was provided by the Software Produc- 
tion Management and Control (SPMC) center, and is dificussed 
below. 

Despite the methodology requirements to produce docu- 
memation before or at least concurrent with code, the delivery 
of formal, validated documentation has consistently lagged 
behind transfer of the code. This lag has, for many subsysi ^ms, 
been as long as a year and is partially due to the volume v f 
documentation produced 

Walston and Felix of IBM (Ref 7) predict on the basis of 
46 software projects of varied purposes and methodologies, 
that about 800 pages of docun^entation is produced for a 
subsystem the size of Command. The values for one standard 
deviation are 400 pages and 12i0 pages of documentation. 
The Command Subsystem produced at least loOO >'ages 
(Table 2). This value is higher than the average found by IBM, 
and indicates that the methodology results in a higher volume 
of documentation. 
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Extrapolation of information supplied by DSN Data 
Systems indicates that documentation support^ for the DCD 
cost $69,000. This is slightly more than the productivity gains 
shown in Section III, and does not include costs incurred after 
tire time*span of this study. Although the documentation cost 
is negating productivity gains in the implementation phase, 
these costs should be offset in the maintenance phase. In the 
absence of hard data, the cost and volume of documentation 
required by the methodolog> is being analyzed. 

VI. Maintenance Phase 

Unfortunately, the maintenance and operations ohase of 
the MDS system has just begun, so there is no quantitative 
cost data. Nevertheless, two pieces of information indicate 
that maintenance of the software will be less costly per LOG 
than previous projects. First, though given extensive testing, 
the number of anomalies was low. Secondly, another project 
produced concurrently with the software methodology 
showed increased software maintainability. There is no reason 
to expect a decrease in maintainability, and every reason to 
expect an increase with the use of the software methodology. 
Thus, it is predicted that there will be additional savings during 
the maintenance phase of the MDS. 

VII. Conclusions from Quantitative Study 

The major quantitatively measurable benefits derived from 
the methodology discussed in this report appear to be three- 
fold: 

(1) The ability to deliver software on time by providing 
greater visibility and allowing management to inter- 
vene before projected slips actually impact delivery 
schedule. 

(2) A lower anomaly ratt than the reported industry 
average, which indicates b’^tter software and predicts 
decreased maintenance costs. 

(3) Small but significant cost savings in manpower result- 
ing in lower implementation costs. These savings will 
increase as the documentation problem is addressed. 


^ Docuiiientalion support, as used here, is defined all die support 
functions performed by the SPMC for th^ Command Subsystem. 
Thus, included in the $69,000 figure are the following: documenta- 
tion production, editing, and reproduction costs; com^^uter operator 
costs for disk generation, job runs, source code updating, diK prepara- 
tion, operaung system copies, and module assembly; and code support 
costs such as the maintenance of the code library. 


Vli;. Qualitative Evaluation— Benefits 

In addition to the quantitative evaluation of the software 
methodology, a qualUative study was performed using inter- 
views and questionnaires of cognizant software individuals. 
Qualitative, rather than quantitative, evaluations were sought 
for several reasons. Record-keeping is a part of the methodol- 
ogy, but quantitative data are often unavailable for older pro- 
jects. It is also difficult to obtain numeric information from 
people recalling old project? : often those who directed old 
projects are no longer at JPL, an^ daca froit? memory may not 
be sufficiently accurate for a valid 'comparison. Thus, the 
best information to be gained from these individuals is 
qualitative. 

A specific questionnaire and open discussion format were 
used in the interview process. The programmers and managers 
were also asked to relate their experiences on projects con- 
ducted prior to addition of the software methodology, and 
specify area3 in which the methodology would have been 
beneficial. The results of the interviews show a number of 
common impressions and experiences, which are discussed 
below. 

The primary problem observed in older software projects 
was a failure to transfer to opeiations on schedule. Often 
accompanying this overrun in schedule was an overrun in 
budget. With these two conditions, and pressure to transfer as 
soon as possible, documentation suffered and maintenance was 
difficult once transfer was achieved. The people interviewed 
noted a number of causes that produced this condition. 

The first was the general difficulty in accurately predicting 
task complexity, required resources, and the final amount of 
code. Cost and schedule overruns occurred and some incom- 
patibilities with hardware resulted. This problem not only 
affected past DSN projects, but still affects the software 
industry as a whole. Nevertheless, it is felt that the meth- 
odology's requirements for modularization and Work Break- 
down Structure, with a strong emphasis on design, strongly 
support better schedule, cost, and resource planning. 

A second cause of overruns noted by the software per- 
sonnel was frequent changes in requirements. This instability 
disrupted schedules and necessitated numerous releases and 
transfer liens. The software methodology addiesses this 
problem by fonnalizing the requirements process and enforc- 
ing approval of requirement changes. 

Occasionally noted was the condition in which the project 
was due for transfer and management then learned it was far 
from complete. Management had not adequately monitored 



the progress of the project. Again, a benefit of the methodol- 
ogy is that it includes the Work Breakdown Structure, design 
reviews, and formal milestones, which all address the monitor- 
ing of project progress. 

The software projects discussed in the interviews all had 
late and/or insufficient documentation, and major transfer 
liens were frequent against these documents. Numerous prob- 
lems resulted: the project was delayed, new personnel found it 
difficult to work on the software, or leiting and debugging 
required excessive resources. It is felt that the concurrent 
documentation and deeign review reqi .rements of the meth- 
odology certainly ease these problems. Nevertheless, some 
I documentation is still late. This is partly attributed lo flow- 

charts that sometimes require as long as a month for machine 
production. By that lime, they must often be resubmitte^ to 
include changes that occurred during that n*«»nth. A faster 
method of documentation p’-oduction was stiongly suggested 
by the software personnel. 

The final cause of p**" .^t ')verruns brought out in the 
interviews was excessive errors in the software. Either major 
pretransfer testing was required, or numerous bugs were found 
after transfer. Neither case is desirable. The personnel inter- 
viewed felt the modularization, emphasis on design, and 
/ formal testing incorporated in the methodology result in less 

testing, fewer anomolies, and easier debugging. 

In summary, projects in the past hav.; often encountered 
“I diffic ulties with schedules and resources, specification changes, 

management visibility, documentation, and testing. The DSN 
Software Methodology provides a good structure with which 
mtrol these problems. 

Qufiiitative Evaluation- 
Suggested Improvements 

In addition to benefits resulting from the software meth- 
odology, the software people suggested the following meth- 


odology improvements that would produce even more cost- 
effective software: 

(1) The methodology requirements are occasionally not 
applicable to a specific project, and obtaining a waiver 
requires time. Perhaps, with a growing experience in 
using the methodology, a list of requirements by 
project size, cost, and type could be developed. 

(2) The projects still find themselves caught between 
documentation requirements and the scheduled trans- 
fer date. Part of the problem is felt to lie in the slow 
turnaround time for finished documentation. Also, 
the requirement of low-level flowcharting results in 
a great deal of production and revision time. It is 
suggested that either (a) only higher-level flowcharting 
be required, (b) a program design language tool be used 
in place of flowcharts, or (c) a quick method of flow- 
chart production be found. In all, the software per- 
sonnel feel that flowcharting currently requires that 
the programmer cater to the machine, rather than the 
reverse. 

(3) The inclusion of programmers, operators, C yiizant 
Development Engineers (CDEs), and customers in the 
requirements and broad design process should be 
encouraged even more than is currently done. 

(4) Too much emphasis is given to the quantity of 
response to the methodology rather than the quality 
of the response. An effort should be made to deter- 
mine levels of quality based on maintainability and 
error rate after transfer. 

(5) The methodology can help standardize the develop- 
ment of software and control the flow of projects, 
but it cannot force a bad manager or programmer into 
good performance. A key also lies in the people who 
perform the job. 
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Sybsy$tem 


Lines of code 


DCD 

Command 

16,000 

DTM 

Telemetry 

30,100 

DTK 

Tracking 

12,000 

DMC 

Monitor 

16,000 

DTT 

Test and Training 

2,100 

CMF 

Conununications Monitor and Formatter 

17,700 

DST 

Host 

5,000 

Total lines of code* 

98,900 


*The Standard Operating System is not included in the total lines of 
code. Size of the Sundard Operating System is 16,500 lines of code. 




Document 

Number of pages 

Project Notebook 

154 

Software Requirements Document 

43 

Software Definition Document 

123 

Software Specifications Document^* 

780 

Software Operations Manual* 

140 

Software Test and Transfer Document* 

250 

Anomaly Reports 

127 

Total Documentation 

16l7 Pages 

* A mount of documentation was estimated because there wa^ no for- 
mal, published document during the time-span under study. 
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MBASIC® Batch Processor Architectural Overview 

S. M. Reynolds 
DSN Data SyttwTw Section 


The MBASIC^** Batch Processor will allow users to run MBASlC^^ programs more 
cheaply. It will be provided as a CONVERT TO BATCH command, usable from du ready 
mode. It will translate the users inogram in stages through several levels of intermediate 
language and optinuzation. Only the fmal stage is machine deperuient: therefore, only a 
mail effort wUl be needed to provide a Batch Processor fora new machine 


L Intfoduction 

Many large, frequently used programs, which are slow and 
expensive to run, currently exist in the MBASIC*^** user 
community. The MBASiC'^** Batdi Processor will reduce the 
time and cost associated with running these programs by 
providing a compiling facility witliin the MBASIC'*'^ environ* 
ment which will allow a user to convert programs to a directly 
executable form. Direct execution typically reduces program 
tun times by at least an order of magnitude below those of the 
same programs executed interpretively. 

Within dw MBASIC^** enviromnent, the user may invoke 
die Batch Processor through the following conuiumds, pre- 
scribed in fUndamenttds ofMBASlC^** , Vd. fl (JPL, Feb. !9, 
1974): 

CONVERT OBename) TO BATCH (fflename) 

CONVERT TO BATCH (fflename) 

(JUEUE BATCH (fflename) 

RUN (fflename) 


II. MBASIC® Batch Processor Internal Vtaw 

When called, the MBASIC^** Batch Comdler receives the 
pseudo-op code (P(K*CODE) buffer, the line Inference Table 
(LRT), and the Name Reference Table (NRT) from the parser. 
The source code of the program to be compiled is never exam- 
ined. The output of the compiler is directly executable code 
for the host machine. 

Within the compiler, a completely machineindependent 
stage is foQowed by a machine-dependent stage. A major goal 
of the compiler design is to minimize the proportion of die 
machine-dependent sti^. 

The machine-independent stage uses the PC^OI^, LRT, 
and NRT to generate a low4evel, machine4ndependent assem 
bly language called MCODE, which is then tninsiated by the 
miKhine-dependent stage into host code (see Fig. 1). 

The first step in the translation of POPCODE to K«X)raE is 
traniiation of POPCOK to an eqdvatent thiee-addiesi code. 
This intermediate notation it easier to handle in the ensuing 
optimization step than the reverie ptfflsh notation of P(V- 
CODE because refeienc' operands are exididt and are 


diwctly attodated with the operation code rather than 
implicit in pushes onto the stack by previous op codes as in 
POKODE. 

The second step in translation to MCODE is the optimiza* 
tion of the intermediate code. Because the input and output of 
this step are intermediate code, expressed in the same internal 
representation, this step will be delayed in the initial imple- 
mentation and added later when the sunounding steps hive 
been verified. 

The third machine^dependent step is translation of the 
intermediate code to MCODE. 

Translation of MCODE to machine-dependent assembly 
language is performed in two steps: direct translation of 
i MCODE to host code followed by machine-dependent optimi- 

: zation. 

f For example, the source statement 

1 A=B+C 

i 

! becomes 

i 

I PUSHB 

I PUSHC 

t ADD 

I POPA 

\ 

I in POPCODE. This is translated into 
I ADD A3,C 

I in the intermediate code, vtdiich is translated into 

FLOAT TEMPI ,B 

ADDA,TEMP1,C 

in MCODE (assuming that B is declared to be an integer), and 
into 

LMJ X10,FLOAT 

B 

TEMPI 

LA0,TEMP1 

AA0,C 

SAO,A 

Befwe translation begins, a declaration processing pass 
through the P(VCODE builds a Symbol Definition Table 
(SDT), which contaiiM type and buu^g informatimi for each 
var‘*l^ explicitly or implicitly declared in the program. The 
translation process is performed one statement at a time in 


order to minimize the size of buffers and the use of external 
mass storage. 

HI. PtMKlfrOpCode 

MBASIC^^ pseudo-op code it a reverse polish notation for 
the MBASIC'*^'* virtual stack machine (VSM). One POPCODE 
string it produced for each MBASIC^^ source statement. The 
MBASIC^^ VSM checks for incompatible types at execution 
time and performs legal type conversions; it can handle arrays 
and strings as atomic data items; and it handles FOR-NEXT 
and GOSUB transfer; implicitly by using auxiliary stacks. The 
reflection of these features in POPCODE makes it a relatively 
hifh'levd language, requiring con^derable translation to reads 
the ievd of a typical red machine. Even syutactically correct 
MBASIC'^^ programs may translate to POPCODE programs 
vdiich generate runtime errors in the VSM. These errors wQl be 
detected and an enor message printed during trandation rather 
than at runtime of the compiled program, whenever possibh. 

IV. kitmuMdiate Languag* 

The Intermediate Language (IL) is a direct mapping of 
POPCODE to an equivalent non-stack-dependent language. All 
implicit references to stack data are replaced by explicit 
references to identifiers, either source defined or temporary 
identifiers defined by the IL processor. 

With one exception, the translation from POPCODE to IL 
involves a one-to-one mapping of POPCODE statements to 
corresponding IL statements. The exception is that POPCODE 
operations vhose only result is to push a variable's value onto 
the data stack do not appear in the IL string, but the 
identifiers are referenced as IL arguments. 

V. POPCODE to IL Translation 

The POPCODE to IL translaticm algorithm is a pseudo- 
execution of the POPCODE string using an identifier stack 
instead of a data stack. The popcode string is scanned, and 
each POPCODE encountered generates the corre^Kmding IL 
If values would have been popped from the data stack during 
executim, identifien are popped from the identifier stack in 
the same order. If values would have been pushed onto the 
stack during execution, temporaries are created and their 
identifiers pushed onto the identifier stack. 

VI. MCODE 

MCODE is a low-level, three-address code for a daai of 
virtual machines with variable storage attributes, such as: 

(1) Differing number and type of registers. 
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(2) DUTering mapping* from main storage to registers (e.g., 

byte addressable machines which map several cells into 

one register and word addressable machines vdiich map 

one cell into one register). 

(3) Differing main storage size. 

The MCODE generator for a specific im(dementation of the 
Batch Processor will be configured by setting its variable 
attributes to match the host system as closely as possible. 

The MCODE virtual machine recognizes the following data 
types: 

(1) Integer - at l<^ast 32 bits precision 

(2) Real - host system floating pmnt format 

(3) Pointer — at least 16 bits precision 

(4) Char — at least 7 bits precision 

MCODE machines have a separate address space for each 
data type and for instruction code, and the unit quantity for 
painter type and for any absolute address represents one ceU 
of address space of the type p<dnted to. 

Registers may be defined to accept one or more different 
data types, e.g., a general purpose register in many machines 
would be represented in the MCODE machine by a register 
accepting real, integer, and pointer data. 

MCODE contains binary, unary, and nonary operations of 
the following kinds: 

(1) Arithmetic 

(a) Add 

(b) Subtract 

(c) Multiply 

(d) Divide 

(e) Unary minus 

(2) Bodean 

(a) And 

(b) Or 

(c) Exclusive or 

(d) Complement 

(e) Conditional jump 

(f) Unconditional jump 


Each instruction is defined separately for each type, and all 
operands mutt be of matching type. Mixed type operations are 
possible through use of the fdlowing conversion operations: 

(1) Integer to real 

(2) Real to integer 

(3) Integer to char 

(4) Char to integer 

(5) Pointer to integer 

(6) Integer to pointer 

Each instruction consists of an operation code and zero to 
three operands. Each operand contains a register/address flag, 
a register number or address, a type, and a boolean indirect 
addressing indicator, v^ch may be true only if the operand 
specifies a pointer register. For arithmetic and boolean 
operations, the P's! operand specifies a destination for the 
result of the operation, and the two other operands specify the 
source locations. For unary operations, only one source 
operand is necessary. Jump instructions have one operand, 
v^ch may only be a pointer register in the case of an indirect 
address reference or an address in the case of a direct address 
reference. 

MCODE also contains the following operations on data 
files: 

(1) Open 

(2) aose 

(3) Remove 

(4) Rename 

(5) Append 

These operations have only one operand, \^ch contains a 
file name. 

The following operations have a file name operand followed 
by a source or destination operand: 

(1) Get char 

(2) Put char 

VII. IL to MCODE Tramilation 

IL to MCODE translation is an in-line expansion performed 
in a single scan of the IL sequence. Each IL is expanded into 
one or more MCODE statements. Code for type checking and 
type ctmversions for legal mixed type operations is generated 
during this process. Error messages ‘or illegal types are 



generated during type checking. The IL arguments generated 
are selected on the basis of the best known itontce for c value 
at compiled code execution time, i.e., a value known to be in 
both register and main storage will be represented by the 
register. 

VIH. MCODE to Host Assembly Language 
TnmsMion 

bi most cases, the host code will be its standard assembly 
language. MCODE to host code translation is a machine- 
dependent process; however, in most cases it can be performed 
in one scan of the MCODE string with one or more host code 
statements generated for each MCODE statement. This is 
possible because MCODE is generated for a virtual machine 
having storage characteristics similar to those of the real 
machine. 

IX- Machine-Independent Optimization 

The machine-independent optimization step is transparent 
(except in efficiency) to the rest of the MIDB processor 
because its input and output have the same representation. 
Optimization will concentrate on expressions and assignment 
statements. Many statistical studies of actual source code and 
the possible optimizations have been made, and all show these 
to be the most effective optimizations in terms of actual 
reduction of the size and execution time of typical object 
code. The selected optimizations are: 

(1) Subexpressions appearing more than once are calcu- 
lated once and saved in a temporary variable; instances 
of the common subexpression are replaced by a 
reference to the temporary variable. 

(2) Subexpressions involving only constants are calculated 
at compile time. 

(3) Multiple assignments to the same variable are recog- 
nized and replaced by the last assignment only. 


X. Machine-Dependent Optimization 

Machine-dependent optimization is also a transparent pro- 
cess (except for efficiency). Some possibilities common to 
most machtoes, but still machine-dependent, are the following: 

(1) Gudns of branches are replaced by a single branch 
instruction. 

(2) Unreachable code is deleted. 

(3) Common code groups preceding a collecting node are 
recognized, and the nod^^ is moved to precede the first 
instance of the common group. The other instance is 
deleted. 

(4) Common code groups heading all branches of a 
program fork are condensed as in (3), above. 

Other completely machine-dependent optimizations might 
involve such things as invoking auto increment mode in place 
of explicit pointer arithmetic or recognizing LOAD-STORE 
groups and replacing them with direct memory to memory 
moves. 


XI. Summary 

The MBASIC^^ Batch Processor is a language translator 
designed to operate in the MBASIC'*'^ envirorunent. It will 
provide a facility with «4uch MBASIC^^ programs which 
have been debugged in the interpretive mode may be trans- 
lated into a directly executable form and stored or executed 
from the MB ASIC’*'” environment. 

The Processor is to be designed and implemented in both 
machine-independent and machine-dependent sections. The 
architecture is planned so that optimization processes are 
transparent to the rest of the system (except for efficiency) 
and need not be included in the Hrst design-implementation 
cycle. 
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GCF HSD Error Control 

C. -K. Mung 
OSN Otu Syttsms Section 


A selective repeat Automatic Repeat Request (ARQ) system has been implemented 
uruier software control in the GCf Error Detection and Conation (EDC) assembly at 
JFL and the Comm Monitor and Formation (CMF) a^embly at the DSSs. The CMFand 
EDC signiflcantfy improved real-time data quality and sigrdficantly reduced the post-pass 
time required for replay of blocks original received in error. Since the Remote Mission 
Operation Centers (RMOO) do not provide compatible error correction equipment, error 
correction will not be used on the RMOC-JPL HSD circuits. The real-time error correction 
capability will correct error bunt or outage of two f2) loop-times or less for each DSS- 
JPL HSD circuit. 


I. bitroduction 

The Ground ConununicatJons Facility (GCF) High Speed 
Data (HSD) Subiyatem coniUtt of the GCF Anembliea 
(Ref. S) used to switch, transmit, record, process, distribute, 
test, and monitor d^ital data. The detailed functional capa- 
bilities of the GCF HSD subqrstem for the Mark III period are 
illustrated in Fig. 1. The HSD drcuits for interchanie of data 
between the DSS/RMOC and JPL are impressed widt serial 
Ut streams at a continuous line bit rate of 7200 bits/sec in a 
1200-bit block size. The blocks are either data blocks or filler 
Mocks. Filler blocks ate inserted when the user data load is 
insufficient to maintain contiguous data blocks on line. 

A selective repeat Automatic Repeat Rerpiest (ARQ) 
qrstem ^ef. 6, 9, 10, 11) in «Mch on^ the particular Mock 
received in error is retransmitted has been impiemented under 
software control in th* GCF Error Detection and Correction 
(EDC) minkomputers at tiie Central Communicatloos Termi- 


nal (CCT) at JPL and the Comm Mcmitor and Formatter 
(CMF) minicomputen at the DSSs. (These minicomputers, 
a major portion of GCF High Speed Data sub^stem, are 
shown in Fig. 2.) The CMF and EIX minicomputers, working 
together, significantly improved real-time data quality and 
significantly reduced the post-pass time required for re^ay 
of blocks originally received in error. Ihe real-time error 
correction system cmrects error burst or outage of two (2) 
loop-times* mless. 

Error correction requites ths* both ends of the HSD line be 
equipped with enor correction capabiUtiet. Error correction 
win not be used on the RMOC-GCP CCT HSD circuits as the 
RMOCs do not provide compatiMe error correction 
equ^ment. 


'a loop-tac is a time between a Mock tiaaimittad and the block 
aeknowMsement leoehed. 
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N. QCF Error Control Algorithm 

In describing the GCF Error Control Algorithm (Ref. 1 , 2 , 
7, 8), only one direction of transmission/rcception will 
be discussed. Ihe reverse direction is identical though inde- 
pendent. A retransmission schematic diagram (Fig. 3) is given 
to aid in understanding the error control algorithm. The 
following buffers are needed for the retransmission scheme 
implemented: 

(1) Acknowledgement Queue (A-Queue) 

This queue is a flrst-in/first-out queue which contains 
data biocks that have been transmitted over the HSD 
line and are awaiting acknowledgement. As soon as 
an acknowledgement indicates that a block has been 
received error-free at the other end, that block is 
dropped from the A-Queue. Thus, at any time, no 
more than a loop-time of blocks are stored in the 
A-Queue (see Table 1, DSS-JPL Loop-Time biocks). 
The top-entry of the A-Queue is called the A-candi- 
date. It is the next block to be acknowledged or to be 
retransmitted if the block is not positively acknowl- 
edged. For each acknowledgement received, the 
A-Queue pushes up thus creating the next A-candidate. 

(2) Transmit Queue (T-Queue) 

This is the data buffer which holds the incoming data 
blocks until they are called for transmission over the 
HSD line. 

(3) Delivery Queue (D-Queue) 

This queur stores all error-free HSD blocks received 
during ret.ansmission until all prior blocks have been 
successfully retransmitted. The data blocks are 
delivered in order. When error correction is inhibited 
in the far-end, the data blocks are delivered in first-in/ 
first-out sequence. 

The error detection encoding and '.ecoding under software 
control may use either a 22-bit or 33-Hit error polynomial. 
Error correction works only with the 22-bit error code. Each 
block has two 8-bit fields reserved for GCF error correction 
use. One field (bits 9-16, word 73) contains the GCF serial 
number of the block. The other field (bits 1-8, word 74) 
contains the serial number of the just-received block which is 
being acknowledged as error-free. These two numbers are 
termed the GCF serial number and the GCF acknowledgement 
number respectively. Since the GCF serial number of a HSD 
block has 8 bits, GCF serial numbers are selected from the 

set: 0,1,2, ,255. lire term “in GCF sequence” means 

that a block is serialized in sequence with the previous trans- 
mission using the set: 1,2, ^55,1,2, ,255,----. 

The zero value is skipped since it is used for filler blocks. 


Data blocks transmitted in the error correction mode will 
carry the non-zero GCF serial numbers. Dat<( blocks (rans- 
mitted in the error correction inhibited mode will have no 
GCF serial number (i.e., zero). Whether or not a data block 
carries a GCF serial number depends on whether or not the 
transmit end is in error correction mode. Whether a data or 
filler block carries a GCF j;cknowledgement number depends 
on whether a data blc with a non-zero GCF serial number 
was just received. Filler blocks normally will not be serial- 
ized but may carry acknowledgement numbers. Each data 
block that is not acknowledged (i.e., is in error) will be cor- 
rected up to N times (normally, N-2). No delay in data flow 
is introduced when the data string is error-free, but delay is 
introduced when error correction by retransmission occurs. 

The GCF error control algorithm is incorporated into the 
receive and transmit functions of the CMF and EDC programs. 
The following is rite logic of the error control a’ tlun. 


A. ftoeciv* Stops 

(1) Accept c high speed data or tiller block from HSD line. 
A normal sequence of interrupts occurs within the 
time of receipt of a block as predicted by the line rate 
(6 blocks/sec). 

(2) Decode. A block may be received in error in two ways: 

(a) Decoding error - Transfer of a complete block, 
but the 22-bit enor polynomial does not corre- 
late with the data in the block. 

(b) Missing block - A normal sequence of interrupts 
does not occur within the time of receipt of a 
block as predicted by the line rate. The time slot 
(166.66 msec in 7.2K bits/sec line rate, 6 blocks/ 
sec) conesponding to this missing block is treated 
the same as a block in enor. 

(3) Save GCF senai number and acknowledgement 
number. 

(4) If near-end is in enor correcoon .node, go to 5; other- 
wise skip to 6. 

(5) If the acknowledgement nvunbet matches the A-candi- 
date (in A-Queue) serial number ^indicating that the 
block was received error-free at the distant end), then 
the A-candidate is deleted from the A-Queue and 
retransmission is not required. If the numbers do not 
match, retransmission is needed. 
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(6) If the far-end is in error correction mode and ihe block 
is received in error, then discard this block and go to 9; 
otherwise go to 1? 

(7) Test block type. 

(a) Discard filler block. 

(b) Data block with GCF serial number (non>zero) — 
consider far-end in error correction mode and 
deliver block to D-Queue foi ordered delivery to 
user. 

(c) Data blocks without GCF serial numbers (zero) - 
consider error correction inhibited at far-end, 
deliver blocks to D-Queue for first-in/iirst-out 
delivery to user. 

(8) Deliveiy to the user (I>Queuc Management). The 
D-Queue is managed by a timer. This queue holds all 
received data blocks. After processing in the enor 
conection mode (far-end), the blocks are delivered in 
GCF serial number sequence (the same order in sidiich 
the blocks Wi e originally transmitted). With error 
correction inhibited, the data is delivered in a first-in/ 
first-out sequence. 

(9) End of receive task. 

B. Transmit Steps 

(1) If near-end is in error conection mode, go to 2; other- 
wise skip to 3. 

(2) If retransmission is required (from receive task, step 5), 
the A-candidate is processed in the following way: 

(a) If the A-candidate is a daU block (GCF serial 
number is non-zero) and if this A-candidate has 
been retransmitted N times, then delete the 
A-candidate and go to 3; otherwise, retransmit the 
A-candidate and go to S. 

(b) if the A-candidate is a filler block (GCF serial 
number is zero), delete the A-candidate and go 
to 3. 

(3) If the T-Queue is net empty, select the top-entry from 
T-Queue for transmission; otherwise select a filler 
block for transmission. 

(4) If the near-end is in the error correction mode, insert a 
GCF serial number if a data block is being transmitted; 
otherwise, insert zeros in the serial number slots. 


^In case of 7, either (1) fareud is error correction inhibit mode, or 
(2) is in enui correction mode and the block is received 

eno^free. 


(5) Insert the serial number from the last received block 
(from receive task, step 3) as the GCF acknowledge- 
ment number. 

(6) Encode in 22-bit rnor code. 

(7y Transmit this block to HSD line and store this block 
in the A-Queue if near-eiid is in error correction mode. 

(8) End of transmit task. 

In the GCF error control algorithm, the transmit task is 
driven by the receive task. That is, each received block induces 
an output block containing its acknovdedgement. 

The GCF error control algorithm results in four (4) possible 
operaticHial modes for each DSS-JPL HSD circuit: 

(1) No error correction m either direction (DSS/RMOC 
^^KJCFCCT). 

(2) Inbound error correction only (DSS‘^GCF CCT). 

(3) Outbound en r correction only (DSSHJCF CCT). 

(4) Error correction both directions (DSS'*KJCF CCT). 


Hi. GCF Eiror Cemtroi hnplemef^ 

The physical configuration and majoi components of the 
GCF HSD Subsystem in the Mark 111 Era is illustrated in Fig 2. 
The GCF error control algorithm has been implemented in 
the CMF and EDC software. Key aspects of the implementa- 
tion are: 

(1) The software has been designed and implemented in 
accordance with DSN software methodology standards 
(top-down implementation, structured programming 
principles, etc.). 

(2) The implementation effort was organized into a work 
breakdown structure, which was then used to monitor 
progress. 

(3) Initial intercomputer protocols (data flow) and timing 
analysis were completed very early in the software 
development. 

(4) Software testing was accomplished in a hardware 
environment which closely matched that of the final 
operational configuration. 
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After the conclusion of acceptance tests of the DSS CMF in 
September 1977 and GCF CCT EDC and HSW software In 
November 1977, the GCF HSD subsystem for the Mark HI 
period became operable and useable to support DSN missions. 
As of July 1978, all Voyager testing, ATRS recalls, and 
Manual DODR replay were used by the GCF HSD error correc- 


tion capabilities. On September 15, 1978, Voyager com- 
menced use of error correction on a routine basis. 


The performance of the GCF HSD error control will be 
described in a subsequent issue of the DSN Progress Report. 
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Tracking Error of 100-m Antenna due to Wind Gust 

M. Massoudi 

OSN Enginaering Stction 

A procedure is shown to derive root-mean-square tracking error for an antenna system 
in response to wind-gpst loading. Example calculations are ittustrated for the design of a 
proposed 100-m antenna. The effect of vnnd-gust correlation is also considered. 


I. Introduction 

Wnd loading is an important factor in the design of many 
above-ground structures. For large antenna structures in 
particular it is necessary to make some allowance for possible 
wind loading. Wind motion may be described as consisting of 
two parts: a mean motion and a superimposed turbulent 
fluctuation. Some studies have provided detailed measure- 
ments of the time history of wind fluctuations, and techniques 
of random process theory have been employed to gafii insist 
into the physical processes controlling the turbulent fluctua- 
tions. Hiis article defines and calculates the tracking error of a 
proposed 100-m-diameter antenna due to wind-gust torque as 
a fimction of equivalent spring constant. 


S^{f) = output spectral density 
^rin» ” tracking error due to wind torque 
/ = the imaginary unit 

The root-mean square (rms) value of the output is obtained 
from 



where S^{f) can be computed by 


II. Afitmina Tricking Error 

Wind*guit torques are analyzed as a random process with 
zero mean. For a system with torque input and angular 
displacement output (m Fig. 1), we define 

- mean torque input due to wind, Nm 
• input spectral density for dynamic wind load 

y(Jf) - aystem admittance to wind torque 


Sjif) - (2) 


in which 



^ -.( 3 ) 





« series combinatiQn of structural stiffness and servo 
drive stiffness 

» servo position loop bandwidth, Hz 
natural frequency of the structure, Hz 
i » damping ratio, relative to critical damping 
/« excitation frequency, Hz 

To define the sptctnl density of the dynamic load, the 
spectrum of the velocity fluctuations is first requited. The 
velocity q>ectral density used in this report was obtained by 
Davenport (Ref. 1), giving the spectrum in the form 


4KUfX^ 


(4) 


in v^ch 


^y(f) ^ velocity spectral density (one-sided) 

= steady speed of wind at reference hei^t of 10 
meters, m/$ 

K - surface drag coefficient for the local terrain 
X * 1200//J/, 


in which 



_ 3.8S /A2T 11.5 /AT 3.8S/AZ 

U U u 


( 8 ) 


and 




(l+«) 



(9) 


in which /iX, AY, and A? - the alongwind, crosswind, and 
vertical distances over which the correlation is being taken, h is 
the hei^t of structure, and a is the exponent in the assumed 
power low-speed profile of the wind vs hei^t of the structure 
relationship given by 


U, 



( 10 ) 


To relate the load spectrum to the velocity spectrum, the 
lelatimi given by Vetlozzi and Cohen (Ref. 2) is used, giving 
the gust-load spectrum in the form: 


A. Calculations 

To develop Eq. (4) specifically for the case at hand, the 
following parameters are used: 


4r* 


a»0.140S ) 

(5) } for Goldstone site (Ref. 3) 

A: = 0.005 I 

Z=48m 


in udiich T~ and are the steady load and mean wind speed 
at height Z (meteit. 

Since the wind gusts are not likely to act simultaneousiy 
over the full height of structure, a correlation function 
(Ref. 2), C^(Ot will be introduced in Eq. (S), whidi is 
dierefoie rewritten as 


dien from Eq. (10), 

/lO\o.Mos 

or£/,=0.80K^ 

with these relationships, the specific representation of Eq. (4) 
becomes 








For examine, with = 4.47 m/s. 


y, =0.80 (4.47) *3,58 m/s 


X * 1200^—= 335,57/ 
3.58 


Table 2 shows the results of Eq. (11) aj^lied to several 
mean speeds using the reference speed and torque from the 
wind tuimel analysis. 


D. Calculation of K«q 

The equivalent spring constant is computed by means of 
the wetl-lmown relationship for springs in series 


and the gust velocity spectrum is given by 


28800/ 


(1 + 112608/*) 

Figure 2 shows wind vdodty spectral density curves for a 
famify of values of mean qiieeds, V^. Table 1 gives the spectral 
density equations derived to {dot the curves. The tabulated 
values of /„„ and 5|,(/g,„) indicate the frequency at which 
the spectra peak and the assodated magnitude, respectively. 


Structure 


From structural and servo system analysis, the following q;»ing 
constants were derived: 




*6.78X10‘®-5^^ 

servo 


B. DorIvoMon of the Wind Torque 

Loading data derived frmn existing antenna wind tunnel 
tests was a{>{died to an ana^cal model of the 100 m antenna. 
The anterma orientation drown in Fig. 3 results in the 
maximum wind torque^pn the antenna dish, at a speed of 
31.29 m/s (70 mph) the resrdt was 

= S34J X 10* Nm at 31.29 m/s 

(3.951X 10*lb-ftat70mph) 

C. CatculaMon of Wind Torque for Different Spood c 

To calculate wind torque for different s{>eeds of wind, the 
fdlowing formula expresses the relationsUp that torque is 
proportional to the square of wind speed, e.g.. 


(fc)' 


vrtiere rr.f and are reference torque and speed, 
respectively. 


substituting these values in Eq. (12) provides 


- 1.99 X 10*®-i^^ 
*1 rad 


E. Calculation of the TracMng Error 

The digital computer was used to integrate Eq. (1) numeri- 
cally and calculate Example results are listed as follows: 

Case 1: 4 * 0.20 H2,4 = 1.00 Hz, f = 0.01 

C*(/) * 1 ( 100 % corr^ation) 

Table 3 shows the constants that can be used to compute the 
tracking error for Case 1. The tracking error is obtained by 
dividing the constant by K^„. The results of 1 are j^otted 
in Fig. 4. 

Case 2: 4 » 0.20 Hz,/, - 1.00 Hz, f » 0.01 


^*(/): ‘ 
I 


h « 104 meters, AAT ■ 50 meters 
AY * 50 meters, AZ * 104 meters 


Noies h and AZ are the overall height of the antenna at stow 
|K>sition. 
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Table 4 ihowa the tracking error conatanta for Case 2. The 
reaulta of Caae 2 are plotted in Fig. 5. 

W. Conduslon 

Hie tracking error of the propoeed 100^*dianieter antenna 
due to wind guat has been defined and calculated. In Caae 1, in 
which wind guat haa 100 percent correlation over the antenna 
structure, the tracking error is about 0.012 deg. at 13.14 m/s 


(30 mph) compared to 0.00043 deg. of Case 2, which is a 
more realistic result in terms of wind-gust effiects on the 
antenna. 

Of course, in considering the above results, one should bear 
in mind that a reasonable toleranoe margin must be taken into 
account, and these results merdy give some understanding to 
the design team about the anteniu performance with respect 
to wind gust. 
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Tsbte 1. rimlty of wind vilocity ipsctral dniltlM 



sv^n. 

■^max* 

^K(/max)» 

m/s 

(m/s)^/Hz 

Hz 

(m/$)2/Hz 

4.47 

28800/ 

(1 + 112608/2)^/’ 

2.31 X 10"3 

35.52 

8.94 

28800/ 

4.62 X 10"5 

71.05 

(1 + 28152/*)^/^ 



13.41 

28800/ 

6.92 X 10"^ 

106.57 

(1 + 12512f^)*l^ 

17:88 

28800/ 

9.23 X 10“^ 

14110 

(1 + 

.7135 

28800/ 

1.15 X 10‘* 

177.63 

(1+4504/*)^/* 

26.87. 

28800/ 

(1+3128/*)^/’ 

1.38 X 10"* 

213.14 

35.76 

28800/ 

1.85 X 10“* 

284.23 

a + ns9fh*^ 

44.70 

28800/ 

2.31 X 10“* 

355.25 

(1 ♦ 1126/*)^/* 

53.64 

28800/ 

(1 + 782/*)^/’ 

2.77 X 10"* 

426.29 




I H HlgiiffWi i»ga 8 Pap ^ ' www > «JMrwy i iiw i w^wij^ ' ■*>■" ‘ 1 —nm i n .! . ^ 


TaM»a. WMtoniwtorMiMmiltiws* 


y^z’ *"/* 

T^* Nm 

4.47 

1.09 X 10‘ 

8.94 

4.37 X 10‘ 

13.41 

9.84 X 10‘ 
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S-Band Maser Phase Delay Stability Tesfs 

J. M. Ursch, F. Alcazar, J. Galvez, and A. Rius 
Madrid Deep Space Station 

C. A. Greenhall 
TDA Enginaering Office 


This article summarizes the results of the S-band travding-wave maser phase delay 
stability measuremmts performed at DSS62, These tests were required for the 
Fioneer-Verm wind experiment 



I. Introduction 

This article summaiizes the results of the S*band traveling- 
wave maser phase delay stability measurements performed at 
DSS 62. In conjunction with earlier group delay measure- 
ments, this article completes the group and phase delay 
measurements required for the Pioneer-Venus 78 wind experi- 
ment. 

All tests were performed for the entire S-band maser 
passband» which includes the 5*MHz passband, from 2290 to 
2295 MHz, used in the entry probe wind experiment. This 
experiment requires that the phase delay variation be less than 
1 deg over a 2-MHz passband (2291 to 2293 MHzj. 

II. rest Configuration 

All the tests were performed on the declination (dec) room 
maser. Hgure 1 shows the physical location of test equipment. 
As can be seen, most of tliis was installed in the empty part of 
receiver rack 9. 

Figure 2 is the test configuration schematic. To have a 
phase delay difference of less than 180 deg between test and 


reference paths, a delay line was fabricated locally. This 
consisted of about 38 meters of soft line encapsulated in a 
plastic box and insulated thermally. 


IN. Pfuwe Delay Stability VorausMaoor 
Paratne^o 

Phase delay varhtion over a given frequency band is defined 
as max AdCT) ~ min Mif)> vdiere is the phase delay 
versus frequency curve, &4(f) is the change in d(/) due to 
some internal or environmental change, and the max and min 
are over all /in the band. Notice that if the (if) curve simply 
moves parallel to itself, the phase delay variation is zero, since 
Ad(/) is constant. If the phases of the signals from all 
Floneer-Venus probes are shifted by the same amount, no 
harm is done becav'e the phases are differenced 

Variations of five maser parameters were simulated eiti ploy- 
ing the same method used in the group delay stabdity 
mejsurements (tee Table 1). 

The expected variations were taken from the technical 
manuab of the corresponding power supplies and checked 
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against actual variations for periods of about 24 houn. 'fhe 
factor is the ratio between the simulated and the 
expected parameter variations. This must be taken into 
account when estimating the expected phase delay stability. 


t*mipeiatures. However, based on the coherency of results 
obtained, it is possible to derive from £q. 1 the phase delay 
stability lA^| ~ ^ ^ group delay stability Ar 

measuted earlier. See Table 3. 


IV. stability Varaus Antenna Ortentallon 


A. RmuNs 

Rguies 3 duoujji 7 are copies of the actual curves obtained 
while peifcxining the above tests. Table 2 summarizes the test 
results. 

As shown, two sets of measurements were performed 
(RgS. 3 dirough 7 correspond to set No. 2) with the dispersion 
between both being 1.5 deg. This is more than acceptable 
taking into account the different factors involved. 

Hie variations were measured for a passLand of 40 MHz 
(2265 to 2305 MHz) as well as foi passband of interest in 

the probe wind experiment (2290 tc 2295 MHz). They should 
be ^vided by the correqionding scale factor to obtain the 
expected variaticais. 

B. Cotwww wWh Group Dwiay StabiNy H — uK a 

B> calculating the group delay stability from the data 
obtained for phase delay stability and comparing it with the 
data derived from measuring the group del^ stability directly, 
a reasonable coherency may be obtained between both group 
and phase delay stability measurements. The relationship 
between group delay stability and phase delay stability is given 
by 

where 

Ar = group delay stabilit>' from phase delay 
stability data 

B - width of passband 

A^i andA ^2 ~ phase delay differences in both extremes 
of the passband, considering their signs 

The value of ^ may be derived for all parameter variations 
and compared with Ar, the group delay stability measured 
earlier. Figure 8 details the result of this comparison for the 
40-MHz passband. It can be seen that tlie difference is less 
than 0.28 ns except for bandwidth current, ^ere the differ- 
ence is less than 1.07 ns. 

C. EtUiiialionofPliaMlMtyStibi 
TMiptriturM 

Due to the heavy tracking load, it was not possible to check 
the phase delay stability versus dec and control room 


The appic^riate range of antenna movement to simulate 
the posstble earth magnetic field variation (due to DSS 43 
antenna position changes) during the entry probe wind 
experiment was estimated using the earth magnetic field model 
of Lang (Ref. 1): 

^ (2m cos 6)/R radial component 
Bq^ (m sin 6)IR^ local meridian component 


From the above we found that the magnitude of the earth's 
magnetic field for DSS 62 is 0.47 gauss and its elevation angle 
is 59 deg. For DSS 43 the magnitude is 0.43 gauss. Therefore, 
by moving the DSS 62 antenna from the az = 0 deg, el - 59 
deg position (HA = 0 deg, dec - 72 deg), where the maser 
magnetic field is perpendicular to the earth's magnetic field, to 
az = 180 deg, el - 40 deg (HA - 0 deg, dec == 340 deg), where 
both magnetic fields are parallel, it is possible to simulate a 
variation larger than that expected in the PV wind experiment. 
In addition, it should be noted that the earth's magnetic field 
at DSS 62 is slightly greater than at DSS 43. 

Although during the actual test the anteima could only be 
moved from HA = 0 deg, dec - 70 deg to HA - 0 deg, dec = 
360 deg due to safety reasons (personnel and test equipment 
in dec room), the above reasoning is still valid. 


Figure 9 illustrates the stability versus anteima orientation 
results. The antenna movement sequence was: Zenith (Z) ^ 
dec 70 deg Zenith (Z") -► dec 360 deg -► Zenith (Z"). 

It should be noted that it is quite impossible to tie down all 
the test equipment and cables so that they do net move during 
the antenna movement. Because of this, after multiple tests 


B^- 0 local parallel compcxient 

v^iere 

m = 8.1 X 10"** gauss X ern^ 

R = 6.37 X 10* cm 

^ = 50 deg for DSS 62 and 125 deg for DSS 43 


A. Rwultt 
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that proved that for very small equipment and cable move- 
ments the phase delay may vary about 5 deg, we conclude that 
the differences between the curves shown in Fig. 9 for all 
zenith and dec 70-deg positions may be due to slight 
movements of the set-up components. 

From the above it may be asserted that phase delay 
stability is better than 25 deg for the entire maser passband, 
40 MHz, and better than 5 deg for the 5-MHz passband of 
interest here. 


The reason for the different shape of the curves in Figs. 3 
through 7 and the curve in Fig. 9 was a maser klystron change 
in between. 


V. Long-Term stability 

The long-term stability of the maser phase delay was 
checked during periods of 8 hours. The test conOguration was 
the one detailed in Fig. 2 except for the X-Y plotter. This was 
replaced by a strip chart recorder HP 7100. 

During the warmup period of the test equipment, the phase 
delay drifted-considexably.* Moreover; for more than- 6 hours 
after this period, the phase delay drifted continuously at a rate 
of 2 deg per hour. 

To check that the phase delay drift observed during the 
first hour and a half was due to the test equipment warming 
up, we performed an additional test in which we replaced the 
maser by a soft line of the proper length. Under these 
conditions, all the drift should be due only to the test 
equipment (cables included). Results showed the same consid- 
erable drift during the first hour and a half. For about 6 hours 
after that, no appreciable drift was observed. 

The conclusion of these long-term stability tests is that a 
phase delay drift of about 2 deg per hour can be expected. 

VI. Conclusions 

A. Qtntral 

Conddeiing the scale factors, one may infer from Table 2 
that the parametric variation yielding the largest phase delay 
instability (in the 5>Mriz passband) is the liquid HE tempera* 
ture (vapor pressure gauge). The same result was obtained 
when measuring the group delay stability. Nevertlieless, this 
phase instability is only about 0.16 deg. The root-sum-square 
(RSS) of the phase instabilities due to all parametric variations 
is 0.17 deg. 
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As far as room temperatures are concerned (Table 3), that 
of the dec room is the one yielding the highest phase delay 
instability, about 0.27 deg. 

The RSS of all phase instabilities from Tables 2 and 3 is 
0.32 deg, and the linear sum is 0.58 deg. 

From the antenna movement tests ;Table 4) it may be seen 
that the instability in the 5-MHz passband is better than 5 deg. 
It should be remembered that the simulated antenna move- 
ment represents almost a complete range of interaction 
between the maser and earth magnetic fields (from perpendic- 
ular to parallel portion). The actual antenna movement during 
the 90-minute wind experiment will be 22.5 deg in hour angle. 
Hence, the change in the angle a between the magnetic fields 
of the eai Ji and the maser is less than 22.5 deg. Assuming that 
the phase delay variation V is proportional to the change in 
cos a we have 


K<5deg(22.5j^)<2de£ 


This bound corresponds to a 4^ error in the difference of 
probe, position components, over 90 minutes (0.74 deg over 
1000 seconds) measured by the Goldstone-Canberra baseline. 
The requirement is for < 1 deg over 1000 seconds, ^ch is 
met. In addition, the wind experiment is intended to measure 
probe velocities. Since error Av in the difference of probe 
velocity components is proportional to dV\dt^ the largest that 
Av can get is 4 m/5400 s = 0.00074 m/s, a negligible error. 
The position error is gradually accumulated over the whole 
90-minute period. The same reasoning applies to the long-term 
drift discussed in Section V. 

B. Curv« Analysis 

When comparing Figs. 3 through 5 with Figs. 6 through 7 it 
was observed that, while the curves in the first figures were 
rather parallel throughout the passband, the curvesin ? igs. 6 
and 7 appear to '^rotate** (or change slopes) near the DSN 
bandwidth (2285 to 2305 MHz) around which the maser 
klystron is tuned. This may be explained by using the Ijeflang 
model for the phase-versus-frequency maser characteristic, 
given in degrees by (Ref. 2) 


6.56 2 G, 




\^ere G/ and ff are the gain in dB and the central frequency of 
tire section, and Aft is the maser material resonance line 
widtli. 
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Assuming »/- /, we have 


#(/) » 6.56 
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bandwidth control current) will yield rather parallel curves. On 
tire other hand, those affecting such as gain control 

«^oltage and cryogenic temperature (vapor pressure gauge), will 
^ve curves with different slopes and rotating around a zone 
determined by the klystron tuning. 


To first approximation. 
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This expression is a stiai^t line with its slope proportional 
tc Gi + (7, and the intersection with the frequency axis being 
the gahi'Weighted mean of the central frequencies of the two 
sectirms. Hence, those parametric variations that mainly affect 
the /, and/, position (such as beam and reflector vdtages and 


As for the shape of curves in Fig. 9, the predominant effect 
for the dec 360^g curve is the total gain s^iile for the dec 
70-deg curve it is the/, ,/, positions. 


VII. Instability due to Cross-Head 
Modulation 

During one set of measurements, the cross-head modulation 
effect appeared for a short period of time and has been 
reflected in Fig. 10. This instability has the same shape as that 
observed vidiile measuring the group delay stablUty. The 
stability may be estimated at less than 5 deg and may be of 
some importance, if present. The shape of this curve can also 
be explained due to the rapid and opposite gain variations 
caused by the cross-head modulaticm effect. 
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Parameter 

Simulated 

variation 

Expected 

variation 

Scale 

factor 

Beam voltage 

±SV 

< 0.75 V 

13.4 

Reflected voltage 

±2V 

<ai V 

40 

Bandwidth control current 

± ImA 

<50 A 

40 

Gain control current 

t 3 mA 

<50A 

120 

Liquid helium temperature 

1.3 N/cm* » 
21b/to*(VPC) 

0.4 N/cm^ « 
0.58 Ib/in^ 

3.2 



Tabto^ RMuto 

OlpMIMM 
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Parameter 

Simulated 

variation 

Measured phase delay 
variation for 
a 40>MKz passband, deg 

Set 1 Set 2 

Measured phase delay 
variation for 
a 5-MHz passband, deg 

Set 1 Set 2 

Scale 

factor 

Beam voltage 

±5 V 

7 

6 

<0.5 

0.5 

13.4 

Reflected voltage 

± 2 V 

9 

8.5 

2 

0.5 

40 

BaiHlwidth control current 

± 1 mA 

2 

2 

a5 

1 

40 

Gain control current 

1 3 mA 

2 

Z5 

<0.5 

<0.5 

120 

Liquid helium temperature 

1.3 nW 

5 

5 

<0.5 

0.5 

3.2 


■> 2 Ib/to* 
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Simulated 

Expected 

Peak-to-peak 

Estimated phase 

Parameter 

variation 

variation 

variation of 
group delay, ns 

delay stability, deg 
B-40 MHz B«5 MHz 

Dec room temperature 

i.7®C 

1.7®C 

0.3 

11 

a27 

Control room temperature 

1.7«C 

0.6®C 

0.2 

1.4 

ai8 




u 



Measured - - 

Pusband. MHz 

phase delay 


stability, deg 

40 MHz (2265-2305 MHz) 

Better than 25 

5 MHz (2290-2295 MHz) 

Better than 5 
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A Preliminary Model for High-Power Waveguide Arcing 

and Arc Protection 

H. C. Yen 

Radio Frequancy and Microwave Subey stems Section 


This is the first article in a series describing the ongoing effort of the Transmitter 
Group to upgrade the arc protection subsystems that are, or will be, implemented in the 
DSN high-power transmitters This article review the status of our present knowledge 
about waveguide arcs in terms of a simple engineering model and discusses a fairly general 
arc detection scheme. Areas where further studies are needed are pointed out along with 
our proposed approaches to the solutions of these problems. 


I- Introduction 

This is the initial report on one phase of the ongoing effort 
of the Transmitter Group to upgrade the arc protection 
subsystems currently used in the DSN hi^-power transmitters. 
The ultimate goal of this effort is to develop, with a known 
margin of safety, an extremely fast and highly reliable arc 
protection system wdiich can be designed, manufactured and 
implemented for the eventual automation of these transmitter 
subsystems. 

High-power microwave breakdown invariably results in 
arcing inside a nitrogen-filled vraveguide system (P~l atm). 
The breakdown and subsequent arc formation interrupts the 
normal hij^-power transmission and causes a very large 
standing wave to b< set up between the microwave source and 
the arc. Such high-level standing waves can be quite destructive 
to the components employed in the hi^-power transmission if 
their peak power ratings are exceeded. In addition, the 
accompanying am once formed was found always moving 
toward the source at a speed depending on the field level 
(Reft. 1,2,3). Such high-power arcs can cause new or addi- 


tional damages through their heating and thermal shocking 
effects as they propagate through the system. Often the arc 
damage is of a very severe nature. As an example, if the arc is 
allowed to travel to reach the output window of a hi^-power 
klystron amplifier before the klystron output can be turned 
oft*, the arc literally hangs up at the window and destroys it, 
which usually leads to a costly destruction of the klystron 
amplifier itself. 

Even though precautions can be taken to reduce the chance 
of microwave breakdown within the designed power rating 
vdien assembling the transmission system (Refs. 4, S), a 
complete prevention can never be assured. To protect system 
components against damages in the event of breakdown, 
protective devices must be employed. Protection is accom- 
plished by removing or reducing ^ microwave energy that 
feeds the breakdown within a very short period of time after 
the breakdown is detected. For a klystron amplifier, this can 
be the turning oK or drastic reduction in the Idystron drive. 
Since the energy level required to sustain an arc is much 
smaller tha.i that to initiate one, once an arc is initiated, the 
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klystron output must be well attenuated to completely 
extinguish the arc. 

The arc and the arc plasma cause hig|i reflected power 
during their formation and propagation (Refs. 1,2,3). A 
reverse power detector, properly located, is a hi^y satisfac- 
tory device for arc detection. However, the reverse power 
detector doesn't provide protection against arcs formed 
between the reverse power detector (directional coupler) and 
the source - the most damaging land. For a positive arc 
protection, other types of devices such as li^t-sensitjve ones 
have to be used in this region, preferably with one device 
viewing the klystron output window directly. 

In the following, we will concern outselves mainly with the 
protection of the output window against arc damage, as arcs 
formed else^ere can be easily detected by monitoring the 
reverse power. A simple model diowing the energy relationship 
associated with microwave arcing and a general optical 
detection scheme is discussed. Due to the complexity of the 
process and lack of applicable experimental data, only a 
semiempirical discussion is attempted at the present time. This 
article will briefly review the status of our present knowledge 
about miCTOwave arcing inside a high-power waveguide system; 
at the same time it will also point out the areas where further 
theoretical and experimental studies are needed as well as our 
proposed approaches to the sdutions of these problems. 

II. Arc Model 

Radiations emitted by nitrogen arc plasma have been 
actively studied both experimentally and theoretically (Refs. 
6,7). Because of the complexity of the arcing process, a 
general theoretical treatment is impossible. Often, simplifying 
assumptions such as local thermodynamical equilibrium, the 
arc thermal-chemical state, arc parameters (pressure, tempera- 
ture, current density, radius, symmetry, etc.) and radiation 
wavelength range have to be made to facilitate the calculations 
of both line and continuum radiations. Experimentally, the arc 
plasma has to be stabilized and stationary for data collection. 
Obviously the applicability of the results from these studies to 
the case of waveguide arcing during high-power transmission is 
not without question, because the latter is essentially a 
time-dependent discharge phenomenon. As shown in these 
studies, arc plasma radiation depends strongly on the thermal- 
chemical state of the plasma, which in turn varies considerably 
with the degree of departure from equilibrium. To make the 
matter more uncertain, these studies often show some degree 
of discrepancy among themselves regarding the interpretation 
of experimental results. Nevertheless, they do provide us an 
important means to understand arc radiation on a microscopic 
level. The model given below will be based on macroscopic 
processes for a first-order characterization of the waveguide 


arc, however. The reason is that the macroscopic characteris- 
tics are only partially understood so far, and a comfdete 
knowledge is needed for any meaningful study of the 
underlying microscopic processes. 

According to the law of energy conservation, the energy 
relationship right before and after the waveguide arcing can be 
represented by the block diagram in Fig. 1. 

Experimentally, the characteristic tiiTie for complete arc 
formation has been found to be of the order of a few 
microseconds. Conceivably, this time must be a function of 
the power level. Unfortunately, to our knowledge, this 
functional dependence has not been systematically examined 
to date. It turns out that the arc protecticn also has to be 
accomplished within approximately the same period of time 
after initial arcing is detected in order to avoid any real 
damage (Refs. 1,8). This microsecond response time is one of 
the important design parameters for an arc detector. The arc 
shape based on the field profile (TEjq mode) would be a 
cylinder with an approximately ellipti^ cross section. The 
actual shape has not been reported previously, probably 
because of the difficulty of direct visual observation of 
waveguide arcs. In some experiments where continuous arc 
tracks were left behind, a spatial periodicity of half guide 
wavelength along the propagation direction can be seen, with 
the spatial extent in the cross-sectional plane being modulated. 
It seems possible to have simultaneous arcings at half guide 
wavelength apart after the initial arcing if the standing wave 
field exceeds the breakdown field. This possibility has not, to 
our knowledge, been explicitly stated or observed in the 
existing literature. 

Once fully fotmed, the arc plasma becomes a very effective 
shorting plane despite its limited extension in space, as 
evidenced by the presence of very little transmitted power and 
a very large amount of reverse power. Typically, the trans- 
mitted power is about 1 percent of the incident power, while 
the reverse power can range from 50 percent (Ptn^^O kW) to 
almost 98 percent (Pj„^ 1 MW) of the incident power (Refs. 
1 , 2 , 8 ). 

The difference between the incident and the reflected plus 
the transmitted power Is the power associated with arcing 
(energy ’^absorbed" by the arc). Although not all of this 
energy is converted into radiation, it is believed that a 
substantial fraction is. A division of the energy according to 
the observable macroscopic i^enomena is given in Fig. 2. Such 
division is somevriiat simplified and not without ambiguity. It 
is done here for the purpose of identifying major phenomena 
which might offer a potential for arc detection. A brief 
discussion of each category is given as follows: 

(1) Sonic wave (pressure wave) is produced along with the 
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breakdown. At high power level, the sonic wave 
becomes loud enough to be audible, and often causes 
the waveguide to shake. It is created in a short burst as 
breakdown is initiated. Conceivably we can use acous- 
tic transducers to pickup such sonic wave as an indirect 
indication of arcing. Signab are then processed to 
activate the protecti(m mechanism. But in view of the 
requirement on the device sensitivity, stability, 
response time, power supply and noise immunity, this 
method seems at best a supplementary one. 

(2) Owing to its surface resistance, the waveguide tempera- 
ture is elevated during power transmission. At equili- 
brium, the temperature is stabilized and the gas inside 
will also come to equilibrium at this temperature. 
Because of the reduction in gas dielectric strength at 
hi^er temperature for a given pressure, the maximum 
power handling capability of the waveguide system is 
lowered. If a breakdown is suddenly initiated at a 
certain point, additional heat is imparted to tfie gas 
molecules surrounding that point, with the effect of 
reducing local dielectric strength even further. With the 
incident power maintained at the same level, it is clear 
that a catastrophic chain breakdown can be triggered. 
This is one of the mechanisms that explains the arc 
movement toward the source. One way to increase the 
ihreshold of arcing is to provide sufflcient cooling 
capability to the waveguide system to remove the 
dissipated heat and to prevent formation of hot spots, 
thus preserving the maximum power rating. The other 
would be to increase the gas pressure, which is limited 
by the strength of the window material and the 
tolerable leakage rate only. 

(3) The waveguide wall heating discussed here is the 
additional heating due to the presence of arcs. During 
arcing, part of the energy is deposited in the waveguide 
wall in the form of heat through the bombardment of 
energetic electrons and ions. For a prolonged arcing, 
the local temperature rise in the waveguide can be such 
that the wall material becomes melted. The turbulent 
nature of arcing can eject small masses of the molten 
material into the waveguide. These ejected masses are 
not directly thermally anchored and get heated up very 
rapidly, causing more arcings because of the reduced 
dielectric strength and local field enhancement around 
these masses. This situation is known as massive 
breakdown and fortunately is more likely to nappen at 
much hi^er power leveb. In any event, arcs once 
initiated have to be quenched as soon as possible to 
avoid any prolonged arcing and its extensive damage. 

(4) The remaining part of the energy is distributed among 
various atomic and molecular processes such as dissoci- 
ation, ionization, reattachment, recombination or even 


chemical reaction with wall material Some cf them are 
radiative, and the others are nonradiative, depending on 
whether photons are released as a product of the 
process. Hie radiation emitted consists of line spectrum 
and continuum. Because the arc is in a highly complex 
state with many chemical species, reabsorption of the 
emission is not ncj^gible, especially for the line 
spectrum. Therefore, what we are interested in is the 
total amount of energy of the net emission coming out 
of the arcing region and how it is spectrally distributed. 
Sucli information for waveguide arcs is rather incom- 
plete and is very difficult to obtain due to its 
time-dependent and random nature. There exist some 
calculations of equilibrium nitrogen plasma radiation 
which agree more or less with experimental results 
(e.g.. Ref. 6). As indicated before, the radiation 
characteristics depend stron^y on the thernial-chemical 
state of the arc plasma. For example, for a difference in 
plasma temperature by a factor of 2 (say 7,000 and 
14,000 K), calculations showed that the net emission 
per unit source volume per unit solid angle can differ 
by more than 3 orden of magnitude. Moreover, these 
calculations covered only a very limited range of 
radiation spectrum while experimental observations 
seemed to indicate a much broader spectral distribu- 
tion, probably ranging from UV to low frequency RF. 
Thus it shows the need that more data pertinent to the 
real arcs have to be gathered experimentally. 

In order to overcome this inherent difficulty' of measuring 
time-dependent arc properties, we plan to carry out measure- 
ments on a sustained arc inside a resonant section of a 
waveguide system. This scheme will not only provide us a 
more or less stable source but needs only lIQ of the power to 
maintain the same level inside the resonant section, where Q is 
the loaded quality factor of the resonant section. This should 
solve the dilemma of trying to initiate an arc on one hand and 
on the other, to protect the klystron amplifier at the same 
time. Both the total radiation intensity and the spectral 
distribution from the sustained arc will be measured. The 
spectral range covered will be determined by the measuring 
instrument. With known arc geometry, detector sensitivity and 
the relative orientation of the arc and the detector, a 
normalized arc emission characteristic can be obtained. Such a 
controlled experimental arc probably is a good approximation 
to an arc hanging at the klystron output window and is 
expected to less resemble a moving arc. The variable param- 
eters in this study would be incident power level, gas pressure, 
gas temperature, arc geometry and pulse length. With sli^t 
modification, the arc formation time, the arc energy absorp- 
tion and arc traveling speed can also be measured. Detail 
experimental considerations, arrangements and results will be 
presented in the future. 
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I III. Arc Detection 

r 

With the characteristics of waveguide arcs approximately 
I known, we could then proceed to design a detection system 

^ with certain predictable performance. Because the detected arc 

I radiation is eventually converted back to electrical energy for 

i control purposes, a fairly general detection scheme can be 

fv outlined as shown in Fig. 3. In addition to the direct arc 

I emissions and those indirectly scattered into the optical 

I entrance aperture, there are also ambient li^t leakage and 

I output window glow. The ambient light can, in principle, be 

reduced to an insignificant level by properly ananging and 
' shielding the microwave components. Ilie window glow on the 

^ other hand presents potenti^y a more serious noise source. 

? The gjiow is a consequence of the dissipation in the output 

I window under normal hi^ power transmission; i.e., it is an 

! integral part of the power transmission process. Its exact 

I nature has not been carefully examined so far, and its effect 

I on the arc detector can only be speculated at this time. Our 

I study will look into this problem to determine its seriousness 

{ and will recommend remedies, if needed, such as better 

I cooling of the window or designing a frequency-sensitive 

detector which will ignore the ^ow. 


T{v) = transmittance of the optical fiber guide 

F(v) = transfer function of the optical signal process- 
ing 

R(ii) - responsitivity of the optical transducer, (A/W) 

Then the optical detector partial electric cunent output due to 
the arc light input at the entrance aperture is 

yi/.w+w+w 

+/ 4 WI mP(p)R(p)dv (2) 

S /j +/j +/^ 

vAim signal current comes from direct and indirect arc 
emissions (7j +/ 2 ). 

The total current of tlie optical transducer is 



Upon entering the aperture, not all the light is transmitted 
due to the Fresnel reflection if there is a difference in the 
indices of refraction across the interface (typically 5 percent). 
Furthermore, if optical fiber is used for transmitting the light 
signal over some distance, only those beams with incident 
angle less than a certain acceptance angle will be launched into 
propagating modes of the optical fiber. 

If we define /| as the spectral density at frequency v 
of arc emission per unit volume per unit solid angle at arc 
spatial coordinates % then 

'i0>) = I Jiiv^dVd- (1) 

I volume 

is the spectral density of the direct arc emission entering the 
entrance aperture that will propagate inside the optical fiber 
guide, nidiere the solid an^e H is subtended by the entrance 
aperture toward the arc. Similarly, we can define 
and J^iy) for the indirectly scattered arc emission, the amUent 
light leakage and the window glow respectively. 

Furthermore, we will define the following quantities: 

■ fraction of Fresnel reflection loss, where »| and 
n, are the indices of refraction for the two 
adjacent media across which the light is tran^ 
mitted. 


^total 


= (/l+/2) + (/.+/4+/j^,) (3) 


where is the noise current of the optical transducer. A 
substantial fraction of the noise current comes from the 
temperature-sensitive dark current of the detector. 

The signal to noise ratio (SNR) at this point is 

, 4 ) 

if ambient light leakage is negligibie. 

The magnitude of this required input SNR is an important 
design parameter to be determined. The final SNR at the 
output of the threshold logic circuit depends, of course, on the 
electrical components used and should be maximized despite 
the presence of power supply ripple, EMI, and thermal drift. 
Its magnitude will be determined based on the required signal 
level, detection probability and acceptable false alarm rate. 

In reality, there are additional losses in the optical path 
when connections have to be made. A general rule of thumb 
here is that about 2*3 dB/connector and 0.1 dB/splice loss are 
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present. For each step of optical path, ambient light leakage is 
also always present. We have assumed they can be minimized 
to be neghgible. Typical fiber transmittance over the visible 
and IR spectrum is about 3dB for a length of 0.5 m. 
Additional fiber input coupling losses may be necessary in 
order to isolate the optical transducer from being exposed to 
intense RF. The dielectric fiber may be advantageous in 
overcoming ground loop problems in the detection circuit, 
however. 

The eventual total cuncnt output is limited by the spectral 
range over which the transducer has appreciable responsitivity. 
The coimnoniy used optical transducers are pin photodiode, 
avalanche diode and photomultiplier tube. The first two are 
solid-state devices which require a bias of about 20 and 200 
volts respectively. The photomultiplier tube is rather bulky 
and requires more than 1 kV supply. Our experimental study 
will primarily concern solid-state devices because, in parallel to 
this study, we will also evaluate a commercial arc detector 
based on photomultiplier tube. Solid-state devices seem to 
have adequate sensitivity and response time for arc detection. 
But their dark currents tend to increase rather rapidly with 
increasing temperature, thus changing the operation point of 
the arc detector. Often this thermal drift is large enough to 
cause false activation of arc protection and to saturate the arc 
detector. These have been the major problems encountered in 
the present arc detectors. The increase in the dark current as a 
function of temperature is inherent to the semiconductor p-n 


junction inside the diode and, therefore, the solution b cither 
to maintain the diode temperature, or to compensate the 
thermal drift or to use ac coupling and/or high-pass filtering 
for the signal and noise. 

The output voltage originated from the optical transducer is 
then compared to a reference voltage to produce a logic 
output. The threshold is determined by the amount of light 
output expected from an arc as well as the safety margin for 
noise and thermal drift. One way to increase the detector 
sensitivity is to increase the load impedance of the optical 
transducer as much as possible. However, the most important 
consideration here is to make sure that the time delay between 
the arcing and the logic pulse output is less than a few 
microseconds for a safe arc protection. This may have to be 
done at slight expense of the detector sensitivity and SNR. 

IV. Conclusions 

We have briefly discussed the waveguide arc according to a 
simple model based on an engineering viewpoint and outlined 
a general scheme for achieving protection against arc damage. 
Several questions have been raised concerning the waveguide 
arc characteristics and the arc protection that we plan to 
answer both experimentally and theoretically as soon as 
possible so that a reliable arc protection system becomes 
available not only to presently attended transmitter subsys- 
tems but also to unattended ones in the future as well. 
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In order to improve the performance and capabilities of the DSN 64~m antennas at 
X-band, extensive modifications to the XRO cone assemblies are in process. The changes 
include a new feed assembly with a dual hybrid mode horn and orthogonal mode 
junction^ dual traveling wave masers, and a new receiver mode selector. 


I. Introduction 

Tlie X-band cone assembly (Ref. 1) on the 64-meter 
antennas ori^naUy provided the capability for RCP reception 
only, with one traveling-wave maser (TWM). The feed was later 
changed to a Mod 11 version (Ref. 2) that provided remotely 
selectable RCP or LCP. 

In order to provide increased performance and additional 
capabilities, especially for the Voyager mission, the XRO cone 
uidergone a complete redesign. It is planned to reconfig- 
.. (he three DSN cones in the field; this is the first time such 
an extensive change to a feed cone has been attempted in 
place. 

After removal of soirte 90 peipent of existing equipment 
within the cone, the following mgjor chatig s will be imple- 
mented: 

(1) A new feed assembly, the Mod III XRO feed. 

(2) Dual X-band traveling wave masers. 

(3) New monitor receiver. 

(4) New TWM support racks. 


(5) New output waveguide runs. 

(6) New output switching configuration for the monitor 
receiver. 

(7) New receiver mode selector. 

To further complicate the changes, the two oversea; 
stations will be implemented with the new MidU! feed 
assembly on a different scheduie from DSS 14. Ihus it was 
necessary to design two new XRO cone coraTgurations: (1) the 
Mod II version for overseas serves as an interim configuration 
until items (1) and (7) abovu are implemented; and (2) the 
Mod 111 version for DS3 14 and, later, DSS 43 and DSS 63. 
Because of the differences in physical size and output 
configuration of the Mod II and Mod III feeds, this required 
extensive differences in such details as the maser support 
stands, maser output waveguides, twitch contrd nib^stem, 
and many support brackets and Iiardware. 

a. ThuMocfllSXROFMAMMnbly 

The new feed assembly (Fig. 1) provides the fcdlowing 
m^ changes: 


m 



(1) A dual hybrid mode feed hom (Ref. 3) for more 
efficient subreflector illumination. 

(2) An orthogonal mode junction to permit the simultane- 
ous reception of RCP and LCP. 

As discussed in Ref. 3, it is expected that the new feed hom 
wi^i increase the antenna gain approximately 0.29 dB. 


III. Dual Traveling Wave Masers 

The primary purpose of implementing dual masers is to 
improve reliability during critical encounter periods by having 
a second maser cooled down and ready to be switched on line 
in seconds should the primary maser fail. Tne presence of two 
operating masers, however, also leads to increased radio 
science capability. Thus, with the orthogonal mode junction, it 
will be *x)ssible to receive two orthogonal polarizations 
simultaneously. Since the polarization selectability of the 
Mod II feed is retained in this new feed, it is possible to 
remotely select which polarization (RCP or LCP) is received 
by TWM 1 , with TWM 2 receiving the opposite polarization at 
all times. 


IV. Receiver Mode Selector 

Because of this dual-channel feature, it was necessary to 
design a new receiver mode selector that provided two inputs 
(rather than a single input as before) an^ four receiver outputs, 
so configured that any receiver can independently connect to 
either maser output. Up to four additional receivers (for a 
total of eight) can also be connected to the mode selector 
panel, but witli some restrictions on independent maser output 
selection. 


V. Configuration Control Subsystem 

Concurrently with the changes to the XRO cone assembly, 
a new configuration control subsystem is being implemented. 
This necessitated additional changes in the cone switch control 
junction box, cable package, and control panel. 

VI. Maser Support Equipment and Output 
Waveguides 

With the removal from the cone of the older configuration 
maser package and monitor receive*^, it was necessary to 
relocate the output waveguide luns and design new brackets 
for the new maser packages and their support racks. This is the 
second in»v.-iation of this equipment after the SXD cone 
assemblies the 34-m antennas (Ref. 4), and some of the 
techniques generated in that project were adapted to the XRO 
cone redesign. A complicating factor in laying out the 
waveguide runs involves the requirements for four waveguide 
couplers. One coupler for each maser output supplies tlie 
monitor receiver with its signal tlirougfi a waveguide switch to 
select TWM 1 or 2. An additional coupler for each maser 
output provides an injection point for calibration signals. 

VII. Fieid Implementation 

At the time of this reporting, DSS 14 has been imple- 
mented with the Mod III configurations, with only a few 
minor details remaining to complete the task. Test data are not 
yet available. 

During October and November the overseas stations will be 
implemented with the Mod II XRO cone configurations (dual 
mast**s only). It is planned to upgrade these stations to tne 
Mod III configuration during late 1 979, in time for Saturn 
encounter. 
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A water vapor radiometer (WVR) has been developed that measures the atmospheric i 

noise temperature at two different frequencies near 22 GHz, These noise temperatures 
are used in empirical-theoretical equations that yield tropospheric range delay,, in centi- 
meters, through the atmosphere along the beam of the H'KR. This range correction is then - 

applied, as needed, to measurements concerning spacecraft range and to VLBI baseline \ 

determinations. This report discusses the WVR design and calibration techniques, j 


r 



I. Introduction 

Development work has continued on the water vapor 
radiometer (WVR). This article wiU discuss progress since the 
last reports (Refs. 1 and 2). Also mentioned are new frequen- 
cies that resulted from analytical support activities and a dis- 
cussion of the calibration load techniques that have evolved 
since this work started. The cold calibration load mentioned 
in Refs. 1 and 2 has been replaced by a hot load. The block 
diagram of the instrument is unchanged since Ref. 1 except for 
this change in load design. Figure 1 shows the WVR in its 
current configuration. 

A formal set of technical specifications and recommenda- 
tions for WVRs has been generated (Ref. 3), with inputs from 
the outside scientific community. 

II. Calibration Loads and Techniquos 

In the water vapor radiometer, the radiometric noise tem- 
peratures of the sl^ at two frequencies near the water vapor 
resonance line, 22.235 GHz, are used to infer the phase delay 


through the atmosphere. To get an accurate measurement of 
these noise temperatures, the radiometer can be accurately 
calibrated by using waveguide calibration loads adiose absolute 
radiometric temperature has been determined from iasertion 
loss and physical temperature measurements (Ref. 4). 

The method of using absolutely calibrated waveguide loads 
to calibrate a radiometer depends on several assumptions: 

(1) The RayleighJeans approximation to the Planck 
blackbody-radiation law must hold at the frequency of 
interest. 

(2) Ihermal gradients between the physical temperature 
sensor and the passive microwave load radiating region 
are negligible. 

(3) Energy radiated by the load is not significantly 
reflected in reaching the radiometer receiver. 

(4) Insertion losses and their physical temperatures 
between the load and the radiometer receiver are 
known. 
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(5) The temperature sensor that measures the physical load 
temperature must be of sufficient accuracy (Ref. 5). 

It is the aggregate effect of all these uncertainties (Ref. 6 iden- 
tifies ten of them) that leads to the total cai-load temperature 
error. It should be noted that an ambient temperature cal- 
load can be designed to have no significant temperature 
gradients. In this case, the radiometric temperature of the 
ambient load is equal to its physical temperature (the load 
looks like a perfect blackbody): 


where for future reference, 

T is temperature in kelvins 
A is ambient 
R is radiometric 
P is physical 


( 1 ) 


The best radiometer measurements are achieved when the 
instrument calibration curve extends on either side of the 
desired measurement re^on, ie., data can be interpolated 
rather than extrapolated. In the case of the WVR, measured 
sky temperatures arc as low as 10 K; it is impractical to try 
to achieve the true interpolation condition for the calibration 
curve. However, having one calibration reference load near 
the coldest expected values would still be desirable to mini- 
mize the strong “leverage” effects that often result when using 
an extrapolated calibration curve. 


To achieve these desired colder temperatures, several 
refrigeration devices were tried. The first was a two-stage 
Peltier junction, which was quite reliable but achieved only 
“40®C (233 K), and then only with a two-gallon-per-minute 
total4o$s, water-cooled heat sink. The second device was a 
reciprocating -expansion refrigerator that achieved -90^C 
(183 K), but required, aiong with its air source, Wgh 
maintenance. 


If the hot cal-load were to be absolutely calibrated, all of 
the previously mer ed uncertainties in its noise output 
would be multiplied by the leverage effect of having both 
cal4oad$ far removed from the sky temperature. For example, 
if the two cal4oads and sky are 400 K, 300 K, and 20 K, 
respectively, cal4oad uncertainties are multiplied by a leverage 
of 2 to 3 for antenna temperature determinations in the 
region of 0 to 100 K. Figure 2 illustrates this condition. How- 
ever, this leverage degradation could be overcome if the cal- 
load could be calibrated using the cosmic background tempera- 
ture of 2.7 K as a calibration source via a technique known as 
‘tipping curves.” This technique does not require accurate 
knowledge of insertion losses between the load and receiver. 

A tipping curve is done by pointing the radiometer antenna 
at elevation angles nominally corresponding to one, two, 
three, and four air masses (elevation angles of 90, 30, 19.47, 
and 14.47 deg respectively) and noting the radiometer output 
in “counts.” This information is used to extrapolate to the 
value of counts at zero air masses. Since this condition of zero 
air masses conesponds to viewing only the cosmic background 
radiation of 2.7 K, the number of counts at 2.7 K can be 
derived. The relationship between counts at zero, one, and two 
air masses is, as an cxai^.plc (to a good approximation), = 
2N^ - N 2 ~ the counts corresponding to the cosmic back- 
ground radiometric noise temperature of 2.7 K (7^). The 
radiometer is then switched to t!ie ambient cal4oad and the 
counts are noted. From these measurements, the radiom- 
eter gain in counts/kelvin is found to be 
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Sky noise temperature can then be determined from 




( 3 ) 


The third cooling device was a multistage, closed-cycle 
fluorocarbon refrigerator that achieved -120’’C (153 K) and 
showed great promise. The fundamental problem with this 
machine was that it was intended for a laboratory environ- 
ment, and after several fleld tesu it failed. The manufacturer 
did major redesign and modification for fleld use, but before 
it could be tried again, the decision was made to abandon the 
cooled caldoad and proceed with the development of a heated 
reference load. This decision was based on a critical examina- 
tion of the complexity and potential unreliability of available 
cryogenic refrigerators when used under fleld conditions. 


where Ng is the counts when the radiometer views the sky. 

It should be noted that these equations ignore several small 
correction terms to avoid complicating the discussion. In 
actual practice, various c'^rrections related to noise tempera- 
ture saturation and instrument parameters are applied. The 
sky temperature can be determined as shown above since the 
system has been designed to be linear. A good feature of this 
tipping curve calibration is that the cosmic background radia- 
tion, used as an equivalent cold calibration load, is actually 
lower than the sl^ temperatures the instrument normally 
measures. When used in conjunction with the ambient load. 
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as is done here, the resulting calibration curve allows sky 
measurements to be interpolated rather than extrapolated. 
This means that sky measurement errors, due to cal-Joad 
uncert.!unties, are less than these cal-load uncertainties them- 
sdves. 

There are some drawbacks to this method. First, the atmo- 
sphere is not always well-behaved enough to allow tipping 
curve calibrations of instrument gain. The possibility of some 
“lumpiness” in water content will prevent a correct extra- 
polation of the tipping curve data to zero air masses. To assist 
in dealing with this problem, a four elevation point tipping 
curve is used to verify proper atmosphere stratification (atmo- 
spheric parameters constant with respect to horizontal move- 
ment). The second drawback of this method is that accom- 
plishing a good tipping curve calibration is relatively time- 
consuming. 

I 

f 

f Since the radiometer may not hold its gain calibration for 

the long periods of time between tipping curve tests, a second 
calibration method is carried out by using the ambient load in 
c conjunction with some nonambient load. In the case of the 

!' present WVR, high reliability could not be achieved with 
' coolers for a cold cal-load, and it was decided to use a hot 

cal-load as the nonambient load. 

It is possible to solve for the radiometnc noise temperature 
of th hot load by using the tipping curve derived gain for the 
ladiometer and switching to the hot cal-load as an additional 
real-time step in the tipping curve sequence. Then the hot load 
radiometric noise temperature can be calibrated with relation- 
ship 

^HR ” ^ ^AP 

1 

In actual practice, the radiometric hot load temperature will 
be somewhat less than (but strictly related to) its physical 
temperature. Since the stability of the hot cal-load, a passive 
j device, is many times better than that of the radiometer, an 
active device with over 100 dB system gain, the calibrated hot 
load can be used in radiometer measurements between infre- 
quent tipping curve cahbrations. 

I This is done by defining a new gain based on the now- 

I calibrated hot load. 


I 
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where the subscript 2 refers to the latest real-time dctermina- f 
tior of counts and radiometric temperatures. I 

f. 

Ideally, in a drift-free instrument, In operation, ; 

the instrument measures counts on sky, hot load, and ambient , 

load, and physical temperatures on hot load and ambient load. ] 

Equation (5) is evaluated each time a sky temperature mea- 
surement is made. The resulting new K 2 (determined every 
5 minutes, or so) is used in real time to determine radio- 
metric sky temperatures: 


2 


The value of ^ ^ continuously compared to (the 
tipping curve deri^'ed gain) and its previous values as a monitor 
of radiometer performance. 


III. Frequency Selection 

The WVR frequency pair of 18.55 and 22.235 GHz were 
originally chosen only on the basis of sensitivity to separation 
of the noise temperature effects of water vapor from water 
droplets (clouds). Frequencies near the resonance line (desig- 
nated probing frequencies) are more sensitive to vapor than 
frequencies well off resonance (designated reference 
frequencies). 

Possible pairs were 26.5/22.2, 18.0/22.2, and 31.4/22.2 
GHz, for reference and probing frequencies respectively, in 
the order of increasing accuracy to perform the measurement. 
Since the 18.0/22.2-GHz frequency pair would achieve almost 
the same accuracy as the 31.4/22.2-GHz pair, and could be 
accomplished in a single waveguide-band system, it was chosen 
for this radiometer. 

More recent advanced models have shown that the relation- 
ship between noise temperature and phase delay, both induced 
by atmospheric water vapor, is excessively altitude-profile 
dependent with the 18.0/22.2-GHz frequency pair. This causes 
an additional tropospheric delay error source due to modeling 
error. However, several other frequency pairs exhibit much 
lower vapor profile dependency (Ref. 7). The two best fre- 
quency pairs are 20.0/26.5 GHz and 20.3/31.4 GHz. Analysis 
predicts that both of these pairs will have an error in measur- 
ing water vapor induced phase delay of less than 2 cm rms 
down to 15 deg elevation (^4 air masses) under nonprecipi- 
tating weather conditions. The relative accuracy of the two 
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pain favon the 203/31.4 pair slightly, about 15 percent 
(1.5 mm) at zenith, decreasing to equality at 15 deg elevation. 

Since the performance figures mentioned are based on an 
analytical model, it would seem prudent to experimentally 
verify them. A N£MS (Nimbus £ Microwave Spectrometer) 
has been refurbished and retuned to 21.0/31.4 GHz. (This 
radiometer would not operate at a frequency less than 21.0 
GHz without serious performance degradation.) NEMS will 
be used to verify the 203/3 1.4-GHz pair. This modification is 
now complete and tests are awaiting completion of software 
uevelopment for the radiometer microcontroller. 

To verify the one-band frequency pair, 20.0/26.5 GHz, it 
was proposed to modify the DSN Developmental WVR pre- 
sently at IS. 55/22.235 GHz. This modification consists of 
changing the local oscillators, mixer, isolators, and horn 
antenna Hie work has been completed except for the 
antenna. Since there exists a large amount of operational 
history with this WVR at the 18.55/22.235-GHz frequency 
pair, the option of restoring this frequency capability has 
been included in the modification package. The installa- 
tion of the modification kit that changes between the 18.55/ 
22.235-GHz pair and the 20.0/26.5-GHz pair requires about 
one day. In fact, this option has been exercised once already 
in support of a VXBl experiment scheduled at DSS 13 during 
August 1978. 

One observation regarding the relative desirability of the 
20.3/31.4- and 20.0/26.5-GHz pair is offered here. A 20.3/ 
31.4-GHz WVR requires two complete and independent 
microwave receivers, their associated antennas and calibration 
loads, while a 20.0/26.5-GHz WVR requires only a single 
microwave receiver and associated parts. A life-cycle cost 
analysis (Ref. 8) shows a 25 to 30 percent advantage in favor 
of the 20.0/26.5-GHz pair when examining the microwave 
portion of the WVR ^stem only, simply due to fewer compo- 
nents, and only a small decrease in performance (^1 .5 mm rms 
additional delay error). These accuracy figures are modeled 
and have not yet been experimentally verified. 


IV. WVR Calibration and Field Tests 

During August and December 1977, a series of WVR 
antenna temperature measurements were made at Edwards 
Air Force Base, about 50 miles north of Los Angeles in the 
Mojave Desert. Radiosondes are launched daily there; and 
simultaneous WVR microwave measurements were made. 
Radiosondes measure temperature, pressure, and relative 
humidity in flight and these can be used to determine the 
tropospheric range delay effect of tire atmosphere along the 
path of the radiosonde. 

Comparison of microwave and atmospheric measurements 
results in a WVR calibration linking antenna noise tempera- 
ture and tropospheric delay. Tipping curves were also 
carried out to determine the radiometric temperature of the 
hot calibration load. 

Data from the Edwards tests are still being examined. 
Preliminary studies indicate that the experimental measure- 
ments encompass primarily clear atmospheres, with water 
vapor only and no liquid water (clouds) present. This cannot 
be considered a sufficient data base from which to calibrate 
the WVR completely; however, certain theoretical considera- 
tions may be invoked to generalize the WVR calibration. A 
complete calibration analysis and description of VLBl support 
tests will appear in a later volume. 

V. Future Work 

The DSN Developmental WVR will be used in conjunction 
with Project ARIES (Astronomical Radio Interferometric 
Earth Surveying) (Ref. 9). Additional work should be done to 
experimentally verify the analytical model used in selecting 
optimum frequencies. This could be done by completing the 
frequency modifications to both the NEMS radiometer and 
DSN Developmental WVR, followed by their use in VLBI 
measurements and direct comparison with a LASER/micro- 
wave ranging instrument (Ref. 10). Additional radiosonde/ 
WVR comparisons should be made to calibrate the WVR over 
a wide range of atmospheric conditions. 
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Appendix 


The hot cal-load used in the DSN water vapor radiometer is a JPL-fabricated wave- 
guide unit. The waveguide housing is a S.08-cm-diameter X 14-cm-long copper bar, into 
which a WR42 waveguide size hole is machined. Die load element is a 10.2-cm-long, 
compound single-taper piece comprised of 60 percent beryllium oxide homogeneously 
mixed with 40 percent silicon carbide. This material mixture combines good thermal 
conductivity (2.2 W/cm “C compared with 3.9 W/cm *C for the copper housii^) and 
excellent microwave qualities (36 dB/cm loss). Reflection coefficient of the completed 
load over the frequencies of intere.it is better than 0.008. The heater is wound on the 
copper housing insulated by a C.08-mm mica sheet. This is then surrounded by a 4-cm 
layer of molded ceramic fiber insulation and assembled into a cylindrical aluminum hous- 
ing to provide mechanical support and mounting capability. Heat leakage from the 
assembly is 40 W at 0**C ambient and ISO^C load temperature, well within the heating 
capability of 100 W. A photograph of the microwave components of this load is shown 
in Fig. 3. This load was installed in June 1977, and has been used for all tests since then. 




